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ABSTRACT 

This report describes a user oriented computer program for analyzing micro- 
wave and millimeter-wave mixers. The program can be used for mixers with a single 
Schottky-barrier diode whose I-V and C-V characteristics are known. The diode 
mount is assumed lossless but may have external loads at any number of sideband 
and LO harmonic frequencies. 

The program first performs a nonlinear analysis to determine the diode conduct- 
ance and capacitance waveforms produced by the local oscillator. A small-signal linear 
analysis is then used to find the conversion loss, port impedances, and input noise tem- 
perature of the mixer, hi the noise analysis thermal noise from the series resistance 
of the diode and shot noise from the periodically pumped current in the diode conductance 
are considered. 

Examples are given, including a study of the effects of the series inductance and 
diode capacitance on the performance of some simple mixer circuits using (i) a conven- 
tional Schottky diode, (ii) a Schottky diode in which there is no capacitance variation, 
and (iii) a Mott diode. It is shown that the parametric effects of the voltage depen- 
dent capacitance of a conventional Schottky diode may be either detrimental or ben- 
eficial depending on the diode and circuit parameters. 
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1. Introduction 


As a result of recent work on the theory of microwave and millimeter-wave 
mixers [ 1, 2] , it is now possible to predict their performance with reasonable ac- 
curacy. The purpose of this report is to make available a user oriented computer 
program for determining the performance of single-diode mixers. The relevant mixer 
theory is reviewed and the implementation of the Fortran code described. To illustrate 
the operation of the program a study of the performance of three simple mixers is 
presented. These are: an ideal Schottky-diode mixer in which the junction capacitance 
has the usual voltage dependence, a Schottky-diode mixer in which the junction capaci- 
tance is assumed independent of voltage, and a Mott-diode mixer for which the junction 
capacitance is an experimentally determined function of applied voltage. 

The mixer analysis is carried out in three parts. First, the diode conductance 
and capacitance waveforms produced by the local oscillator are determined using a non- 
linear circuit analysis. Next, a linear small-signal analysis is used to find the input 
and output impedances and the conversion loss between the mixer ports. Finally, the 
shot and thermal noise contributions from the diode are determined. 

The most difficult step in analyzing a mixer is finding the diode waveforms pro- 
duced by the local oscillator. In earlier work by Torrey and Wliitmer [3 ] and others, 
a sinusoidal driving voltage across the diode junction was assumed, with the implication 
that the harmonics of the local oscillator were short circuited. Fieri and Cohen [4] 
used a numerical integration algorithm to obtain the voltage waveforms for a diode with 
nonlinear conductance and capacitance in a simple lumped element embedding network. 
However this method can not be used if the embedding circuit contains distributed ele- 
ments, which is visually the case for microwave mixers. Egami [5] , using a harmonic 
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balance technique, was able to consider arbitrary embedding impedances, but conver- 
gence was difficult to achieve when more than three local oscillator harmonics were 
considered. In another approach, developed by Gwarek [ 6] , the embedding network 
is represented as a simple lumped element circuit in series with a string of voltage 
sources one at each harmonic of the local oscillator. The amplitudes and phases of 
these generators are adjusted to keep the apparent terminal impedance of the circuit n 

equal to that of the actual embedding network. Although the scheme works well for 
many mixers, it is strongly dependent on the guessed values of the lumped circuit ele- 
ments and does not converge for all embedding impedances. A final technique, based on 
an earlier one developed in our laboratory [ 7 ] , solves the nonlinear problem as a 
series of reflections between the diode and embedding network. The algorithm operates 
in the time domain when considering the diode and in the frequency domain when dealing 
with the embedding network. Although this multiple reflection technique often requires 
more time than Gwarek's method, solutions have been obtained for all the mixer circuits 
we have tried and no initial guesses are required. As a result tins method is used in the 
large-signal section of the mixer analysis program described in this report. 

The small-signal analysis follows that of Held and Kerr [ l] which is an exten- 
sion of the original theory of frequency conversion put forward by Torrey and Whitmer 
[3 ] . The small-signal properties of the mixer are derived from a knowledge of the 
large-signal waveforms at the diode terminals and the impedance of the embedding 
network. 

In analyzing the noise properties of a mixer it is important to take into account 3 

the partial correlation of the periodically varying shot noise at the various sideband fre- 
quencies. Following earlier work by Strutt [8], van der Ziel [9], Kim [10], 
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Dragone [11] , and Uhlir* [12] , the theory of noise in Bchottky diode mixers was inves- 
tigated hy Held and Kerr [1] and put into sl form suitable for computer analysis. This 
noise theory is used in the computer program presented in this report. 

The mixer analysis program allows arbitrary embedding impedances at the 
harmonics of the local oscillator and the sideband frequencies, and any Sehottky diode 
can be used. The diode series resistance includes a frequency dependent component 
due to skin effect. The program cannot handle diodes with voltage dependent series 
resistance or in which charge storage effects are significant, i. e. it is assumed that 
the carrier recombination time is small compared with the period*of the local oscillator. 
In all cases the diode mount is assumed lossless and reciprocal. The run time for the 
program as listed in Appendix I is 32 seconds on an IBM 360/95 computer. 


* 


The authors thank Dr. Uhlir for bringing his excellent paper to their attention. 
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2. Outline of the Necessary Theory 

The parameters most often used to characterize the small signal performance of 
a mixer are: the conversion loss between the signal and intermediate frequencies, the 

input and output impedances, and the equivalent input noise temperature. These quanti- 
ties can all be determined if the Fourier coefficients of the diode conductance and capa- 
citance waveforms and the embedding impedance at each sideband frequency are known. 
The first step in the analysis is therefore to determine the conductance and capacitance 
waveforms produced at the diode by the local oscillator. 

2. 1 Solution of the Large-Signal Problem 

The equivalent circuit of the mixer is shown in Fig. 1. Expressions for the 



Fig. 1: The equivalent circuit of the mixer. 
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diode conductance current and capacitance current are well known? 


i g. " 1 s Cexp(av d ) - 1] , 


( 1 ) 


■ - fid 

S " Cj dt * 


( 2 ) 


where 


a = q/rjkT , 


( 3 ) 


and 


c j =c o ( 1 ' v d /0 )“ Y 


(4) 


The differential conductance is obtained from (1): 


di. 


% - ar = ai s = “<V + V - “V 

j J 


(5) 


The voltage and current at the diode can be written in terms of their Fourier coefficients 
as follows: 


-f* CO 


T d<*> 


V V 3 XP (J „» t) , V = < 
n *■ -n n 


L, 

Xi = -oo 


( 6 ) 
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Ut) = i_ (t) + i (t) = ) I d exp(jnu) t) , I =1, (7) 

d S J 5 n t-„ * P * 

where 0 )^ is the local oscillator frequency. 

The conditions imposed on the steady state response of the mixer circuit by the 
embedding network can more conveniently be expressed in the frequency domain. Re- 
ferring to Fig. 1: 

+ R (niu )] , n = 2, 3, ... “> (8a) 

s p 

= - [> e (w p ) + R s (u> p )] , (8b) 

= - [Z e (0) + R s (0)] , (8c) 

where and V Q are the Thevenin equivalent LO and dc voltages seen by the diode. 

The frequency dependence of R is due to the skin effect which we will assume to be 

s 

represented by 


v . 


n 


n 




V d 1 - V 1 


Vj - V„ 
d o 0 
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B S W =B s (d0) +R skin ft • 


( 9 ) 


where ( OhmsA/Hz ) is a constant for a particular diode. The effects of diode 

heating on the series resistance* will not be included in this report. 


If (6) and (7) can be solved for Vj(t) and i^(t) then {l)-(5) can be used to find 

i (t), g.(t) and c.(t) which along with Z (u>) determine the small signal behavior of 
gj J J e 

the mixer. As mentioned in the introduction, a method of solution which works well for 
a broad range of embedding impedances is the multiple reflection technique [7] which 
is described below. 

The circuit of Fig. 1 is modified by inserting a transmission line of arbitrary 
characteristic impedance Z Q between the diode and the embedding network as shown in 
Fig. 2. By making the transmission line an integral number of wavelengths long at the 
LO frequency, and hence at the harmonics of the LO, the steady-state waveforms in the 
modified network will be the same as those in the original circuit. The desired voltages 
and currents can then be obtained by alternately solving the two circuits shown in Fig. 3 
in an iterative procedure. 


o 
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Decker and Weinreb {see ref. [1]) have shown that the dc series resistance of the q 

diode determined from a dc log I-V plot includes a negative component due to heating 
of the diode by the test signal. When the resistance is determined at microwave fre- 
quencies this effect is no longer present and thus the deduced value of de resistance 
will be larger than the apparent value obtained from the dc log I-V plot. 


I. 1 

embedding! 
* network ”*1 


aioae_^ 
'series 
resistance 


(.transmission . 
line 


intrinsic 

diode 



0 


0 


Fig. 2: The equivalent circuit of the mixer modified by the insertion of a long trans- 

mission line of length m^ (m -* <») and characteristic impedance Z Q . The 

steady state v^t) and i^(t) of this circuit will be the same as in the circuit of 
Fig. 1. V Q and Vj are the Thevenin equivalent dc and LO source voltages. 
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Fig. 3.’ The two circuits which must be solved iteratively to obtain the steady state 
solution to the network of Fig. 2. The linear network (a) is analyzed in the 
frequency domain while the nonlinear circuit (b) is solved in the time domain. 
V Q and V are the Thevenin equivalent dc and LO source voltages. 
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Initially a wave propagates along the transmission line towards the diode in 
Pig. 2. Tins wave contains voltage components only at the local oscillator frequency 
and do. The initial do wave incident on the diode is 


% = V 0 V (Z 0 + Z e<°> +R s< 0 >)’ 


and the initial incident wave at frequency tu^ is 


\ = V 1 V < Z 0 + Z e« + E W) ) • 


whore Vq and are the Thevenin equivalent dc and LO source voltages seen by the 
diode and Z (n) represents the embedding network impedance at harmonic n (frequency 
nil) ) of the local oscillator. 

Jr 

After the steady state has been reached between the transmission line and the 
diode, the wave reflected from the diode contains components at all harmonies of the 
local oscillator. Denoting the incident and reflected voltages at frequency n U) by 

ir 

V. and V we have at the diode: 
i 

n n 


V d -\ + V r > 

nun 


\ - ( vv*« 


r? 
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V, and are the Fourier coefficients of the diode voltage and current as given 

in (6) and (7), They can he determined by solving the state equation of the network in 
Fig. 3b in the time domain and then performing a Fourier analysis. Inspection of 
Fig. 3b yields the state equation: 

'If 



Note that v (t) = + V, on the first cycle, 

sol 

Equation (14) can be solved for v^(t), and then i^t), , and 1^ can be 

n n 

found using (6) and (7), Now (12) and (13) can be used to determine : 

n 





(15) 


or 





Z 0 )/2 . 


(13) 


The wave represented by V travels back towards the embedding network where it is 

r n 

reflected and a new wave is launched towards the diode. The amplitude of this newly 
reflected wave is determined by the embedding network reflection coefficient at each 


v 


P 


LO harmonic: 
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Z e (n) + R s( n) -z 0 
p n " Z e (n) + R s (n) + z 0 


(17) 


and is 





(18) 


The wave now incident on the diode is the sum of the reflected wave given in (18) and 

that of the previous reflection cycle, i. e. V. -» V. + p V . At the diode a 

x n *n r n 

new state equation can be written and solved to obtain the voltage and current in the time 
domain. The Fourier coefficients can be calculated and the wave reflected from the diode 
found from (15) and (16). The n<°w incident wave is again determined using (18) and the 
cycle is repeated until the voltages and currents at the two ends of the transmission line 
are equal, that is (referring to Fig. 2): 


and 




V, = V ( for n > 1 
d_ e„ 



I 


n 


(19) 


for n > 1 . 


( 20 ) 
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At this point the solution has completely converged and 





I Z (n) + R_(n) | , for n > 1 . 

‘ 6 S 


(21) 


In practice perfect convergence will not be achieved in a finite time. However, 
a convergence parameter can be formed to indicate when the waveforms at the two ends 
of the transmission line are sufficiently close to one another to merit halting the itera- 
tion procedure: 


6(n) = { I (V d / I, ) | / | Z e (n) + p <n) | } - 1 , for n >1 . 
n n ' 


( 22 ) 


The solution is considered to have converged when 6(n) ^ e for all n > 1, where e 

is an arbitrarily set limit. 


2.2 Small-Signal Analysis 


2.2. 1 Frequency and Subscript Notation 


If a mixer is pumped at frequency uj and has an intermediate frequency uj , 

p u 

the only small-signals which can pi'oduce an IF response are at the sideband frequencies 
(«j 0 + nujp, n=0, ±1, ±2 ...). Following Saleh [13] it is useful to define the side- 
band frequencies by: 


uj = ui n + nuu 
n 0 p 


n = 0, ±1, ±2, ±3, . . . 


(23) 
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For n < 1 the sideband frequencies are seen to be negative, A brief comment on the 
meaning of these negative frequency terms is given in the footnote. * 

Saleh's frequency notation leads to a considerable simplification of the mixer 
theory. Using this notation all upper sideband frequencies (u) Q + j n | tu ) are con- 
sidered positive, while all lower sideband frequencies (u> n - jn| oj ) are negative. 

The sideband frequency index n is used as a subscript with the various electrical quan- 
tities and hence the upper sideband, intermediate, and lower sideband frequencies are: 
u> + i ' « tu 0 + ujp, u) Q , and « (Uq - oi^; and V +1 , V Q and V_ 1 represent 
voltages at these frequencies. 

2.2.2 Small-Signal Conversion Matrix 

Using tlie sideband notation described in the previous section let §1 and §V 
denote the vectors of the small-signal sideband currents (6l n ) and voltages (6V n ) at the 
terminals of the intrinsic diode (the diode excluding its series resistance). Then 

2 = [..., 6I r 6I 0 , 6I_ r ... ]* 


Electrical quantities are frequently described by a single complex quantity associated 
with some frequency, assumed positive. For example, a voltage of frequency uj may 
be described simply by its complex half-amplitude V, implying an instantaneous volt- 
age v(t) = V + V* ^ . It is just as meaningful to work with a negative fre- 

quency (- tu) and the conjugate of the complex half-amplitude (V*), provided the conven- 
tion is clearly understood. Impedances and admittances are then simply the conjugates 
of their conventional positive frequency values, i. e. Z.(-io) = V*/I* = Z *(u>). 
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and 


6V = flVj, 6V 0 , 6V_ X , ...] 


(25) 


r^j r^j 


Torrey and wnitmer [3] have shown that 6 1 and 6V are related via a conversion ad- 
mittance matrix ¥ defined by 


61 = Y 6V 


(26) 


If the row and column numbering of the square matrix Y correspond with the sideband 
numbering, Y can be written out as 


Y = 


Y 

Y 


• 

• 

• 

# 

4 

9 

a 

* 

a 


11 

Y 10 

Y l-1 

a ■ a 

01 

Y oo 

Y o-i 

a a a 

-11 

Y -10 

Y -l-l 

a a • 


(27) 


with element values given by [3 ] : 


Y = G + j(ui n + mu> ) C 
mn m-n 0 p m-n 


( 28 ) 


G m „ and C are the (m-n) th Fourier coefficients of the diode conductance and 
m-n V&-11 

capacitance waveforms, g,(t) and c.(t), and are defined by: 

J J 


CO 

gj(t) = Y, G n = G n » (29) 

n = “Co 


•j» 


I C n op [Jump t] , c_ n = C* 
n = -» 


m-n 


and C „ are given by: 
m-n 


(30) 
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m-n 
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m-n 


J n T/2 

T J 

T 

2 


u T/2 

_ T 
2 


gj(t) exp [ -j (m-n) u> p t] dt , 


e.(t) exp [ - j (m-n) u> t] dt , 


(31) 


(32) 


where the integration is taken over one period T of the local oscillator. 

The matrix Y can be regarded as the admittance matrix of a multifrequency 
multiport network, as shown in Fig. 4, in which there is one port for every sideband 
frequency uj . If the embedding impedances 2 and diode series resistance R 

JlL 6 S 

n n 

corresponding to the sideband frequencies tu n are now connected in parallel with the 

intrinsic diode an augmented network is formed as shown by the broken line in Fig. 4. 


IQ 


Sv >! zt 


Ugimoip 


SV t JP 


intrinsic diode 


augmented network 


The small signal representation of the mixer as a multifrequency linear multi- 
port network. The voltage and current fiV m and 6l m are the small-signal com- 
ponents at frequency (Wq + m Uj^) at the intrinsic diode. Each port represents 

one sideband frequency. The conversion matrix Y is the admittance matrix of 
the intrinsic diode. The augmented network includes all the sideband embedding 


impedances Z 


and is characterized by the augmented admittance matrix 


Y' . 61^ is the equivalent signal current generator which is connected at port 1 
during normal mixer operation, the other ports being open circuited. During the 
noise analysis equivalent shot and thermal noise current sources 6Ig and filih, 

m i 

are connected to all ports. 
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The ports of the augmented network correspond to the terminals of the intrinsic diode 
at the various sideband frequencies and do not represent physically accessible ports in 
the real mixer. The augmented network can be described by the admittance matrix Y f , 
defined by 


ar = y* 6V 

where 

61* 61^ , 6Iq j 5ll r ... ]* 


(33) 


and 


6V = [..., 6V 1( 6V q , 6V_ r ...]* 


(34) 


6V and 61* are the small signal voltage and current, at sideband U) = uj_ + 
mm m 0 

m uu (port m), of the augmented network. The elements of the augmented ad- 
P 

mittance matrix Y* are given by 


Y T = Y 
mn mn 


m ^ n 


(35a) 


and 


Y* = Y 
mm mm 


IX 


+ ] 


-l 


m = n 


m 


m 


(35b) 
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Inverting (33) gives 


6V = Z' 6 1' , 


where 


Z' = (Y r ) A . 


The impedance matrix Z' enables us to calculate the conversion loss and the input and 
output impedances of the mixer and is also needed in computing the noise temperature. 


2. 2.3 Mixer Port Impedances 


The impedance of any port of the intrinsic diode (see Pig. 4) can be found 

by open cii-euiting the corresponding embedding impedance Z and then forming the 

m 

Z' matrix defined by (36). The desired port impedance is given by the mm-th element 
of the newly formed Z 1 matrix: 


Z =5 Z’ 
m mm, 05 , 


where the subscript “ indicates that Z' has been formed with Z open circuited. 

e m 

The corresponding mixer input impedance seen by the embedding circuit includes the 

A- 

diode series resistance and is 


+ JR = Z T m + R 
s mm, » s 
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In particular tlie IF output impedance is given by 


J IF 88 2 in rt 
out 0 


= Z n + E_ 


Z 00, » + R s t 


(40) 


Throughout the remainder of this report it will be assumed that the IF load im- 

pedance Z . is conjugate-matched to the IF output impedance of the mixer, thereby 
e 0 

minimizing the conversion loss. Once the mixer performance with matched IF is 
known it is a simple matter to calculate the performance with any other IF termination. 
The value of the conjugate-matched IF load impedance is, using (40), 


= Z 


IF 


out 




00 , 


+ R s ) 

s o 


Z 00, « + R s 0 * 


(41) 


where Zl n is the center element of the Z' matrix with Z = <=> . Eather than 
00, « e Q 

reforming the Z' matrix each time an input impedance is calculated, the intrinsic 

diode port impedance Z m can be found from (referring to Fig. 4): 


Z' 

mm 



+ E 


m 



(42) 


where Z' is the mm-th element of the mixer impedance matrix formed with the IF 
mm 

load impedance conjugate-matched to the IF output impedance. The corresponding 
mixer input impedance is then 
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B 


(Z -f R ) Z' 

' e s ; mm 

m m 


(Z + R ) - Z' 
m e m s m mm 


(43) 


2. 2. 4 Conversion Loss 


The conversion loss from sideband j to sideband i in a mixer is 


power available from the signal source at sideband ui. 

J ij i ^ converted power from the signal source, delivered to 

the load at sideband cu . 

x 


L.. A 


) • 


Consider for the moment only the intrinsic diode shown in Fig. 4. The power available 
from impedance (Z + R ) at the j-th sideband is 

6 a 3 a 

3 J 


Available = * l 5l jl Be[Z^ + R. 3 . 


(44) 


The power delivered to a load impedance (Z + R ) at sideband i is, using (36) 

6* S a 

1 1 


l z ij | 2 leijl* Belize +R Sj ] 

IV + E ». I 2 


delivered 


» Be [6V. 6I£ ] - 


(45) 
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Dividing (44) by (45) gives the conversion loss L!. of the intrinsic diode: 


available _ , 


Z + R 
I e. s. 


|z 


delivered 


ij 


+ R_ 

fj fi 


4 ! Z ij I +E Sj ] Se C Z e. +H J 


(46) 


To find the conversion loss of the actual mixer, the additional loss in the series 
resistance at the input and output frequencies must be included. This is accounted for 
by two terms: 


K. 

J 


A 


fQ 

Power Available from Z 

e. 

1 

Power available from Z + R 

e* s. 

J J 


Re[Z e + R g ]/Re[Z e ] (47) 

j j 3 


K. 

x 


A 


Power delivered to Z + R„ 

e. s 

Power delivered to Z„ 

e. 

i 


Re[Z e + R ]/Re[Z ] . (48) 

i i i 


Multiplying K. and Eh by Ljj gives the expression for the loss of the mixer in convert- 
ing from sideband j to sideband i: 



(49) 


where Z!. is the ij** 1 
il 

defined in (36), 


element of the impedance matrix Z' 


of the augmented network as 


2.3 Mixer Noise Theory 


{ 

The noise observed in a Schottky diode comes mainly from three sources, 

(i) shot noise due to the statistical nature of the current flow across the depletion 
layer, (ii) thermal noise due to the random motion of the charge carriers in the un- 
depleted semiconductor material and (iii) noise due to lattice scattering, which occurs 
in regions of high current density. At room temperature the noise contribution due to 
lattice scattering is usually small enough to be approximated by a slight increase in the 
temperature of the diode series resistance to T [ 1, 2] . In cryogenic mixers how- 
ever, scattering noise may account for a more substantial part of the overall noise and 
a more complex analysis than is performed here is required to take account of its par- 
tially correlated components. 

The equivalent circuit of the Schottky diode, including noise sources, is shown 

in Fig. 5(a). T g ^ is the equivalent temperature of the series resistance and includes 

the effects of lattice scattering and pump heating, k is Boltzmann's constant and q the 

electronic charge. In Fig. 5(b) the thermal and shot noise are both represented as 

9 2 

equivalent current sources in parallel with the intrinsic diode. fii T and 6ig are the 
mean-square values of the thermal and shot noise currents in the frequency range f to 
f+ Af. These current sources can be regarded as generating a multitude of quasi- 
sinusoidal frequency components, each with its own amplitude and phase. In the mul- 
tifrequency multiport equivalent circuit of the mixer, Fig. 4, the noise sources can be 
included by connecting a noise current source at each sideband frequency to the appropri- 
ate port of the augmented network. 
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(a) 



Fig. 5: (a) The equivalent cii'cuit of the diode including noise sources. T is 

the equivalent temperature oi the series resistance and takes into account 
pump heating and lattice scattering noise, (b) The equivalent circuit with 

.2 *2 

the thermal noise source fiv^, transformed into a current source . The 
sideband components ol’ the noise sources are treated in the same way as 
the equivalent small-signal current sources 61^ applied at the ports of the 
augmented network of Fig. 4. 


i 
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As seen by the intrinsic diode, the thermal noise current source in Fig. 5(b) is 



4 k T 


eq 



Af 



(50) 


At the intermediate frequency (m = 0) the mixer is conjugate matched with Z - 

e 0 

(Z n + R )* - As seen by the IF load Z„ the equivalent thermal noise current 
0 s 0 e 0 

source connected to the IF port of the augmented network is 



4k T 


'eq 




Af 


4 k T eq 





Af 

2~ 


(51) 


Since the noise performance of the mixer is determined from the total noise power de- 
livered to the IF load Z , (51) rather than (50) must be used to account for IF noise 

0 

in the series resistance that is delivered to the IF load. 

The mean square value of the shot noise current source in Fig. 5(b) is given by 
the usual expression: 



2 q i„ Af • 

s j 


( 52 ) 
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2.3.1 Shot Noise 

The shot noise in a mixer arises from the current produced in the diode con- 
ductance by the local oscillator and dc bias. It can be considered as white (Gaussian) noise, 
amplitude modulated by the LO waveform. Draggne [11] and Uhlir [12] have investigated 
\ the properties of this modulated noise and shown that there is partial correlation between 

the quasi-sinusoidal components at the various sideband frequencies. The correlated 
components at these sidebands are down converted in the diode to the intermediate fre- 
quency where they add vectorially. 

Let 6Ig represent the quasi-sinusoidal component, at frequency (u n , of the 
n 

periodically pumped shot noise current source in Fig. 5(b). Each of the sideband com- 
ponents 5lq is connected to the complete equivalent circuit of the mixer, Fig. 4, at 
n 

~ rs/ 

the appropriate sideband port. We define &Ig and sVg as the vectors of the input 
shot noise currents and voltages at the ports: 


6Ig = [... 5Ig Slg 51* 




and 


4 






Using (36) the output noise voltage at the IF is 


6V S 


0 




( 53 ) 


z o 6I S 


(54) 


where Zjj is the zeroth or center x'ow of the augmented impedance matrix Z’ defined 
in equation (36). It follows that 

6 \ •« V S 0 -?i •(%%)* - ZJ #1 Z'o • < 5S ) 

where t indicates the conjugate transpose of a matrix. 

Talcing the ensemble average* of (55) yields 

< I 6V Sq | 2 > = Z«j < 6Ig 61’g > Z{j f . (56) 

(dig dig ) is the shot noise current correlation matrix and has the general element 

<611, 611* ) . It can be shown that [11], [12] : 

b 



6 ¥ > " 2 <5 Vn 

n 


Af 


(57) 


where I is the (m-n)th Fourier coefficient pf the diode conductance current- is 
m-n 

defined by 


Talcing the ensemble average is equivalent to considering a small but finite bandwidth 
as must be used in any physical measurement. The finite bandwidth contains a multi- 
tude of quasi-sinusoidal noise components with random amplitudes and phases. 


(58) 


V < £ > " I J n ®!PCio“>p t] , I_ n = I* . 

J n = -o 

* I m-n is tIien given kys 

o 

t 

i r 2 

l m-n = t" J t i g _ (t) exp [ “j (m - n) «) p t] dt , (59) 

2 

where the integration is taken over one LO cycle. 

2.3.2 Thermal Noise 

Thermal noise generated in the diode series resistance has components which 

are uncorrelated at the various sideband frequencies. Let 81^, represent the quasi- 

n 

sinusoidal component at sideband frequency a>^ of the thermal noise current source in 

Fig. 5(b), and let be the sideband noise voltage produced by filfp . The noise 

n n 

voltage produced at the IF port of the augmented network (Fig. 4) by the thermal noise 

at all the sidebands can be found using equation (36): 

9 


6V T 0 =2' Mfc 


0 


J 


(60) 
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rw 

where 6IL = [..., filL, * film > film , ...] is the vector of input thermal noise 
T T 1 T 0 T -1 

currents at the sideband ports of Fig. 4 and Z^ is the center row of the augmented 
impedance matrix Z' . From (60) 


6V T q * 6V T q = ^T ‘ <*0 ^)* = *0 bl' T fil£ Z'* . (61) 

Talcing the ensemble average gives 


< i6V T J 2 > = Z{j < Slfp 6Bj>Zj, t . (62) 


r"'wi i-O ,j+ 

The square matrix ( 6IJp 61^. ) is the thermal noise current correlation matrix. 

Since the thermal noise components at the various sideband frequencies are uneorre- 

* 

lated the elements (6IL 6IL ) =0 unless m = n, i. e. the matrix is diagonal, 

m n 

From equations (50) and (51): 
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< 61 * film > 
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(63b) 
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<611, filjT > = 0 , m n 

m ii 


(63c) 


2.3.3 Total Mixer Noise 

The total output noise voltage of the mixer is obtained by combining the ther- 
mal and shot noise components. From (56) and (62): 


< |v No | 2> - zj [<e% 6i' s f > + <6 r T Si^>] z'J . 


(64) 


It follows that the noise power delivered to the matched IF load Z from the mixer 

e 0 

itself is 


P 0 * <I V N | 2 > R^1/\Z * B 


(65) 




The equivalent input noise temperature T^ of the mixer is the temperature to 
which the signal source conductance must be heated to give the same output noise from 
a noiseless but otherwise identical mixer as the actual mixer would produce when its 
signal source conductance was maintained at absolute zero temperature. Thus, 


t m = p o L oi /rkA± ' 


( 66 ) 
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Using (65) and equation (49) for L 01 
to the IF) 


(the conversion loss from the signal frequency 



<|V N | 2 > 

0 


2 



Ke[Z ] 
®1 


(67) 


where is an element of the augmented impedance matrix Z' defined in (36). 

Tj^ is the single-sideband input noise temperature of the mixer, meaning that 
all the noise of the mixer is attributed to a hypothetical source at the signal frequency. 
It is related to the single sideband noise figure by: 


F 


SSB 


= 1 + 


M 

290 6 K 


( 68 ) 


When describing the performance of a mixer whose physical input port is coupled 
to both the signal and image frequencies, it is more convenient to talk in terms of a 
double-sideband noise temperature is the temperature to which the 

signal and image source conductances must be heated to give the same output noise 
from a noiseless but otherwise identical mixer as the actual mixer would produce when 
its signal and image source conductances were maintained at absolute zero temperature. 
For mixers in which the conversion loss from the signal and image frequencies and 

Lq_i are equal 
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t dsb t ssb /2 


(69) 


and when L Q1 ^ 


2.4 Summary of Mixer Theory 

The performance of a mixer can he characterized by its conversion loss and 
equivalent input noise temperature. These quantities depend on the large-signal wave- 
forms at the diode and on the embedding impedances of the mixer at the small-signal 
sideband frequencies. The diode waveforms can be found using the multiple reflection 
technique described in section 2.1. Once a steady state solution has been obtained, the 
Fourier coefficients of the conductance and capacitance waveforms can be extracted and 
used to find the conversion admittance matrix which relates the small-signal sideband 
currents and voltages of the intrinsic diode. An augmented admittance matrix can then 
be formed, which describes the multiport network consisting of the intrinsic diode, the 
diode series resistance, and the sideband embedding impedances. The inverse of this 
matrix is the augmented impedance matrix Z r , whose elements are used to calculate 
the conversion loss at the various sideband frequencies and the input impedances of the 
mixer ports. It is assumed throughout that the IF load impedance is conjugate-matched 
to the IF output impedance. 


J 0-l 


l ssb 


DSB 


1 + 


•'01 

J 0-l 


(70) 


In the noise analysis, two components are considered: shot noise in the junc- 
tion, and thermal noise in the series resistance. These are represented by equivalent 
noise current sources in parallel with the intrinsic diode. The periodically varying 
shot noise has correlated components, while the thermal noise does not. Correlation 
matrices are formed and evaluated for both shot and thermal noise sources. The shot 
noise correlation matrix has elements related to the Fourier coefficients of the large - 
signal diode conductance current, while the thermal noise correlation matrix depends 
upon the embedding impedances at the sideband frequencies. The two matrices to- 
gether yield the total output noise voltage from which the equivalent input noise tem- 
perature of the mixer can be calculated. The analysis is complete at this point. 

2. 4. 1 Comment on the Effect of the Number of Frequencies Considered 


In transforming these procedures into a workable computer program there is a 
practical limit on the number of harmonics of the local oscillator which can be used. 

This means that the small-signal admittance matrix Y will be truncated above some 
finite harmonic number, which is equivalent to short circuiting the intrinsic diode at all 
higher sideband frequencies, hi the nonlinear analysis the restriction on the n um ber of 
harmonics is equivalent to terminating the intrinsic diode in an impedance Z^, the char- 
acteristic impedance of the hypothetical transmission line, at all higher frequencies. 

As long as the number of harmonics considered is not too small this approximation is a 
reasonable one. 


In the next section the computer program based on this analysis is described. 


3. Description of the Mixer Analysis Program 


Using the theory given in the previous sections, a user oriented computer pro- 
gram was written which can readily accommodate a variety of mixer problems. 

The program requires as inputs (i) the embedding impedances seen by the diode 
at each harmonic of the local oscillator and at the harmonic sidebands, (ii) the diode 
characteristics, including an arbitrary capacitance-voltage dependence, and (iii) the 
operating conditions for the mixer, i, e. the bias voltage applied to the diode, rectified 
current, pump and intermediate frequencies, and the equivalent diode temperature. 
Other valuables which can be input to change specific program operations will be dis- 
cussed later in this section. 

The output includes (i) the large signal current and voltage waveforms at the 
diode, (ii) the Fourier coefficients of the diode conductance and capacitance, (iii) the 
conversion loss between every pair of sideband frequencies, (iv) the IF output imped- 
ance, (v) the input impedance at each sideband, and (vi) the equivalent single-sideband 
input noise temperature. 

The remainder of this section explains the program in detai 7 and illustrates the 
steps used for running it, A complete annotated listing of the program appears in 
Appendix I and a general flowchart is given in Fig. 6. A list of the main program 
variables and their counterparts in the theory developed in parts one and two of this 
report can be found in Fig. 7. 

3.1 Program Implementation 

The program begins with a call to subroutine LGSIG to do the nonlinear analysis 
using the multiple reflection method described in section 2. 1. The embedding network 



Fig. 6: Chart illustrating the general flow of the mixer analysis program and the 

operations of the subroutines. 
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7: Alphabetical list of the main variables used in the theory in this report, and 

their counterparts in the mixer analysis program. 


38 


impedances Z (n) (ZER + jZEI) and the sideband impedances Z (ZEMBSB) are 
6 6 

m 

input in the BLOCK DATA program or formed, up to the highest harmonic of the LO 
(NH), assumed even, in subprogram Z EMBED. The real and imaginary parts of the 
embedding impedance Z g (n) at harmonic n become elements n of the arrays ZER 
and ZEI. The dc term is considered separately. The embedding impedance Z 

e m 

at sideband m becomes, in the notation of section 2. 2. 1, array element 


C 


NH 


+ 1 - m 


^),m = 0, ±1, ±2, ... ± of ZEMBSB. The diode series re- 


NH 

2 


sistance at each LO harmonic (RSLO) and at the sideband frequencies (RSSB) is also 

formed in subprogram Z EMBED using equation (9). The array element notation for 

RSLO is the same as that used for ZER and ZEI, that is R at harmonic n becomes 

s 

array element n of RSLO. Similarly, the notation used for RSSB is the same as that 
of ZEMBSB, R g at sideband m becomes array element + 1 - m^) of RSSB, 


The reflection coefficients, (RHO) of the embedding network at the various LO 
harmonics are found using (17) for a given value of the transmission line characteristic 
impedance Zq(ZO), After initializing the incident voltage at the diode from (10) and 
(11) , the input data are printed through subroutine PRINT1. 

The next step in the nonlinear analysis is to find the local oscillator voltage 
V^(VLO) at which the desired rectified diode current (IDBIAS) is obtained. An outer 
loop (JLO) is established which zeros in on the required LO voltage. To do this, ten 
reflection cycles of the large signal analysis are run after which the dc current 
(IDCOS(l), the first Fourier coefficient of the total current in the diode) is compared 
with the desired value IDBIAS, and VLO is appropriately changed. The loop continues 
until the desired and calculated currents agree within a preset accuracy (IDCACC). 

The time domain calculation of the current and voltage at the diode is accom- 
plished by using the IBM SSP routine DRKGS which solves the network state equation 


39 


(14) for the circuit of Fig, 3b using a Runge-Kutta algorithm. The problem is time 

sealed so that the LO period is 2 tt seconds, Eaoh LO cycle contains a fixed number 

of points (KPTS) and DEKGS is called separately to integrate the cirouit equations 

over each of the (NPTS-1) intervals per LO cycle. Subroutine FCT, called by DRKGS, 

supplies the state equation (14). The resulting values of the total diode ourrent i^ 

(IGJ + ICJ), voltage v^ (Y(l)), oonductanoe g. (GJ), and capacitance c. (CJ) for 

j j 

each point are stored in arrays. The integration loop is repeated NLO times to allow 
a steady state to be reached and then the IBM SSP routine FORIT is used to find the 
Fourier coefficients of the final voltage and current waveforms. At tills point the volt- 
ages reflected from the diode at each LO harmonic can be found from (15) or (1G) and 
the new incident voltages determined from (17) and (18). The convergence of the solu- 
tion is checked using (22) and if more iterations are required (ZQFLAG ^ 0) the pro- 
cedure is repeated from the beginning of the DRKGS integration. However, if the im- 
pedances at both ends of the transmission line are equal within some fixed accuracy e 
(ZQACC) at all harmonics of the LO above the first, then the large signal analysis is 
complete. The results of the last iteration are printed via subroutine PRINT2 and sub- 
routine POWER is called to calculate and print the required LO power, * The Fourier 


We have found that convergence of the nonlinear analysis can generally be made more rapid 
if the embedding impedances at dc and the LO frequency, Z (0) and Z (1), are artiiically 

set equal to the characteristic impedance Zq of the hypothetical transmission line. Provi- 


ded the LO and dc bias sources are appropriately adjusted this has no effect on the steady- 
state waveforms. In computing the LO power required to drive the mixer (i. e. available 
power from the LO source), the correct value of Z e (l) must be replaced and the LO source 

voltage adjusted accordingly. It follows that, in terms of the LO source voltage V^q, when 


Z Q (0) = Z Q , the LO power required to drive the actual mixer (with Z (1) at its original 
value) is: 



> V LO 
SRelZ e (l)] 


R s (l) + Z e (l) + Z Q * ZQ(1) j 2 
Z Q * (1+ ZQ(1))| 2 


(71) 


where ZQ(1) is the convergence parameter given in (22) when n = 1, 
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coefficients of the diode conductance and capacitance are now found using FORIT, after 
having un-tims-scaled the frequency dependent terms. Finally, the diode conductance and 
capacitance, total diode current, and voltage across the diode are plotted over one LO 
cycle using subroutine PLOT. 

After the large-signal analysis is finished subroutine SMSIG is called to perform 

the small -signal analysis described in section 2. 2. The Fourier coefficients of the diode 

conductance and capacitance (GJCOS, GJSIN, CJCOS, CJSIN) obtained in the large-signal 

section are converted to complex form (FG, FC) and printed in subroutine PRINTS. 

Calculation of the conversion loss matrix (XLMAT) and input and output impedances 

then begins by forming the small-signal admittance matrix Y (A in the program) of 

equation (26) using (28). The IF load impedance is open circuited at this stage and the 

augmented admittance matrix of (33) Y' is formed in subroutine YPRIME by adding 

(Z m b H ) to the diagonal terms of Y (A)* Y T (A) is inverted to obtain 
m m 

Z' (A) using the IBM SSP routine MINV (matrix inversion), slightly modified to 
handle a complex matrix (CMINV). The IF output impedance (ZIFOUT) is the center 
element Z^ (A(NHD2P1, NI-ID2P1)) of this matrix plus the diode series resistance R g 

(see eq. (40)). If IF load impedance is set equal to the complex conjugate of ZIFOUT 
and the augmented admittance matrix is again formed and inverted to give Z' (A) for the 
mixer with a matched IF load. The elements of die convei'sion loss matrix (LIJ) and 
the input impedances at the sideband ports (ZIN) are then found using (43) and (49). 


The noise analysis proceeds according to section 2.3 and begins with the forma- 

~ rU J. 

tion of the shot noise correlation matrix ( 6 If 5 If ) in subroutine CORBEL, using 
(57). The equivalent input shot noise temperature (SHOT) referred to the signal port is 
calculated in subroutine TMTX using equation (67) with < 1 6V" N | 2 ) replaced by 


n 

< I 6V„ j ) as given by (56) and (57). The total mixer output noise voltage is found, 

0 

using subroutine TNOISE, by adding the shot noise correlation matrix to the thermal 
noise correlation matrix given by (63) - see equation (64). Subroutine TMIX is then 
used to find the total equivalent input noise temperature (TM) at the signal port. The 
thermal component (THERM) is found by subtracting the shot noise contribution (SHOT) 
from the total noise temperature (TM). The results of the conversion loss and noise 
analyses are printed using subroutine PRINT4. This completes the program. 

For a more detailed description of the operation of the program the reader is 
referred to the comments in the program listing and the general flow chart of Fig. 6. 

3.2 Running the Mixer Analysis Program 

A listing of the mixer analysis program in its simplest form appears in Appen- 
dix I along with the output of a run. Using the IBM Fortran IV-H compiler, the execu- 
tion time for this particular listing is 32 seconds on an IBM 360-95 computer. The 
co mm ents in the listing provide a step by step description of the Fortran coding and ex- 
plain the meanings of the variables. To run the program the following information must 
be supplied by the user via the BLOCK DATA subprogram: 

(1) The embedding network impedances at the LO frequency and the higher harmonies 
as real and imaginary parts (ZER, ZEI) in Ohms. 

(2) The sideband impedances in complex form(ZEMBSB)in Ohms, where sideband 
m corresponds to array element (NH/2 + 1 - m) and there are NH+1 array 
elements in all. Note that for all lower sidebands (m<0) the sideband frequency 
is negative and the ZEMBSB(NH/2 + 1 - m) is therefore the conjugate of the 
positive frequency impedance. 
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(3) The LO frequency (Fp) and the intermediate frequency (IF) in Hertz. 

(4) The dc bias voltage across the diode (VDBIAS) in Volts. 

(5) The desired rectified current (IDBIAS) in Amperes. 

(6) The equivalent temperature of the diode series resistance (TEQ) in degrees Kelvin. 

(7) The diode series resistance at dc (RS) in Ohms. 

(8) The diode skin resistance constant (RSKIN) in Ohms A/Hertz. 

(9) The reverse saturation current of the diode (IS) in Amperes. 

(10) The diode contact potential 0(PBI) in Volts. 

(11) The diode capacitance at zero volts (CO) in Farads. 

(12) The diode capacitance law exponent V (GAM). 

(13) The diode I-V law exponent a £ q/r)kT (ALP) in (Volts)" 1 . 


The values of the remaining variables are somewhat arbitral and have been 
optimized for the example given in Appendix I. The following information may prove 
useful in choosing values for these variables when running other examples. 

The characteristic impedance ( ZO) of the hypothetical transmission line inserted 
between the diode and the embedding network for the nonlinear analysis has a significant 
effect on the number of reflection cycles required for convergence. Values near 200 
Ohms result in a fairly rapid rate of convergence for the example in Appendix I and 
also for the examples in section 4. 1 in which the embedding impedances above the first 
harmonic are open circuited. However a lower value (50 Ohms) gave somewhat faster 
convergence when the embedding impedances above the first harmonic were short 


circuited. 
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Tlxe initial value of the local oscillator voltage (VLO) and the initial increment 
(VLOINC) used to zero in on the desired dc rectified current can be chosen so as to 
avoid many time consuming loops in the large signal analysis. If, as in the examples 
in the next section, many runs are desired with only slight variations in the circuit 
parameters VLO will change very little between each rim and VLOINC should be made 
fairly small. 

The number of LO cycles needed to reach a steady state (NLO) for the circuit 
of Fig. 3b in the nonlinear analysis, need not be greater than one for the examples given 
in this report (bear in mind that the solution will continue settling in successive reflec- 
tion cycles) but some circuits may require additional settling time. NLO can be in- 
creased by changing the appropriate statement in the BLOCK DATA routine. 

If the dc current reached after ten reflection cycles in the nonlinear analysis is 
not sufficiently close to its ultimate value it may be necessary to increase the reflection 
cycle number at which the LO voltage is adjusted. This is accomplished by changing 
the program step which follows statement label 10 in subroutine LGSIG. 

The results of any of the reflection cycles in the nonlinear analysis can be 
printed by adjusting the parameter NPRINT in the BLOCK DATA subprogram which 
causes results to be written out every NPRINT cycles. Upper limits on other program 
loops such as the total number of allowed nonlinear analysis cycles (NITER) or VLO 
adjustments (NVLO) can be increased or decreased as desired by changing the valuables 
appearing in the common block labelled LOOPS. 

The local oscillator cycle was divided into 50 parts (51 points) in these examples 
to yield a reasonable number of data points for plotting the diode current, voltage, 
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capacitance and conductance waveforms and to avoid aliasing. * If the number of points 
(NPTS, assumed odd) is altered some of the array dimension's must also be changed. 

In subroutine LGSIG the variables VDDATA, IDDATA, IGJDAT, ICJDAT, CJDATA and 
GJDATA all have dimension NPTS. 

If other than eight harmonics of the local oscillator are to be considered the 
variable NH (assumed even) must be set to that number in the BLOCK DATA subpro- 
gram. Also the following array dimensions must be changed in subroutines LGSIG, 
SMSIG, FCT, and BLOCK DATA to the value NH, if they represent LO harmonics, or 
to NH + 1, if they refer to the sidebands: 

ZER, ZEI,RSLO, AY, BV, ZEMB, RHO, VI, VR, DVR, ZQMAG and ZQPHA must be 
dimensioned NH and ZEMBSB, RSSB, CJCOS, CJSIN, GJCOS, GJSIN, VDCOS, VDSIN,IDCOS, 
ID8IN, A, COR, FG, FC, T, ZIN, XLMAT, GJMAG, GJPHA, CJMAG, CJPHA, WK1 and WK2 
must be dimensioned NH + 1. In addition certain print formats will also have to be 
altered. 

The relationship chosen to represent the variation of diode series resistance with 
fi equency due to the skin effect involves the determination of a constant (RSKIN) which 
will be different for every diode. An increase in the resistance at 100 GHz. of 30% of 


The sampling theorem indicates that if NH harmonics are considered, it should only be 
necessary to consider 2 * NH + 1 points in the diode waveforms. This would be true if 
the wavewms produced by the Rimge-Kutta integration contained only NH harmonics. 

However the integration solves the circuit of Fig. 3(b) quite faithfully and, because of ° 

the exponential nonlinearity of the diode, harmorirs above NH are present in the wave- 
form. These are ignored in successive reflection cycles of the nonlinear analysis. If 
only 2* NH + 1 points are considered in the waveforms, the phenomenon of aliasing 
will occur, by which higher frequency components are "mixed” with harmonics of the 
sampling frequency thereby causing errors in the computed Fourier coefficients. 
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the dc value has been assumed in the example in Appendix I of this report. RSKIN 
can be changed via the BLOCK DATA subprogram. A frequency dependent reactive 
term may also be added to the diode series resistance by mating RSKIN, RSSB, and 
RSLO complex. 

It is also possible to use a different capacitance-voltage law for the diode than 
that given in equation (4). This is illustrated by the Mott diode examples discussed in 
the next section. 


4. Examples Using the Mixer Analysis Program: A Study of the Effects of the Series 
Inductance and Diode Capacitance on Some Simple Mixer Circuits 

In this section the mixer analysis program is used to study the effects of the 
series inductance and diode capacitance on the performance of the simple mixer circuits 
shown in Pigs. 8 and 9. The embedding networks were chosen to simulate mixers in 
which there is inductance due to the diode package or contact whisker. Higher harmon- 
ics are either short circuited outside the series inductance, or open circuited as shown. 
The mixer performance is investigated for three different diodes: (i) a Sehottky diode 
with Y = 1/2, (ii) the same Sehottky diode but without a voltage-dependent capacitance 
(Y =0), and (iii) a Mott diode with an experimentally determined capacitance- voltage 
relationship {see Pigs. 10 and 11). The program in Appendix I was modified to allow 
multiple runs with varying values of series inductance (LS) or diode capacitance (CO ). 
As these changes are typical of the alterations which might be made by other users of the 
program they are detailedin the next section. 

4. 1 Program Alterations for Running the Examples in this Section 

The main driver program was altered to allow repeated runs with different values 
of diode capacitance CO or series inductance LS and to enable the result's to be stored 
on a disk file. The impedances at the LO harmonics and the sideband frequencies are 
formed at the start of each new run in subroutine Z EMBED rather than Inputting them 
through the BLOCK DATA program,and the diode series resistance is taken to be inde- 
pendent of frequency. (RSKIN=0). The capacitance-voltage relationship of the Mott diode 
(Pig. 11) was determined from the doping profile* (Pig. 10) using a piece-wise linear 

The data on the doping profile of the Mott diode was supplied by M. V. Schneider of 

Bell Telephone Laboratories, Holmdel, N. J. 
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Fig. 11: Piecewise-linear approximation to the capacitance-voltage relationship of the 

Mott diode whose doping profile appears in Fig. 10. The slope and C axis inter- 
cept for each of the linear regions is supplied to the mixer analysis program for 
the calculation of the diode capacitance at any given voltage. 
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approximation. The C-V dependence was incorporated into the program by defining an 
internal function C(VD) whose value is the capacitance at a given diode voltage Y(l). 

The function routine is handed the value of the diode voltage by subroutine FCT and a 
search is conducted through a series of predefined voltage regions until the one in which 
Y(l) lies, is found. Then, data for the slope and intercept of the C-V straight line ap- 
proximation in that region is used to calculate the value of C(VD). Appendix n contains 

Appendix n contains a list of the modifications which must be made to the program 
in Appendix I to run each of the examples. 

4.2 Computed Results and Discussion 

The simple mixer circuits of Figs. 8 and 9 were analyzed with three different 
diodes: a Ga As Schottky diode (including nonlinear capacitance), a Ga As Sehottky 
diode with constant capacitance, and a practical Ga As Mott diode. In each case the 
diode was forward biased at 0. 4 V and the LO level was adjusted to give a rectified 
dc current of 2.0 mA. The signal frequency (taken as the upper sideband) was 119 GHz 
while the LO and intermediate frequencies were 115 GHz and 4 GHz respectively. The 
conversion loss (upper sideband to IF), equivalent input no se temperature, ard real 
part of the IF output impedance were plotted as functions of the series inductance, 
which was allowed to vary from 0. 01 nH to 0. 25 nH. These results are given in Ap- 
pendix IE, section A III. 1. Note that skin effect was not taken into account in these 
examples. 

The effect of the junction capaci tance on the mixer perfromance was then investi- 
gated for the Schottky diode ( with and without nonlinear capacitance) in the same two 
circuits-Figs. 8 and 9. To isolate the effects of the diode capacitance and series in- 
ductance, the series resistance was assumed constant at 4.4 Ohms. It should be noted 
that this leads to some unrealistically high (and low) diode cutoff frequencies, CO 
was varied from 1 to 20 fF for each of nine values of LS between 0. 04 and 0. 2 nH. 
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Again the noise temperature, conversion loss and output impedance were plotted, this 
time as functions of the junction capacitance for each value of LS (Appendix HI, 

Section A 111-2). 

In Appendix ni, section A HI. 1 , the performance of the Schottky , Mott, and Y = 0 
diode mixers operating with an LO frequency of 115 GHz and an IF of 4 GHz are com- 
pared. It is assumed that each diode has a zero-bias capacitance (CO) of 11- 8 fF. No- 
tice that the minimum noise temperature is achieved with the constant capacitance 
diode; however this should not be assumed to be a general result, as will be demon- 
strated below. Except for the Schottky diode with constant capacitance, the minima 
in the noise temperature and conversion loss for each mixer circuit do not occur at 
the same value of LS. 

A broader view of the performance of these mixers is obtained from the graphs 
in Section A m. 2, in which LS is held constant while CO is varied. In all cases the 
effect of increasing series inductance is to sharpen the noise temperature and conver- 
sion loss minima and shift them towards smaller values of CO. Better performance is 
obtained for larger values of LS, with a corresponding increase in the IF output im- 
pedance. 

Amongst the results there are some points, particularly for low values of series 
inductance, which appear to be randomly scattered. This is due to the fact that for low 
values of the series inductance each increment in LS causes a large change in the reso- 
nant frequency of the diode with the external circuit. If these resonances fall near har- 
monics of the LO the diode waveforms can be strongly affected, and so can the embed- 
ding impedance seen by the small-signal sidebands near these harmonics. 


The LO power required was between 1 and 7 mW in every ease. 
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An interesting result of this analysis is that the parametric effects of the junction 
capacitance do not necessarily degrade the mixer performance. 
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5. Summary and Concluding Remarks 

In this report a user oriented computer program for the analysis of microwave 
and millimeter-wave mixers has been presented. The necessary mixer theory was 
briefly described and used as the basis of the mixer analysis program. The operation 
of the program was discussed as well as some useful modifications. 

As an example of the program's use, a study of the effects of the series induct- 
ance and diode capacitance on the performance of two simple mixer circuits was under- 
taken. Three different diodes were used: a Sehottky diode, a Schottky diode with no 
capacitance variation, and a Mott diode. 

The results show clear minima in the conversion loss and noise temperature as 
the inductance and capacitance are varied. The parametric effects of the voltage depen- 
dent diode capacitance are found to have either beneficial or detrimental effects on the 


mixer performance, depending on the circuit and diode parameters. 
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APPENDIX I 

Listing and Sample Run of the Mixer Analysis Program 

This program calculates the performance of a mixer, given the embedding 
network im p e danc es at the harmonics of the LO and the sideband frequencies. The 
theory is described in section 2 of this report and section 3 outlines the program op- 
eration. The flow chart (Fig. 6) and table of variables (Fig. 7) may be found helpful 
in understanding the workings of the program and are repeated in this appendix for 
convenience. 
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VARIABLE 


k 

L iJ 

q 

n s 

njuci 


n s mujp) 


skin 


V 

Y' 


PROGRAM NAME 
CO 
CtT 

FC(m-n+l) 

GJ 

FG(m-ntl) 

1CJ 

IDCOS{n+l)-JIDSIN(n+I) 
IGJ + ICO 
IGJ 

FGlm-n+U/ALP 

IS 

BOLTZ 

LIJ 

QEL 

US 

Its 

RSLO 

RSSB 

RSKIN 

TEQ 

TM 

VDC 

VLO 

VDCOS (ml) - JVDSIN <n+I> 
Y(l) 

VI DC 

VI (n) « AV (n) - jBV (n) 

VRDC or DVRDC 

VR or DVR 

A 

A 


VARIABLE 
Z*(o/o. IF) 

Z 1 (matched IF) 

Z Q (T. lino imp, ) 

(0) * lydo) 

V n V * V l V 
V"V 


in 


mm 


mm, o° 

a 


5 
e 

IT 

P 

p 0 

6 


*0 

lU 

P 

<&%«£> 

<6i!p 5I’j> 

<|sv N | 2 > 

<! 6V s/> 

tlv d / dt 


PROGRAM NAME 
A 

A 

20 

ZEMBDC 
ZEMB 
ZER + JZEI 
ZEMBSJ3 
ZIFOUT 

A(NHD2Pl-i , NHD2P1-J) 
ZIN 

A(NHD2Pl-m t NHD2Pl-mj 
A(NHD2Pl-m, NIJD2Pl-m) 
ALP 

GAM 

ZQ-X 

ZQACC 

PI 

RHO 

RHODC 

PHI 

VVIF 

WP 

con 

COR 

VSQ 

VSQ 

VSQ 

DERY (1) 


Alph ab etical list of the main variables used in the theory in 
this report, and their counterparts in the mixer analysis program 
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1 . 

2 . 

3. 

MIXER ANALYSIS PROGRAM 4. 

5. 

6 . 

7. 

P.H. SIEGEL AND A. R. KERR 8. 

GODDARD INSTITUTE FOR SPACE STUDIES 9. 

NASA-GODDARD SPACE FLIGHT CENTER 10. 

2880 BROADWAY II. 

NEW YORK, N, Y . 10025 12. 

13. 

14. 

15 . 

GENERAL INFORMATION 16. 

17. 

THIS Pf SRAM ANALYZES MIXERS WITH A SINGLE SCHOTTKY“BARRIER DIODE 18. 
WHOSE I -V AND C-V CHARACTERISTICS ARE KNOWN. ARBITRARY EMBEDDING 19. 

IMPEDANCES AT THE SIDEBAND AND LO HARMONIC FREQUENCIES ARE ALLOWED. 20. 

THE DIODE MOUNT IS ASSUMED LOSSLESS AND RECIPROCAL. 21. 

22 . 

THE PROGRAM IS ORGANIZED INTO TWO MAIN SECTIONS. THE FIRST PER- 23. 

FORMS A NONLINEAR ANALYSIS TO DETERMINE THE DIODE WAVEFORMS PRO- 24. 

DUCED BY THE LOCAL OSCILLATOR. THE SECOND PERFORMS A SMALL-SIGNAL 25. 

AND NOISE ANALYSIS TO COMPUTE THE CONVERSION LOSS, PORT IMPEDANCES, 26. 

AND NOISE TEMPERATURE OF THE MIXER. 27. 

28 • 

THE NONLINEAR ANALYSIS IS BASED ON THE MULTIPLE RELECTION METHOD 29 ! 

(A.R. KERR, IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, 30. 

VOL. MTT-23.N0. 10, PP. 828-831, OCT. 1975) DEVELOPED EARLIER IN OUR 31. 

LABORATORY, MODIFIED TO TAKE INTO ACCOUNT THE NONLINEAR CAPACITANCE 32. 
OF THE DIODE. 33. 

34. 

THE SMALL-SIGNAL ANALYSIS IS BASED ON THAT OF D.N.HELD AND 35. 

A.R. KERR, IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, 36. 

V0L.MTT-26,N0.2,PP .49-61 ,FEB. 1978. 37. 

38. 

39 . 

PROGRAM NOTES 40 ! 

41 . 

TWO MAIN SUBROUTINES CONTROL THE ANALYSIS: LGSIG WHICH PERFORMS 42. 

THE NONLINEAR ANALYSIS AND SMSIG WHICH COMPUTES THE SMALL-SIG- 43. 

NAL AMD NOISE PROPERTIES OF THE MIXER. EACH OF THESE CALLS A NUM- 44. 

BER OF SECONDARY SUBROUTINES WHICH PERFORM SPECIFIC CALCULATIONS 45. 

OR CONTROL THE OUTPUTTING OF RESULTS. 46. 

DATA IS INPUT VIA THE BLOCK DATA SUBPROGRAM. THE FOLLOWING 47. 

INFORMATION MUST BE SUPPLIED BY THE USER: 48. 

49. 

1) THE EMBEDDING IMPEDANCES AT THE LO FREQUENCY AND THE HIGHER 50. 

HARMONICS AS REAL AND IMAGINARY PARTS (ZER.ZEI) IN OHMS. 51. 

2) THE SIDEBAND IMPEDANCES IN COMPLEX FORM ( ZEMBSB ) IN OHMS, WHERE 52. 

SIDEBAND M IS ARRAY ELEMENT (NH/2+1-M) AND THERE ARE NH+1 ARRAY 53. 
ELEMENTS IN ALL. NOTE THAT, BECAUSE ALL LOWER SIDEBANDS ARE TREATED 54. 

AS NEGATIVE FREQUENCIES , VALUES OF ZEMBSB FOR LOWER SIDEBANDS MUST 55. 

BE THE CONJUGATES OF THEIR USUAL POSITIVE FREQUENCY VALUES. 56. 

3) THE LO FREQUENCY ( FP ) AND THE INTERMEDIATE FREQUENCY l IF) IN HZ. 57. 

4} THE DESIRED RECTIFIED CURRENT UDBIAS) IN AMPERES. 58 
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C 5) THE DC BIAS VOLTAGE ACROSS THE DIODE (VDBIAS) IN VOLTS, 59. 
C 6) THE DIODE SERIES RESISTANCE AT DC ( RS > IN OHMS. 60. 
C 7) THE DIODE SKIN RESISTANCE CONSTANT {RSKINJ IN OHMS/SQRT( HZ > . 61. 
C 8) THE DIODE EQUIVALENT TEMPERATURE (TEQ) IN DEGREES K. 62. 
C 9) THE DIODE REVERSE SATURATION CURRENT CIS) IN AMPERES. 63. 
C 101 THE DIODE CAPACITANCE AT ZERO VOLTS (C0) IN. FARADS. 64. 
C 11) THE DIODE CONTACT POTENTIAL (PHI) IN VOLTS. 65. 
C 12) THE DIODE CAPACITANCE LAW EXPONENT (GAM). 66. 
C 13) THE DIODE I-V LAW EXPONENT (ALP=Q/NKT) IN 1/VOLTS. 67. 
C 68. 
C THERE ARE SEVERAL OTHER VARIABLES WHICH MAY BE ADJUSTED TO 69. 
C CONTROL THE OPERATION OF THE PROGRAM. THEIR VALUES HAVE BEEN 70. 
C OPTIMIZED FOR THE LISTING WHICH FOLLOWS AND MAY BE ALTERED 71. 
C WHEN THE PROGRAM IS USED FOR OTHER PROBLEMS. THESE VARAIBLES ARE: 72. 
C 73 . 
C ACCJTHE ACCURACY OF THE RUNGE KUTTA INTEGRATION USED TO SOLVE THE 74! 
C STATE EQUATION OF THE DIODE NETWORK. 75. 
C IDCACCJTHE ACCURACY WITH WHICH THE DC CURRENT MUST APPROACH THE 76. 
C DESIRED VALUE (IDBIAS). 77. 
C NLOsTHE NUMBER OF LO CYCLES NEEDED TO REACH A STEADY STATE IN THE 78. 
C NONLINEAR ANALYSIS ROUTINE. SINCE SETTLING OCCURS IN SUCCESSIVE 79. 
C REFLECTION CYCLES NLO CAN USUALLY BE SET TO ONE. 80. 
C NPRI,NT:THE NUMBER OF CYCLES BETWEEN PRINTOUTS OF THE INTERMEDIATE 81. 
C RESULTS IN THE NONLINEAR ANALYSIS. 82. 
C NPTSs THE NUMBER OF INTERVALS+1 INTO WHICH THE LO CYCLE IS DIVIDED 83. 
C FOR THE INTEGRATION AND STORAGE OF DATA POINTS. TO AVOID 84. 
C ALIASING NPTS SHOULD BE CHOSEN CONSIDERABLY LARGER THAN 85. 
C ( 2*NH+1 ) THE VALUE REQUIRED BY THE SAMPLING THEOREM. 86. 
C VLOJTHE INITIAL VALUE OF THE LOCAL OSCILLATOR VOLTAGE. 87. 
C VLOI NC :THE INITIAL INCREMENT SIZE USED TO ZERO IN ON THE DESIRED 88. 
C DC RECTIFIED CURRENT (IDBIAS). 89. 
C ZQACC : THE DEGREE OF CONVERGENCE OF THE FINAL SOLUTION IN THE NON- 90. 

c linear Analysis. 91 . 

C Z0:THE CHARACTERISTIC IMPEDANCE OF THE HYPOTHETICAL TRANSMISSION 92. 
C LINE INSERTED BETWEEN THE DIODE AND EMBEDDING NLTWORK . Z0 HAS A 93. 
C SIGNIFICANT EFFECT ON THE RATE OF CONVERGENCE OF THE NONLINEAR 94, 
C ANALYSIS. 95. 
C 96. 
C . IF THE DC CURRENT REACHED AFTER TEN REFLECTION CYCLES IN THE 97. 
C NONLINEAR ANALYSIS. IS NOT SUFFICIENCY CLOSE TO IDBIAS, THE CYCLE 98. 
C NUMBER AT WHICH THE LO VOLTAGE IS ADJUSTED MUST BE • INCREASED. THIS 99. 
C IS ACCOMPLISHED BY CHANGING THE STEP AFTER STATEMENT LABEL 10 IN 100. 
C THE LGSIG ROUTINE. 101. 
C 102. 
C THE USER MAY FIND IT NECESSARY TO ALTER OTHER PROGRAM VARIABLES 103. 
C FOR SPECIFIC PROBLEMS. FOR THIS REASON A LIST OF ALL THE VARIABLES 104. 
C (EXCEPT THOSE INTERNAL TO THE IBM'SSP ROUTINES USED IN THE PROGRAM), 105. 
C SUBROUTINES AND COMMON BLOCKS USED IN THE PROGRAM FOLLOWS. 106. 
C 107. 
C 108. 
C LIST OF VARIABLES 109. 
C 110. 
C A; THE SMALL-SIGNAL AUGMENTED ADMITTANCE (Y*) OR IMPEDANCE (Z 1 ) HI. 
C MATRIX OF THE MIXER. 112. 
C ACC: THE INTEGRATION ACCURACY USED IN DRKGS. 113. 
C ALP: THE DIODE I-V LAW EXPONENT ( Q/NKT) . 114. 
C AUX: DRKGS STORAGE ARRAY OF DIMENSION (8,NDIM>. 115. 
C AV: REAL PART OF THE VOLTAGE WAVE INCIDENT ON THE DIODE AT EACH 116. 
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C HARMONIC. 117. 
C BLANK: A NUMERIC USED FOR PLOTTING A BLANK. 118. 
C BOLTZ: BOLTZMANN'S CONSTANT (1.38 E-23). 119. 
C BV: IMAGINARY PART OF THE VOLTAGE WAVE INCIDENT ON THE DIODE AT 120. 
C EACH LO HARMONIC. 121. 
C CJ: FREQUENCY SCALED DIODE JUNCTION CAPACITANCE USED IN LGSIG. 122. 
C COCOS: FOURIER COSINE COEFFICIENTS OF THE DIODE CAPACITANCE. 123. 
C CODATA: STORAGE ARRAY CONTAINING THE DIODE CAPACITANCE FOR EACH 124. 
C OF THE NPTS POINTS IN THE LOCAL OSCILLATOR CYCLE. 125. 
C COMAG: THE MAGNITUDES OF THE FOURIER CAPACITANCE COEFFICIENTS. 126. 
C COPHA: PHASES OF THE FOURIER CAPACITANCE COEFFICIENTS (IN DEGREES). 127. 
C COPOS: POSITION OF CJ IN THE PLOT OF DIODE CAPACITANCE. 128. 
C COSIN: FOURIER SINE COEFFICIENTS OF THE DIODE CAPACITANCE. 129. 
C COR*. THE NOISE CURRENT CORRELATION MATRIX USED IN THE NOISE ANALYSIS. 130. 
C C0: THE DIODE CAPACITANCE AT ZERO VOLTS (IN FARADS). 131. 
C C0POS: POSITION OF C0 IN THE GRAPH OF THE DIODE CAPACITANCE. 132. 
C DERY: INITIALLY THE RKGS ERROR PARAMETER AND LATER THE DERIVATIVE 133. 
C IN THE NETWORK STATE EQUATION (DY(1)/DT). 134. 
C DET: DETERMINANT OF THE A MATRIX AS RETURNED BY THE CMINV ROUTINE. 135. 
C DOT: A NUMERIC USED FOR PLOTTING A DOT. 136. 
C DVR: DIFFERENCE IN THE REFLECTED VOLTAGE BETWEEN TWO CYCLES. 137. 
C DVRDC : DIFFERENCE IN’ THE DC REFLECTED VOLTAGE BETWEEN TWO CYCLES. 138. 
C FC: COMPLEX FOURIER COEFFICIENTS OF THE DIODE CAPACITANCE (HALF 139. 
C ’AMPLITUDE ) . 140. 
C FG: COMPLEX FOURIER COEFFICIENTS OF THE DIODE CONDUCTANCE (HALF 141. 
C AMPLITUDE). 142. 
C FP: THE LOCAL OSCILLATOR OR PUMP FREQUENCY IN HERTZ. 143. 
C GAM: THE DIODE CAPACITANCE EXPONENT. 144. 
C GJ: THE DIODE CONDUCTANCE. 145. 
C GJCOS: FOURIER COSINE COEFFICIENTS OF THE DIODE CONDUCTANCE. 146. 
C GO DATA : STORAGE ARRAY CONTAINING THE VALUES OF THE DIODE 147. 
C CONDUCTANCE FOR EACH OF THE NPTS POINTS IN THE LO CYCLE. 148. 
C GOMAG: MAGNITUDES OF THE FOURIER CONDUCTANCE COEFFICIENTS. 149. 
C GJPHA: PHASES OF THE FOURIER CONDUCTANCE COEFFICIENTS (IN DEGREES). 150. 
C GJSIN : FOURIER SINE COEFFICIENTS OF THE DIODE CONDUCTANCE. 151. 
C I: ARRAY ELEMENT OR COUNTING VARIABLE. 152. 
C ICJ i THE CURRENT THROUGH THE DIODE CAPACITANCE. 153. 
C ICUDAT: STORAGE AFRAV FOR ICJ AT EACH POINT IN THE LO CYCLE. 154. 
C IDBIAS: DESIRED RtCTIF IED CURRENT AT WHICH THE MIXER IS TO BE 155. 
C OPERATED ( IN AMPS). 156. 
C IDCACC : DESIRED ACCURACY OF THE CALCULATED DC CURRENT, MEASURED AS 157. 
C THE MAXIMUM TOLERABLE DEVIATION FROM THE DESIRED DC CURRENT 158. 
C , IDBIAS. 159. 
C IDCOS : FOURIER COSINE COEFFICIENT OF THE TOTAL DIODE CURRENT. 160. 
C IDDATA : STORAGE ARRAY CONTAINING THE VALUES OF THE TOTAL DIODE 161. 
C CURRENT FOR EACH OF THE NPTS POINTS IN THE LO CYCLE, 162. 
C IDPOS : POSITION OF ID ON THE GRAPH OF TOTAL CURRENT IN THE DIODE. 163. 
C IDSIN: FOURIER SINE COEFFICIENT OF THE TOTAL DIODE CURRENT. 164. 
C IER: THE ERROR MESSAGE CODE OF SUBROUTINE DFORIT . 165. 
C IF: THE INTERMEDIATE FREQUENCY IN HERTZ. 166. 
C IGJ : THE CURRENT THROUGH THE DIODE CONDUCTANCE. 167. 
C IGODAT: STORAGE ARRAY CONTAINING VALUES OF THE DIODE CONDUCTANCE 168. 
C CURRENT FOR EACH OF THE NPTS POINTS IN THE LO CYCLE. 169. 
C IGJPOS : POSITION OF IGJ ON THE GRAPH OF DIODE CONDUCTANCE CURRENT. 170. 
C IHLF : DRKGS PARAMETER GIVING THE NUMBER OF TIMES THE INTEGRATION 171. 
C INTERVAL HAS BEEN HALVED OR DOUBLED. 172. 
C IMZIN: IMAGINARY PART OF THE INPUT IMPEDANCE AT EACH SIDEBAND. 173. 
C IS: DIODE SATURATION CURRENT (IN AMPS). 174. 
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C ITER: LOOP REPRESENTING THE REFLECTION CYCLES IN THE NONLINEAR 17S. 
C ANALYSIS 176. 
C IVLO: COUNTING VARIABLE USED TO SIGNAL THE END OF THE JLO LOOP. 177. 
C J: ARRAY ELEMENT OR COUNTING VARIABLE. 178. 
C OH: LOOP OVER THE LO HARMONICS. 179, 
C JLO: LOOP OVER THE LO CYCLES OF INTEGRATION TO OBTAIN STEADY STATE. 180. 
C OPRINT: KEEPS TRACK OF THE CYCLES WHOSE RESULTS ARE TO BE PRINTED. 181. 
C OPT: LOOP OVER THE INTERVALS FOR WHICH DRKGS WILL INTEGRATE. 182. 
C OVLO: LOOP FOR ZEROING IN ON THE DESIRED RECTIFIED CURRENT {IDBIAS). 183. 
C K: ARRAY ELEMENT OR COUNTING VARIABLE. 184. 
C LIO: CONVERSION LOSS FROM SIDEBAND 0 TO SIDEBAND I. 185. 
C LOFLAG: A FLAG TO KEEP COUNT OF THE NUMBER OF TIMES VLO HAS BEEN 186. 
C LOWERED TO TRY AND OBTAIN THE DESIRED RECTIFIED CURRENT. 187. 
C LOPWR: THE REQUIRED LO POWER FOR A RECTIFIED CURRENT OF IDBIAS. 188, 
C LOVLO: THE! LOWER BOUND OF THE VLO ADJUSTMENT DONE IN SUBROUTINE 189. 
C ADJLO. 190. 
C MAXID : MAXIMUM VALUE OF ID IN AN LO CYCLE {USED FOR PLOTTING). 191. 
C MAXIGJ : MAXIMUM VALUE OF IGJ IN AN LO CYCLE. (USED FOR PLOTTING). 192. 


C MAXCJ : MAXIMUM VALUE OF CJ IN AN LO CYCLE (USED FOR PLOTTING). 193. 
C MAXVD: MAXIMUM VALUE OF VD IN AN LO CYCLE {USED FOR PLOTTING). 194. 
C MINCO: MINIMUM VALUE OF CO IN AN LO CYCLE (USED FOR PLOTTING). 195. 
C NDIM: THE NUMBER OF DIFFERENTIAL EQUATIONS TO BE SOLVED BY DRKGS. 196. 
C NH: TOTAL NUMBER OF LO HARMONICS CONSIDERED . IN THE ANALYSIS. 197. 
C NHARM: USED IN PRINT1 AS THE TOTAL NUMBER OF LO HARMONICS. 198. 


C NHD2: THE TOTAL NUMBER OF LO HARMONICS DIVIDED BY TWO. 199. 
C NHD2P1: THE TOTAL NUMBER OF LO HARMONICS DIVIDED BY TWO, PLUS ONE. 200. 
C NHP 1 : THE TOTAL NUMBER OF LO HARMONICS PLUS ONE {# OF SIDEBANDS). 201. 
C NHP2 : THE NUMBER OF LO HARMONICS PLUS 2. ' 202. 
C NITER: TOTAL NUMBER OF ITERATIONS ALLOWED FOR ACHIEVING FULL CONVER- 203. 
C GENCE AFTER THE DESIRED RECTIFIED CURRENT HAS BEEN OBTAINED. 204. 
C NLO: THE TOTAL NUMBER OF LOCAL OSCILLATOR CYCLES TO BE INTEGRATED 205. 
C THROUGH IN ORDER TO REACH A STEADY STATE. 206. 
C NPRINT : CONTROLS THE REFLECTION CYCLE NUMBER AT WHICH THE PRINTING 207. 
C OF THE INTERMEDIATE RESULTS OCCURS. 208. 
C NPTS: THE NUMBER OF INTERVALS + 1 OVER WHICH DRKGS INTEGRATES. 209. 
C NSB : SIDEBAND TO WHICH THE EQUIVALENT NOISE TEMPERATURE OF THE 210. 
C MIXER IS REFERRED TO IN SUBROUTINE TMIX. 211. 
C NVLQ: THE NUMBER OF TIMES THE LO VOLTAGE IS ALLOWED TO BE ADJUSTED. 212. 
C PHI: DIODE CONTACT POTENTIAL {IN VOLTS). 213. 
C Pi: THE CONSTANT PI. 214. 
C PRMT: ARRAY USED BY RKGS CONTAINING THE FOLLOWING VARIABLES: 215. 
C PRMT(l): THE LOWER BOUND o’F THE INTEGRATION INTERVAL. 216. 
C PRMT(2): THE UPPER BOUND OF THE INTEGRATION INTERVAL. 217. 
C PRMT( 3 ) : THE INITIAL INTEGRATION INTERVAL STEP SIZE. 218. 
C PRMTC4): THE INTEGRATION ACCURACY. 219. 
C PRMTC5): ALLOWS HALTING OF THE RKGS ROUTINE IF SET TO 0. 220. 
C QEL : THE CHARGE OF THE ELECTRON (1.6 E-19 COULOMBS). 221. 
C REZIN: REAL PART OF THE. INPUT IMPEDANCE AT EACH SIDEBAND. 222. 
C RHO: THE COMPLEX VOLTAGE REFLECTION COEFFICIENT OF THE EMBEDDING 223. 
C NETWORK (INCLUDING RS) AT EACH LO HARMONIC. 224. 
C RHODC : THE DC REFLECTION COEFFICIENT OF THE EMBEDDING NETWORK {IN- 225. 
C CLUDING RS). 226. 
C RS: LOW FREQUENCY DIODE SERIES RESISTANCE MEASURED AT A HIGH ENOUGH 227. 
C FREQUENCY THAT THERMAL EFFECTS ARE NEGLIGIBLE BUT SKIN EFFECTS 228. 
C HAVE NOT YET BECOME’ SIGNIFIGANT. 229. 
C RSKIN : CONSTANT WHICH DETERMINES THE SKIN RESISTANCE FOR A 230, 
C PARTICULAR DIODE, WHERE RS=RS+RSKI N*SQRT( F ) AT FREQUENCY F. 231. 
C RSLO: VALUES OF RS AT THE LO HARMONICS. 232. 


c RssB: values of rs at the sideband frequencies.* 233 . 

C SHOTi THE SHOT NOISE CONTRIBUTION TO THE TOTAL MIXER TEMPER.' TURE . 234. 
C STAR: A NUMERIC USED FOR PLOTTING AN ASTERISK. 235. 
C T: INTERMEDIATE VARIABLE IN THE MIXER TEMPERATURE CALCULATION. 236. 
C TEQt THE EQUIVALENT TEMPERATURE OF THE DIODE SERIES RESISTANCE 237. 
C INCLUDING THE EFFECTS OF LATTICE SCATTERING AND LO HEATING. 238. 
C THERM: THE THERMAL NOISE CONTRIBUTION TO THE TOTAL MIXER TEMPERATURE. 239. 
C TM: THE TOTAL EQUIVALENT INPUT NOISE TEMPERATURE OF THE MIXER. 240. 
C UPFLAG: A FLAG TO KEEP COUNT OF THE NUMBER OF TIMES VLO HAS BEEN 241. 
C RAISED TO TRY AND GET THE. DESIRED RECTIFIED CURRENT IDBIAS. 242. 
C UPVLO: THE UPPER BOUND OF VLO DURING ADJUSTMENTS DONE IN ADJIO. 243. 
C V: THE SUM OF THE INSTANTANEOUS SOURCE VOLTAGES AT ALL LO HARMONICS. 244. 
C VDBIAS: THE DC BIAS VOLTAGE ACROSS THE DIODE AND SERIES RESISTANCE. 245. 
C VDC : THE DC VOLTAGE SEEN BY THE DIODE. 246. 
C VDCOS: FOURIER COSINE COEFFICIENTS OF THE VOLTAGE ACROSS THE DIODE. 247. 
C VODATA: STORAGE ARRAY CONTAINING YU), THE VOLTAGE ACROSS THE DIODE 248. 
C AT EACH OF THE NPTS POINTS IN THE LO CYCLE. 249. 
C VDINIT: INITIAL GUESS FOR THE VOLTAGE ACROSS THE DIODE. 250. 
C VDPOS : POSITION OF VD ON THE GRAPH OF THE VOLTAGE ACROSS THE DIODE. 251. 
C VDSIN: FOURIER SINE COEFFICIENTS OF THE VOLTAGE ACROSS THE DIODE. 252. 
C Vis COMPLEX VOLTAGE WAVE INCIDENT ON THE DIODE AT EACH HARMONIC. 253. 
C VIDC : DC VOLTAGE INCIDENT ON THE DIODE. 254. 
C VLO: AMPLITUDE OF THE LOCAL OSCILLATOR VOLTAGE. 255. 
C VLOINC: AMOUNT BY WHICH THE LO VOLTAGE IS INCREMENTED WHEN TRYING 256. 
C TO OBTAIN THE DES'RED DC CURRENT IDBIAS. 257. 
C VR: COMPLEX VOLTAGE WAVE REFLECTED PROM THE DIODE AT EACH HARMONIC. 258. 
C VRDC : DC VOLTAGE WAVE REFLECTED FROM THE DIODE AT EACH HARMONIC. 259. 
C VSQ: MEAN SQUARE OUTPUT NOISE VOLTAGE. 260. 
C VX: INTERMEDIATE VARIABLE FOR THE COMPLEX VOLTAGE WAVEFORM. 261. 
C VXDC : INTERMEDIATE VARIABLE FOR THE DC VOLTAGE WAVEFORM. 262. 
C WIF: 2*P I*INTERMEDIATE FREQUENCY. 263. 
C WKl: WORK SPACE USED IN THE MATRIX INVERSION ROUTINE CMINV. 264. 
C WK2 : WORK SPACE USED IN THE MATRIX INVERSION ROUTINE CMINV. 265. 
C WP: 2*P I *PUMP FREQUENCY. 266. 
C X: THE DEPENDENT VARIABLE IN DRKGS ( X=2*P I*FP*TIME ) . 267. 
C XL MAT : THE CONVERSION LOSS MATRIX WHICH GIVES THE CONVERSION 268. 
C LOSSES BETWEEN PAIRS OF SIDEBANDS. 253. 
C Y: DRKGS VARIABLE TO BE FOUND ( Y=VOLTAGE ACROSS THE DIODE WITHOUT 270. 
C THE SERIES RESISTANCE). 271. 
C YCPOS : USED FOR PLOTTING THE DIODE CAPACITANCE. 272. 
C YGPOS: USED FOR PLOTTING THE CURRENT THROUGH THE DIODE CONDUCTANCE. 273. 
C YIDPOS: USED FOR PLOTTING THE TOTAL DIODE CURRENT. 274. 
C YPT: A DO LOOP VARIABLE USED FOR PLOTTING POINTS ACROSS A PAGE. 275. 
C YVDPOS : USED FOR PLOTTING THE VOLTAGE ACROSS THE DIODE. 276. 
C ZD: THE EFFECTIVE DIODE IMPEDANCE AT EACH HARMONIC (VD/ID). 277. 
C ZEI: IMAGINARY PART OF THE EMBEDDING IMPEDANCE AT EACH LO HARMONIC. 278. 
C ZEMB: COMPLEX EMBEDDING IMPEDANCE AT EACH LO HARMONIC INCLUDING 279. 
C THE DIODE SERIES RESISTANCE. 280. 
C ZEMBSB : THE EMBEDDING IMPEDANCES AT THE SIDEBAND FREQUENCIES. 281. 
C ZER: REAL PART OF THE EMBEDDING IMPEDANCE AT EACH LO HARMONIC. 282. 
C ZERO: A NUMERIC USED FOR PLOTTING AN 0. 283, 
C ZIFOUT: IF PORT IMPEDANCE INCLUDING THE DIODE SERIES RESISTANCE. 284. 
C ZIN: INPUT IMPEDANCES AT THE SIDEBANDS INCLUDING THE DIODE 285. 
C SERIES RESISTANCE.. 286. 
C ZQ: IMPEDANCE QUOTIENT (ZD/ZE) AT EACH HARMONIC (USED AS THE LARGE 287. 
C SIGNAL CONVERGENCE PARAMETER). 288. 
C ZQACC : DESIRED DEGREE OF CONVERGENCE MEASURED AS THE DEVIATION FROM 289. 
C UNTIY OF (THE IMPEDANCE AT THE DIODE/EMBEDDING IMPEDANCE). 290. 
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ZQFLAG: THE NUMBER OF HARMONICS WHICH HAVE NOT VET CONVERGED IN A 291, 
PARTICULAR CYCLE OF THE NONLINEAR ANALYSIS. 292. 

ZQMAG ! MAGNITUDE OF THE IMPEDANCE QUOTIENT <ZQ). 293, 

ZQPHA: PHASE OF THE IMPEDANCE QUOTIENT (ZQ) IN DEGREES. 294, 

Z0: CHARACTERISTIC IMPEDANCE OF THE HYPOTHETICAL TRANSMISSION LINE 295. 
INSERTED BETWEEN THE DIODE AND THE EMBEDDING NETWORK. 296. 

297. 

298 

LIST OF SUBROUTINES 299! 

300 , 

ADOLOs SUBROUTINE FOR ADJUSTING THE LOCAL OSCILLATOR VOLTAGE UNTIL 301. 

THE DC CURRENT IS CLOSE ENOUGH TO THE DESIRED VALUE , IDBIAS . 302. 

CMINV : IBM SSP PROGRAM SLIGHTLY MODIFIED TO INVERT A COMPLEX MATRIX. 303. 

CORRELs SUBROUTINE FOR FORMING THE NOISE CORRELATION MATRIX. 304. 

DFORITs DOUBLE PRECISION VERSION OF AN IBM SSP PROGRAM WHICH 305. 

PERFORMS A FOURIER ANALYSIS ON A PERIODIC WAVEFORM. 306. 

DRKGSJ AN IBM SSP PROGRAM WHICH SOLVES A SYSTEM OF DIFFERENTIAL 307. 

EQUATIONS USING A RUNGE-KUTTA ALGORITHM. 30B . 

FCTs SUBROUTINE, FOR USE WITH DRKGS .CONTAINING THE NETWORK STATE 309. 

EQUATION. 310, 

LGSIG: CONTROL PROGRAM FOR THE LARGE-SIGNAL , NONL I NEAR ANALYSIS. 311. 

OUTP: OUTPUT ROUTINE REQUIRED BY DRKGS BUT NOT USED IN THIS PROGRAM. 312. 

PLOTS SUBROUTINE FOR PLOTTING THE LARGE-SIGNAL DIODE WAVEFORMS. 313. 

POWER: CALCULATES THE LO POWER REQUIRED FOR THE ACTUAL MIXER. 314. 

PRINT1 : CONTROLS THE PRINTING OF INITIAL VALUES. 315. 

PRINT2 : CONTROLS THE PRINTING OF VARIABLES FOR A PARTICULAR CYCLE 316. 

OF THE NONLINEAR ANALYSIS. 317. 

PRINT3 : WRITES THE FOURIER COEFFICIENTS OF THE DIODE CONDUCTANCE 318. 

AND CAPACITANCE, 319. 

PRINT4: PRINTS THE RESUi THE SMALL-SIGNAL AND NOISE ANALYSIS. 320. 

SMSIG: CONTROL PROGRAM F " HE SMALL-SIGNAL AND NOISE ANALYSIS. 321. 

TMIX: SUBROUTINE FOR CALL . ..TING THE EQUIVALENT NOISE : TEMPERATURE 322. 

OF THE MIXER. 323. 

TNOISE: ADDS THE THERMAL NOISE COMPONENT TO THE CORRELATION MATRIX. 324. 

YPRIME : SUBROUTINE FOR FORMING THE AUGMENTED MATRIX Y'. 325. 

ZEMBED: FORMS THE EMBEDDING IMPEDANCES AT THE HARMONICS OF THE LO 326. 

AND AT THE HARMONIC SIDEBANDS AND THE SERIES RESISTANCE 327. 

AT THESE FREQUENCIES. 328. 

329. 

330. 

LIST OF COMMON BLOCKS 331. 

332 

COMMON /AD JVLO/ : CONTAINS VARIABLES USED IN THE ADJUSTMENT OF THE LO 333*. 

VOLTAGE TO GIVE A DC CURRENT OF IDBIAS. 334. 

COMMON /CON ST/ : CONTAINS CONSTANTS USED IN THE PROGRAM. 335. 

COMMON/D I OD.E/S CONTAINS DIODE PARAMETERS. 236. 

COMMON /FOR ITS/ : CONTAINS VARIABLES RETURNED BY DFORIT. 337. 

COMMON/ IMP ED/ : CONTAINS INPUT EMBEDDING IMPEDANCES AND THE DIODE 338. 

SERIES RESISTANCE AT THE VARIOUS FREQUENCIES. 339. 

COMMON / I VMAG/: CONTAINS THE FOURIER COEFFICIENTS OF THE DIODE CON- 340. 

DUCTANCE AND CAPACITANCE IN MAGNITUDE AND PHASE.. 341. 

COMMON/LOOPS/ : CONTAINS THE LIMITS OF THE VARIOUS PROGRAM LOOPS. 342. 

COMMON/RKG/ CONTAINS THE INITIAL VALUES FOR THE DRKGS INTEGRATION. 343. 

COMMON/TLINE/ :CONTAINS PARAMETERS EFFECTING THE CONVERGENCE OF THE 344. 

NONLINER ANALYSIS. 345. 

COMMON /VOLTS/ : CONTAI NS VALUES OF THE CIRCUIT VOLTAGES AND CURRENTS. 346. 

347. 

348. 
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BEGIN THE MIXER ANALYSIS, 


— -'CALL LGSIG TO DO THE LARGE SIGNAL ANALYSIS 
CALL LGSIG 

C-— CALL SMSIG TO DO THE SMALL-SIGNAL AND NOISE ANALYSIS 
CALL SMSIG 
STOP 
END 


SUBROUTINE LGSIG 

LGSIG PERFORMS THE NONLINEAR ANALYSIS OF THE MIXER TO DETERMINE 
THE DIODE CONDUCTANCE AND CAPACITANCE WAVEFORMS PRODUCED BY THE 
LOCAL OSCILLATOR* THE MULTIPLE REFLECTION METHOD IS USED , ( A * R . KERR , 
IEEE TRANS .ON MICROWAVE THEORY AND TECH * , VOL .MT T -23 , NO * 10 , PP . 828-83 1 , 
OCT. 1975) MODIFIED TO INCLUDE THE NONLINEAR CAPACITANCE OF THE DIODE. 
ALL INITIALIZED VARIABLES ARE INPUT THROUGH THE BLOCK DATA PROGRAM 
AND ARE TRANSFERRED THROUGH COMMON STATEMENTS. 

THE OUTPUT INCLUDES: 

1) VALUES OF THE INITIALIZED VARIABLES AND INPUT DATA. 

2) THE EMBEDDING IMPEDANCES AND DIODE SERIES RESISTANCE AT THE 
LO HARMONICS AND SIDEBAND FREQUENCIES. 

3) THE RESULTS OF SOME OR ALL OF THE REFLECTION CYCLES. 

4) THE RESULTS OF THE FINAL REFLECTION CYCLE. 

5 > THE REQUIRED LO POWER OF THE MIXER. 

6) THE FOURIER COEFFICIENTS OF THE DIODE CONDUCTANCE AND CAPACITANCE. 

7) PLOTS OF THE DIODE VOLTAGE, TOTAL CURRENT, CONDUCTANCE CURRENT, 

AND CAPACITANCE OVER A SINGLE LO CYCLE. 

TO CHANGE THE DIMENSIONS OF ANY ARRAY SIMPLY CHANGE THE APPROPRIATE 
VARIABLE IN THE TYPE STATEMENT IN THIS ROUTINE , SUBROUTINE SMSIG, AND 
SUBPROGRAM BLOCK DATA. (AV AND BV MUST ALSO HAVE THEIR DIMENSIONS 
ADJUSTED IN SUBROUTINE FCT ) . 

THE SUBROUTINES DRKGS t DFQRIT AND CMINV ARE IBM SSP PROGRAMS AND 
ALL EXCEPT DRKGS HAVE BEEN SLIGHTLY ALTERED FOR USE IN THIS ANALYSIS. 


-THE VARIABLES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 

-FOR COMMON /A DO VLO/ : 

REALMS LOVLO,UPVLO t VLOINC, IDCACC 
INTEGER UPFLAG , LQFLAG 

C FOR COMMON/CONST/: 

REAL*8 QEL , BOLTZ, P I ,TEQ 
C---FOR COMMON/DIODE / % 

REAL *8 ALP , PH I , GAM , C0, I S , RS , RSKIN , FP , IF t IGO , ICJ , GO ,CO 
C — FOR COMMON/FORITS/: 

REAL *8 GOCOSC 9 ) ,G0 S 7, N < 9 ) , COCOS ( 9 ) ,CUSIN( 9 ) , VDCOS ( 9 ) , VDS I N{ 9 ) 
REALMS I DCOS{ 3 ) , I OS I N ( 9 ) 

INTEGER IER 

C— -FOR COMMON/IMPED/: 

COMPLEX.* 16 ZEMBSB ( 9 ) 

REAL *8 ZER( 0 > , ZE I < 8 ) , ZEMBDC , RSLO{ 8 ) , RSSB( 9 ) 

C — FOR COMMON/LOOPS/: 

INTEGER NH, NLO, NVLO , NPTS , NP R1 NT t NITER 

C FOR COMMON/ R KG/ : 

REALMS ACC , VDINIT 
INTEGER NDIM 
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C FOR COMMON/TUNE/: 407. 

REAL*8 Z0,ZQACC 408. 

INTEGER ZQFLAG 409. 

C FOR COMMON/VOLTS/: 410. 

REAL*8 AV( 3 ) ,BV( 8 ), VIDC , VLO, VDBIAS, I OB IAS 411. 

C— FOR VARIABLES NOT IN ANY COMMON BLOCKS: 412. 

COMPLEX* 1 6 RHO( Q ) ,ZEM&( 8 ) , VI ( 8 ) , VR( 8 > , DVR( 8 ) , VX ,ZD ,ZQ 413. 

REAL*8 Y<1),DERY(1> , PRMT< 5 ) , AUX ( 8 , 1 ) 414. 

REAL*8 VRDC , OVRDC , VXDC , VDC , RHODC , LOPWR 415. 

REAL*8 ICJDATt 51 > , IGJDATt 51 ) , CO DATA! 51) , GJDATAI 51 ) 416. 

REAL*8 VDDATAt 51 ) , IDDATA( 51 ) ,ZQMAG( 8 ) ,ZQPHA( 8 ) 417. 

INTEGER IHLF , ITER, IVLO, JVLO, JLO.JPT, JH , NHP 1 , NHD2, NHD2P1 418. 

C— THE COMMON BLOCKS USED ARE: 419. 

COMMON/ADJVLO/LOVLO.UPVLO, LOFLAG , UP FLAG , VLO INC , IDCACC 420. 

COMMON/CONST/QEL , BOLTZ, P I , TEQ 421 . 

COMMON/ DIODE /ALP , PHI ,GAM , C0, IS , RS ,RSKIN , FP , IF , IGJ , ICJ ,GJ ,CJ 422. 

COMMGM/FORITS/GJCOS,GJSIN,CJCOS,CJSIN,VDCOS,VDSIN, IDCOS,IDSIN,IER 423. 
COMMOjv /IMPED/ZER.ZE I ,ZEMBDC,ZEMBSB, RSLO, RSSB 424. 

COMMON /LOOP S/NH,NLO,NVLO,NPTS,NPR I NT, NITER 425. 

COMMON/RKG/ACC,VDINIT,NDIM 426. 

COMMON/TL I NE /Z0, ZQACC , ZQFLAG 427 . 

COMMON/VOLTS/AV.BV, VIDC, VLO, VDBIAS, I DB IAS 428. 

C— -SINCE THE FCT AND OUTP SUBPROGRAMS ARE CALLED BY DRKGS THEY MUST BE 429. 
C— - DEFINED EXTERNALLY 430. 

EXTERNAL FCT, OUTP 431. 

C— DEFINE SOME USEFUL CONSTANTS 432. 

NHP1=NH+1 433. 

NHD2=NH/2 434. 

NHD2P l=NH/2+l 435. 

C— CALL ZEMBED TO FORM THE EMBEDDING IMPEDANCES 436. 

CALL ZEMBEDIZER.ZEI , ZEMBDC ,ZEMBSB , RS, RSLO. RSSB, RSKIN.FP, IF 437. 

1 , P I .NH.NHPl .NHD2P1 ) 438. 

C— -SET THE IMPEDANCE AT DC AND THE FIRST HARMONIC TO Z0 TO SPEED THE 439. 
C-— ANALYSIS. THIS DOES NOT EFFECT THE DIODE WAVEFORMS PROVIDED THE 440. 

C— BIAS VOLTAGE IS ADJUSTED TO GIVE THE DESIRED DC VOLTAGE AT THE 441. 

C— DIODE TERMINALS. THIS IS DONE 5 LINES BEFORE STATEMENT LABEL 2. 442. 

ZEMBt 1 )=DCMPLX( Z0.0.0D0 ) 443. 

ZEMBDC=Z0 444. 

C— FORM THE SET OF COMPLEX IMPEDANCES WITH THE SERIES RESISTANCE ADDED 445. 
DO 1 JH=2 , NH 446. 

1 ZEMB(JH>=DCMPLX(ZER<JH>+RSLOCJH>,ZEI(JH>> 447. 

C— CALCULATE THE REFLECTION COEFFICIENT OF THE EMBEDDING NETWORK AT 448. 

C EACH LO HARMONIC 449. 

RHODC=< ZEMBDC-Z0)/( ZEMBDC +Z0 ) 450. 

DO 13 JH=1 ,NH 451. 

13 RHO< JH)=xZEMB( JH>-Z0)/{ZEMB< JH)+Z0> 452. 

C— INITIALIZE THE VARIABLES FOR THE VLO ADJUSTMENT LOOP 453. 

OVLO= 1 454. 

IVLO=NVLO 455. 

UPFLAG=0 456. 

LOFLAG=0 457. 

C-— INITIALIZE VARIABLES FOR THE INTEGRATION BY DRKGS 458. 

PRMT( 4 )=ACC 459. 

PRMTI 3 )=2 .0D0*PI /DFLOATt NPTS ) 460. 

PRMT{ 2,)fPRMT( 3 ) 461. 

PRMTI 1 1-0.0D0 462. 

YU)=VDINIT 463. 

DERYt 1 1 = 1 .0D0 464. 
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C— CALCULATE THE DC SOURCE VOLTAGE FROM THE GIVEN BIAS VOLTAGE , VDB IAS , 465, 

C— -ACROSS THE DIODE PLUS SERIES RESISTANCE 466. 

VDC=VDBIAS+IDBIAS*(ZEMBDC-RS ) 467. 

C— SET THE INITIAL INCIDENT AND REFLECTED VOLTAGES 468. 

DO 2 OH»l , NH 469. 

VI (OH >=DCMPLXC0.0D0,0.0D0> 470. 

VR( OH >=DCMPLX<0. 008 , 3 . 800 ) 471. 

AV(OH )=0.0D0 472. 

2 BV(JH)=0.0D0 473. 

C— THE DC TERMS 474. 

VRDC=0.0l>0 475. 

VIDC=VDC*Z0/<Z0+ZEMBDC > 476. 

C— RETURN HERE IF THE LO VOLTAGE HAS BEEN ADOUSTED 477. 

15 ITER=0 478. 

C— SET THE INCIDENT VOLTAGE AT THE LO FREQUENCY AND STORE THE RESULTS 479. 
VX=VL 0*20/ ( ZEMB ( 1 ) +20 ) 4 B0 . 

AV( 1 )=DREAL(VX> 481. 

BV( 1 )=-DIMAG< VX ) 482. 

VI(1)=VX 483. 

C CALL PRINT 1 TO WRITE THE INITIAL CONDITIONS 484. 

IF(OVLO.NE.l) GOTO 3 485. 

CALL PRINT1 ( ZEMB .ZEMBDC ,ZER ,ZE I , ZEMB S'B , PRMT , Y ,DERY , VLO , VDB IAS 485 . 

I,IDBIAS,RSSB,RSL0,NH,NHP1,NHD2) 487. 

C— START THE REFLECTION CYCLE 488. 

3 ITER=ITER+1 489. 

C— PRINT ONLY AFTER MULTIPLES OF NPRINT CYCLES HAVE BEEN COMPLETED 490. 

OPRINT=MOD( ITER.NPRINT) 491. 

C BEGIN THE LO CYCLE LOOP FOR REACHING STEADY STATE BETWEEN THE 492. 

C DIODE AND TRANSMISSION LINE AND SET THE INTEGRATION INTERVAL 493. 

DO E OLO=1 ,NLO 494. 

PRMT( 1 )=0.0D0 495. 

PRMT( 2 J=PRMT( 3 ) 496. 

DO 5 0PT=1,NPTS 497. 

CALL DRKGS(PRMT,Y,DERY,NDIM, IHLF , FCT, OUTP ,AUX ) 498. 

C— RESET THE DRKGS ERROR WEIGHT FOR THE NEXT INTERVAL 499. 

DERYt 11=1 .0D0 500. 

I F( OLO. NE .NLO } GOTO 14 501. 

IF ( JPT.EQ.NPTS) GOTO 4 502. 

C— STORE THE RESULTS FOR EACH POINT STARTING AT TIME T=0 503. 

VDDATA( JPT+1 )=Y( 1 ) 504. 

ID DATA ( OPT+1 )=IGO+ICO 505. 

IGODATlOPT+l ) = IGO 506. 

IC0DAT(0PT+1 ) = ICO 507. 

GJDATA( OPT+1 )=G0 508 . 

CODATA(OPT+1 }=C0 509. 

GOTO 14 510. 

4 VDDATA( 1 >=Y< 1 ) 511 . 

IDDATA< 1 ) = IGO + ICJ 512. 

GODATA( 1 )-G0 513. 

CODATAt I >=C0 514. 

IGODAT( 1 ) = IGJ 515. 

ICODATC 1 ) = ICJ 516. 

C GO ON TO THE NEXT TIME INTERVAL 517. 

14 PRMTl i )=PRMT( 2 ) 518. 

PRMT( 2 )=PRMT( 2 )+PRMT< 3 ) 519. 

5 CONTINUE 520. 

6 CONTINUE 521 . 

C— CALL DFORIT TO FORM THE FOURIER COEFFICIENTS OF THE DIODE CURRENT 522. 
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C AND VOLTAGE 523, 

CALL DFORITI VDDATA, NPTS/2 ,NH,VDCOS,VDSIN,IER) 524. 

CALL DFORITI I DDATA , NPTS/2 , NH , I DCOS , IDSI N , IER) 525. 

C — SET THE FLAG FOR THE CONVERGENCE TESTS 526. 

2QFLAG°0 527. 

C---FIND THE DIODE IMPEDANCES AND CALCULATE THE REFLECTED VOLTAGE 528, 

C-“ -WAVE FORMS 529. 

DO 7 OHM ,NH 530. 

ZD d DCMPLX( VDCOSI JH+ 1 ) , -VDSINI OH+1 >)/DCMPLX( IDCOSIOH+l ), 531. 

1-lDSINtJH+U) 532. 

VX=DCMPLX( VDCOSI OH+1 ) , -VDS I N( JH+ 1 > )-DCMPLX< AV( OH > , -BV< JH ) ) 533. 

C FIND THE CHANGE IN THE REFLECTED VOLTAGE SINCE THE LAST CYCLE 534. 

DVRIOH KVX-VRvOH ) 53S. 

VR<OH)“VX 536. 

C— CALCULATE THE IMPEDANCE RATIOS AT EACH LO HARMONIC 10 DETERMINE 537. 

C — -THE DEGREE OF CONVERGENCE 538. 

ZQ=ZD/2EMB( OH ' 539. 

ZQMAGI OH ) a CDABS( ZU ) 540. 

ZQPHAf OH )“DATAN2( D IMAGI ZQ ),DREAL(ZQ) )*57. 29577951 D0 541 . 

IFIOH.EQ.l) GOTO 7 542. 

I F < ZQMAGI OH ) . GT. 1 .0D0+2QACC ) ZQFLAG n ZQFLAG+ 1 543. 

IF(ZQMAG<OH).LT.l .0D0-ZQACO ZQFLAG=ZQFLAG+ 1 544. 

7 CONTINUE 545. 

C THE DC REFLECTED VOLTAGE AND ITS CHANGE SINCE THE LAST CYCLE 546. 

VXDC=0. 5P0*t VDCOSI 1 >-Z0*IDCOS< 1 ) > 547. 

DVRDC«VXDC-VRDC 548. 

VRDC “VXDC 549. 

IFIOPRINT.NE.SO GOTO 9 550. 

C— CALL PRINT2 TO WRITE THE RESULTS OF THIS REFLECTION CYCLE 551. 

IFIOVLO.NE. IVLO) GOTO 9 552. 

CALL PRINT21RH0.VI , VR , DVR , VDCOS , VOS I N , I DCOS, I DS IN , ZQMAG , 2QPHA, 553. 

1VRDC . DVRDC , VIDC , RHODC ,AV , BV, ITER ,ZQF LAG, OVLO, NH , NHP 1> 554. 

9 CONTINUE 555. 

C F I ND THE NEW VOLTAGE INCIDENT ON THE DIODE AT THE LO AND HARMONICS 556. 

DO 10 0H= 1 , NH 557. 

VX°VI ( OH )+RHO( OH )*DVR1 OH > 558. 

VI < OH }=VX 559. 

AVION >«DREAL< VX > 560. 

10 8V< OH ) n -D IMAGt VX ) 561 . 

C-— FIND THE NEW DC VOLTAGE WAVE INCIDENT ON THE DIODE 562. 

VIDC=VIDC+RHODC*DVRDC 563. 

C DON ' T ADOUST THE DC CURRENT UNTIL WE HAVE RUN FOR ENOUGH CYCLES TO 564. 

C REACH A STEADY STATE 565. 

IF! ITER. NE. 10) GOTO 11 566. 

C ADOUST THE DC CURRENT TO THE DESIRED VALUE BY CHANGING VLO 567. 

CALL ADOLOIOVLO, IVLO, VLO, IDCOS, IDBIAS.NNPl > 568, 

C WAS THIS THE LAST VLO ADOUSTMENT LOOP? 569. 

I F < OVLO . EQ . IVLO ) GOTO 1 1 570, 

C REPEAT THE ANALYSIS WITH A NEW VALUE OF VLO 571. 

OVLO»OVLO+ 1 572. 

GOTO 15 573. 

C WAS THIS THE LAST REFECTION CYCLE ALLOWED? 574. 

11 IF( ITER. EQ. NITER > GOTO 12 575. 

C HAS THE SOLUTION CONVERGED? 576. 

IFIZOFLAG.EQ.0) GOTO 12 577. 

C GO ON TO THE NEXT REFLECTION CYCLE 578. 

GOTO 3 579. 

C CALL PRINT2 TO WRITE THE RESULTS OF THE FINAL REFLECTION CYCLE 590. 
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■12 CALL PRINT2C RHO, VI , VR, DVR, VDCOS * VDSIN , IDCOS, IDS1N*ZQMAG,ZQPHA, 

1 VRDC f DVRDC , VIDC , RHODC, AV , BV, ITER,ZQFLAG , JVLO, NH , NHP 1 > 

C CALL POWER TO FIND THE REQUIRED La POWER 

CALL POWERCZQMAGC 1 ),ZQPHAC 1 >,ZER{ 1),ZEI( 1 >,RSLO( 1 ) t VLO,Z0,LOPWR> 
C -UNSCALE THE CAPACITANCE VALUES (THEY WERE SCALED IN SUBROUTINE FCT 

C WHICH IS CALLED BY THE DRKGS INTEGRATION ROUTINE ) . 

DO 19 0PT=1 r NPTS 

1 9 C\1DATA( OPT >=C0 DATA! OPT ) / C 2 .0D0*P I*FP i 
Cr — -FINISH THE ANALYSIS BY OBTAINING THE FOURIER COEFFICIENTS OF THE 

C D I0 DE CONDUCTANCE AND CAPACITANCE, 

CALL DFORITC GODATA, NPTS/2, NH , GJCOS , GJSIN , IER ) 

‘CALL D FOR I T< CO DATA, NPTS/2, NH ♦ COCOS, CJS IN » IER ) 

C CALL PLOT TO PRINT THE DIODE WAVEFORMS IN THE TIME DOMAIN 
CALL PLOT t IGJDAT , CO DATA, VD DATA , ID DATA, NPTS , ITER,C0 ) 

RETURN 

END 


SUBROUTINE ZEMBEDC ZER, ZE I , ZEMBDC , ZEMBSB , RS , RSLO, RSSB , RSKIN , 

1FP , I F , P I , NH ,NHP 1 , NHD2P1 ) 

C 

C ZEMBED FORMS THE EMBEDDING IMPEDANCES AT THE HARMONICS OF THE 
C LO “AND AT THE. SIDEBAND FREQUENCIES (ASSUMING THEY HAVE NOT BEEN 
C INPUT VIA THE BLOCK DATA PROGRAM). IT ALSO FORMS THE SERIES 
C RESISTANCE, INCLUDING SKIN EFFECT, AT THESE FREQUENCIES. 

C NOTE THAT IF YOU WISH TO INPUT THE SIDEBAND EMBEDDING IMPEDAN- 
C CES THROUGH THE BLOCK DATA SUBPROGRAM THE SIDEBAND FREQUENCY NO- 
C TATION MUST BE TAKEN INTO ACCOUNT. ALL LOWER SIDEBAND EMBEDDING 
C IMPEDANCES (ZEMBSBCI) , I . GT. CNH/2 ) + l > SHOULD BE FORMED AS THE 
C COMPLEX CONOUGATES OF THEIR POSITIVE FREQUENCY VALUES. THIS IS 
C CONSISTANT WITH THE USE OF NEGATIVE FREQUENCIES FOR ALL LOWER 
C SIDEBANDS. NOTE THAT SIDEBAND I IS ARRAY ELEMENT CNH/2 + 1-1) 

C IN THIS FREQUENCY NOTATION, 

C 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 

COMPLEX*! 6 ZEMBSBC NHP 1 ) 

REAL*8 ZERC NH),ZEI{ NH), ZEMBDC ,RS, RSKIN, FP , IF, PI 
REAL*8 RSSB{NHPn,RSLO(NH} 

INTEGER NH , NHP 1 , NHD2P 1 , K, I 

C™ IN THIS EXAMPLE THE EMBEDDING IMPEDANCES AT THE LO HARMONICS 

C AND AT THE SIDEBAND FREQUENCIES HAVE BEEN INPUT VIA THE BLOCK DATA 

C SUBPROGRAM AND THUS THEY WILL NOT BE FORMED IN THIS SUBROUTINE 

C 

C FORM THE SERIES RESISTANCE OF THE DIODE AT THE LO HARMONICS AND 

C THE SIDEBAND FREQUEr.’' ,ES. THE FREQUENCY DEPENDENCE DUE TO THE 

C SKIN EFFECT IS ASS.^.D TO HAVE THE FORM RS< F )=RS+RSKI N*SQRT< F ) 

DO 40 1 = 1 ,NHP 1 
K=NHD2P1-I 

RSSB C I )=R5+RSKI N*DSQRT( DAB SC FP*K+IF ) ) 

I F ( I .EQ.NHP1 ) GOTO 40 
RSLOC I )=RS+RSKIN*DSQRT( I*FP > 

40 CONTINUE 
RETURN 
END 
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639. 

SUBROUTINE POWER! ZQMAG .ZQPHA, ZER.ZEI , RSLO ,VLO ,Z0, LOPWR ) 64jEf . 

641. 

POWER CALCULATES AND PRINTS THE REQUIRED LO POWER USED BY THE 642. 

MIXER WITH THE ORIGINAL VALUE OF ZE(1) C BEFORE IT WAS SET TO 213), 643. 

644. 

—THE VARIABLE TYPES USED'IN THIS SUBROUTINE ARE AS FOLLOWS: 645. 

COMPLEX*I6 ZQ,ZE 646. 

REAL*8 ZER.ZEI, ZQMAG, ZQPHA, RSLO, VLO,Z0, LOPWR, D 647. 

C-— TRANSFORM ZQMAG AND ZQPHA INTO THE REAL AND IMAGINARY PARTS 64B. 

C OF A COMPLEX NUMBER 649. 

D=DSQRT! 1+! DTANt ZQPHA/ 5 7 . 2957795 1D0) )**2 ) 650. 

ZQ=DCMPLX( ZQMAG /D , ZQMAG *DTAN( ZQP HA/57 .29577951 D0)/D ) 651 . 

ZE=DCMPLX( ZER,ZE I ) 652. 

C CALCULATE THE LO POWER 653. 

LOPWR=!CDABS!VLO*{.RSLO+ZE+ZQ*Z0)/iZ0+Z0*ZQ) )**2>/( B.0D0*ZER) 654. 

C— PRINT THE RESULTS 655. 

WRITE! 6 , 100) LOPWR 656. 

IM FORMAT(//2X, 'REQUIRED LO POWER : ' , 1PE 10. 3// ) 657. 

RETURN 658. 

END 659. 

660. 

661 . 

662. 

SUBROUTINE ADJLO! OVLO , I VLO , VLO , I DCOS , IDB IAS , NHP 1 ) 663. 

C 664. 

C SUBROUTINE ADJLO ADJUSTS THE LOCAL OSCILLATOR VOLTAGE UNTIL 665. 

C— -THE RECTIFIED CURRENT IS WITHIN IDCACC OF THE DESIRED VALUE , IDBIAS. 666. 

C 667. 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 668. 

C— FOR COMMON/ADJVLO/: 669. 

REAL*8 LOVLO.UPVLO.VLOINC, IDCACC 670. 

INTEGER UPFLAG , LOFLAG 671. 

C FOR VARIABLES NOT IN ANY COMMON BLOCKS: 672. 

REAL*8 IDCOSINHP1 ), VLO, IDBIAS 673. 

INTEGER OVLO, IVLO.NHPl 674. 

C THE COMMON BLOCKS USED ARE: 675. 

COMMON/ADJVLO/LOVLO, UPVLO, LOFLAG, UPFLAG, VLOINC, IDCACC 676. 

C IF THIS IS ALREADY THE LAST VLO LOOP THEN DON'T OUTPUT 677. 

IF(OVLO.EQ. IVLO) GOTO 25 678. 

WRITE! 6 , 100) JVLO, IDBIAS , 679. 

100 FORMAT! /’ VALUES OF THE DC CURRENT AND LO VOLTAGE FOR CYCLE 680. 

112,' OF THE LOOP FOR ADJUSTING VLO TO G IVE ',F8.6,' AMPS ARE:') 681. 

WRITF! 6 , 1 10) IDCOS! 1 ) , VLO 682. 

110 FORMAT! 10X,' IDCOSf 1 >=' ,F8.6,T35, ' VLO BEFORE ADJUSTMENT: ', FB .5 ) 683. 

IF! IDCOS! 1 ) ,GT. IDBIAS+! I DB IAS*IDCACC > > GOTO 10 684. 

IF ! IDCOS! 1 ).LT. IDBIAS-! IDB IAS* IDCACC ) > GOTO 15 685. 

I VLO=JVLO 686. 

GOTO 20 687. 

10 UP VLO=VLO 688. 

C KEEP TRACK OF THE NUMBER OF TIMES VLO IS GREATER THAN ITS DESIRED 689. 

C— VALUE 690. 

UPFLAG=UPFLAG+1 691. 

C IF WE HAVE NOT YET PASSED THE DESIRED VLO CHANGE VLO 692. 

IF! LOFLAG. EQ.0) GOTO 11 693. 

VLO=VLO-( UPVLO-LOVLO1/2.0D0 694. 

GOTO 20 695. 

11 VLO=VLO-VLOINC 696. 
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IFiVLO.LT .0.0) VLO=0.0 Dj 0 697. 

GOTO 20 698. 

IB LOVLQ=VLO 699. 

C KEEPING TRACK OF THE NUMBER OF TIMES VLO IS LESS THAN ITS DESIRED 700. 

C- — VALUE 701. 

LOFLAG=LOFLAG+I 702. 

C- IF WE HAVE NOT YET PASSED THE DESIRED VLO, CHANGE VLO 703. 

IF1UPFLAG.EO.0) GOTO 16 704. 

VLO=VLp+( UPVLO“LOVLO }/2 ,0D0 705. 

GOTO 20 706. 

16 VLO=VLO+VLOINC 707. 

20 WRITE! 6, 120) VLO 708. 

120 FORMAT! T35, ' VLO AFTER ADJUSTMENT: ' ,F8.5) 709. 

25 RETURN 710. 

END 711. 

712. 

713. 

714. 

SUBROUTINE DRKGS! PRMT, Y , DERY , ND IM« IHLF , FCT.OUTP ,AUX ) 715. 

C 716. 

C DRKGS IS AN IBM SSP PROGRAM WHICH SOLVES A SYSTEM OF DIFFERENTIAL 717. 

C— EQUATIONS BY THE RUNGE-KUTTA ALGORITHM. IT HAS NOT BEEN ALTERED 718. 

C— FOR THIS ANALYSIS. 719. 

C-— 720. 

DIMENSION Y! 1 ),DERYt 1 ) , AUX! 8, 1 ),A!4),B!4),C!4), PRMT! 1 •) 721. 

DOUBLE PRECISION PRMT, Y , DERY , AUX ,A, B , C , X , XEND , H ,AJ , BJ ,CO , RI ,RZ , 722. 

1DELT 723. 

DO 1 1=1 ,NDIM 724. 

1 AUX! 8 , I )= .066666666666666 667D0*DERY( I ) 725. 

X=PRMT! 1 ) 726. 

XEND=PRMT! 2 ) 727. 

H=PRMT! 3 ) 728. 

PRMT! 5 )=0. D0 729. 

CALL FCTC X, Y , DERY ) 730. 

IF!H*!XEND-X) )38 , 37 , 2 731. 

2 A! 1 )= . 5D0 732. 

A! 2 )=. 292 B 9 32 188 13 4524 8 D0 733. 

A! 3 )= 1 . 707 106781 1865475D0 734. 

A! 4 )=. 16666666666666667 D0 735. 

B! 1 >=2 . D0 736. 

B!2)=1.D0 737. 

B!3)=1.D0 738. 

B! 4 )=2 . D0 739. 

C! 1 )=. 5D0 740. 

C! 2 >=.29289321881 345248 D0 741. 

C( 3 )=1 . 707 106781 1865475D0 742. 

C! 4 )=. 5D0 743. 

DO 3 1=1 , ND IM 744. 

AUX! 1 , 1 )=Yt I ) 745. 

AUX! 2 , I >=DERYC I > 746. 

AUX! 3,1 )=0. D0 747. 

3 AUXC6, I )=0.D0 748. 

IREC=0 749. 

H=H+H 750. 

IHLF=-1 751. 

I STEP=0 752. 

IEND=0 753. 

4 IF! !X+H-XEND)*H)7,6,5 754. 
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5 H=XEND-X 

6 IEND=i 

7 CALL OUTP £ X, Y* DERY, IREC, NDIM, PRMT ) 
I P ( P RMT C 5 ) ) 40 * 9 , 40 

8 ITEST=0 

9 ISTEP= ISTEP+ 1 
0=1 

10 AO=A( 0 ) 

BO = B t 0 ) 

C0=C<0) 

DO 11 1=1 , NDIM 
R1=H*DERY ( I ) 

R2=AJ*(Rl-B0*AUX<6 t I)) 

Y ( I )=Y( I )+R2 
R2= R2 + R2+R2 

11 AUXt6,l)=AUX<6, I)+R2-CJ*R1 
IP 1 0-4 >12,15,15 

12 0 = 0+1 

I F 1 0-3 >13,14,13 

13 X=X+,5D0*H 

14 CALL FCTC X,Y,DERY) 

GOTO 10 

15 IF C ITESTJ16, 16,20 

16 DO 17 1=1, NDIM 

17 AUXU, I )=Y( I ) 

ITE ST= 1 

I STEP = I STEP + 1 STEP -2 

18 IHLF=IHLF+1 
X = X-H 

H= . 5D0*H 
DO 19 1 = 1, NDIM 
Y( I }=AUX( 1,1) 

DERYt I }=AUXC2,I) 

19 AUXt 6 , 1 )=AUX C 3,1) 

GOTO 9 

20 I MOD= I STEP/2 

IP( I STEP “ IMOD- IMOD >21,23,21 

21 CALL FCTt X, Y , DERY ) 

DO 22 1=1, NDIM 
AUXt 5, I ) = Y( I ) 

22 AUX(7, I ) = DERY t I > 

GOTO 9 

23 DELT=0.D0 

DO 24 1 = 1, NDIM 

24 DELT=DELT+AUX{ 8, I >*DABSC AUXt 4 , I )~Y 
IF ( DELT-PRMTt 4 ) >28,28,25 

25 IF { IHLF-10J26, 36,36 

26 DO 27 1=1 , NDIM 

27 AUX ( 4 , I ) =AUX( 5,1) 

I5TEP = ISTEP + I STE P-4 
X = X“H 

I E ND=0 
GOTO 18 

28 CALL FCTt X , Y , DERY ) 

DO 29 1=1, NDIM 
AUX(1,I> = YU) 

AUXt 2 , I > = DERY ( I ) 

AUX ( 3 , I )=AUXt 6,1) 
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Y<I)=AUX<5,I) 813 . 

29 DERYC I )=AUX£ 7,1) 814. 

CALL aUTP£X-H,Y,DERY,IHLF,NDIM,PRMT> 815. 

I F £ P RMT£ 5 ) ) 40 , 30 , 40 816. 

30 DO 31 1=1 ,NDIM 817. 

YtI)=AUX{I,n 818. 

31 DERYt I )=AUX(2,I) 819. 

IREC=IHLF 820, 

IF ( IEND )32 ,32 |39 821. 

32 IHLF=IHLF“1 822. 

ISTEP= ISTEP/2 823. 

H=H+H 824. 

IF £ IHLF )4 , 33 , 33 825. 

33 IMOD=I STEP/2 826. 

IF( ISTEP-IMQD-IMOD 14,34,4 827. 

34 IF( DELT- .02D0*PRMTC 4 > )35 ,35,4 828. 

35 IHLF=IHLF-1 829. 

ISTEP=ISTEP/2 830. 

H=H+H 831. 

GOTO 4 832. 

36 IHLF=1 1 833. 

CALL FCT( X , V ,DERYJ 834. 

GOTO 39 835. 

37 IHLF=12 836. 

GOTO 39 837. 

38 IHLF=13 838. 

39 CALL GUTP( X , Y »DERY , IHLF , NDIM.PRMT) 839. 

40 RETURN 840. 

END 841. 

842. 

843. 

844. 

SUBROUTINE FCT( X , Y , DERY ) 845. 

C 846. 

C FCT IS REQUIRED BY DR KGS AND SUPPLIES THE NETWORK STATE EQUATION TO 847. 

C BE SOLVED. NOTE THAT THE JUNCTION CAPACITANCE HAS BEEN FREQUENCY 848. 

C SCALED BY 2*PI*FP SO THAT ONE LO CYCLE OCCURS IN 2*PI SECONDS 849. 

C 850. 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 851. 

C FOR COMMON/CONST/ : 852. 

REAL*8 QEL , BOLTZ , P I ,TEQ 853. 

C FOR COMMON/DIODE/: 854. 

REAL*8 ALP , PHI , GAM ,C0, IS , RS , RSKIN , FP , IF , IGJ , ICO ,GJ , CJ 855. 

C FOR COMMON/LOOPS/: 856. 

INTEGER NH.NLO.NVLO.NPTS.NPRINT, NITER 857. 

C FOR COMMON/TLINE/: 858. 

REAL *8 Z0,ZQACC 859. 

INTEGER ZQFLAG 860. 

C FOR COMMON/VOLTS/: 861. 

REAL*8 AV( 8 > ,BV( 8 ) , VI DC , VLO , VDB IAS , IDBIAS 862. 

C— F oR VARIABLES NOT IN ANY COMMON BLOCKS; 863. 

REAL*8 X , Y( 1 ) , DERYt 1 ) , V 864. 

INTEGER JH 865. 

C THE COMMON BLOCKS USED ARE: 866. 

COMMON/CONST/QEL .BOLTZ, PI, TEQ 867. 

COMMON/D IODE/ALP, PHI, GAM, C0, IS, RS, RSKIN, FP, IF, IGJ, ICO, GO, CJ 868. 

COMMON/ LOOP S /NH, NL 0, N VLO , NPTS, NPRINT.NITER 869. 

COMMON /TL I NE /Z0, ZQACC .ZQFLAG 870. 
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COMMON/VOLTS/AV, BV,VIDC,VLO , VDBIAS , IDBIAS B7I. 

C— CALCULATE THE TOTAL VOLTAGE ON THE TRANSMISSION LINE INCIDENT ON 872. 

C THE DIODE. 873. 

V=VIDC 874. 

DO 1 JH=1 , NH 87S , 

I V=V+AV< JH)*DCOSCUH*XHBVCOH)*DSIN{JH*X> 876. 

C— MULTIPLY BY 2 TO FIND THE EQUIVALENT SOURCE VOLTAGE ASSOCIATED 877. 

C— WITH THE TRANSMISSION LINE IMPEDANCE Z0 878. 

V=V*2.0D0 879. 

C FIND THE FREQUENCY SCALED JUNCTION CAPACITANCE. B80. 

C NOTE THAT THE DIODE VOLTAGE Ytl) IS ARTIFICIALLY CLAMPED IF IT 881. 

C— EXCEEDS . 999*PNI . THIS AVOIDS ERROR MESSAGES DURING THE INITIAL 882. 

C— TRANSIENT IN THE NUMERICAL SOLUTION. 883. 

CO =2 . 0D0*P I*FP*C0/( 1 .0D0-DMI N1 { Y C 1 >,0.999*PHI J/PHI >**GAM 884. 

C— FIND THE CURRENT THROUGH THf. DIODE CONDUCTANCE 885. 

IG0 = IS* t DEXPl ALP*DMAX 1( -174 . J0TDJ3 1 / AL P , DM INKY! 1) , 174 .0D0/ALP ) ) ) 886. 

1-1. 0D0) 887. 

C DVD/DT 888. 

DERY ( 1 )=( { V-Y{ 1 >>/Z0-IGJ)/CJ 889. 

C FIND THE CURRENT THROUGH THE DIODE CAPACITANCE AND THE DIODE 890. 

C— -CONDUCTANCE 891. 

ICJ=DERY( 1 )*C0 892. 

GU=IGJ*ALP+IS*ALP 893. 

RETURN 894. 

END 895. 

89G . 
897. 
89 8 * 

SUBROUTINE OUTP ( X, Y ,DERY , IHLF , ND IM, PRMT ) 899! 

900 . 

—OUTP IS REQUIRED BY DRKGS BUT IS NOT USED IN THIS PROGRAM 9J0T1 . 

902. 

REAL*8 X,Y(1 J.DERYU >,PRMT(5> 903. 

INTEGER IHLF.NDIM 904. 

RETURN 905. 

END 906. 

907. 

908. 

909. 

SUBROUTINE DFORITt FNT.N ,M,A, B , IER ) 910. 

911. 

DFORIT IS A DOUBLE PRECISION VERSION OF FORIT.AN IBM SSP ROUTINE 912. 

—THAT PERFORMS A FOURIER ANALYSIS. ON A PERIODICALLY VARYING FUNCTION. 913. 
-“-IT COMPUTES THE COEFFICIENTS OF THE TERMS IN THE SERIES WHICH 914. 

— IS GIVEN BY :A( 1 )+SUM( A( N )COS< {N-1)X) + B(N)SIN((N“1 >X)> N=2,3,4... 915. 

916. 

—THE PARAMETERS USED ARE: 917. 

FNT/: TABULATED VALUES OF THE FUNCTION TO BE ANALYSED 918. 

NOTE THAT FNT( I ) CORRESPONDS TO TIME T=0 919. 

— M/: THE MAXIMUM ORDER OF THE HARMONICS TO BE FITTED 920. 

— N/s DEFINES THE INTERVAL OVER WHICH THE POINTS ARE TAKEN. THE 921. 

INTERVAL GOES FROM 0 TO 2*PI AND 2N+1 POINTS ARE TAKEN AS DATA. 922.. 
—A/: THE FOURIER COSINE COEFFICIENTS 923. 

B/: THE FOURIER SINE COEFFICIENTS 924. 

— IER/: THE' RESULTANT ERROR MESSAGE CODE WHERE IER=0 MEANS NO ERROR, 925. 

IER=1 MEANS N IS LESS THAN M, AND IER=2 MEANS M IS LESS THAN 0 926. 

] 927 

REAL*8 A< 1 ) , Bt 1 } , FNTt 1 >-,CONST 928 ! 
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REALMS COEF , C,S,C1 t Sl , AN # FNTZ,U0, U1 , US,G 929, 

INTEGER N ,M 930. 

IER-0 931 , 

20 I F ( M > 30,40, 40 932. 

30 IER=2 9 33 , 

RETURN 934 . 

40 IF(M-N) 60,60,50 935. 

50 I E R— X 936. 

RETURN 937, 

60 AN=N 938. 

COEF=2.0D0/(2.0D0*AN+1 ,0D0> 939. 

CONST-3* 1 4 1 5926535B979D0*CQEF 940. 

S1=DSI N( CONST ) 941 . 

Cl =DCOS( CONST } 542. 

C=1.0D0 943. 

S-0. 0D0 944. 

J = 1 945. 

FNTZ=FNT( 1 > 946. 

70 U2=0 . 0D0 947. 

U 1 -0 .0D0 948. 

I=2*N+ 1 949. 

75 U0=FNT( I )*2 .0D0*C*U1-U2 950, 

U2=U 1 951. 

U1=U0 952. 

1=1-1 953. 

I F ( 1-1 > 90,80,75 954. 

80 A(0)=COEF*CFNTZ+C*U1-U2> 955, 

B(D )=C0EF*S*U1 956. 

IF{0“*(M+1 > ) 90,100,100 957. 

90 G=C1*C-S1*S 958. 

S=C1*S+SI*C 959. 

C=Q 960. 

+ 1 961. 

GO TO 70 962. 

100 A( 1 )=A( 1 ) *0 . 5D0 963. 

RETURN 964. 

END 965. 

966 . 

967. 

968, 

SUBROUTINE SMSIG 969. 

C 970* 

C SMSIG PERFORMS THE SMALL-SIGNAL AND NOISE ANALYSIS OF THE MIXER 971. 

C TO DETERMINE THE CONVERSION LOSS BETWEEN ALL PAIRS OF SIDEBANDS, 972. 

C THE INPUT AND OUTPUT IMPEDANCES , AND THE EQUIVALENT INPUT NOISE 973. 

C TEMPERATURE. THE THEORY IS BASED ON THAT OF D * N * HELD AND A. R. KERR, 974. 

C IEEE TRANS. ON MICROWAVE THEORY AND TECH . , VOL . MTT-26 , NO . 2 , PP . 49-6 1 , 975. 

C FEB. 1978, 976. 

C ALL INITIALIZED VARIABLES ARE INPUT THROUGH THE BLOCK DATA SUB- 977. 

C PROGRAM AND ARE TRANSFERRED THROUGH COMMON BLOCKS. 978. 

C THE OUTPUT INCLUDES: 979. 

C 1) THE CONVERSION LOSS BETWEEN ALL PAIRS OF SIDEBANDS (PRINTED 990, 

C AS A CONVERSION LOSS MATRIX). 981. 

C 2) THE INPUT IMPEDANCES OF THE MIXER AT EACH SIDEBAND* 982. 

C 3) THE OUTPUT IMPEDANCE AT THE IF. 983. 

C 4) THE EQUIVALENT INPUT NOISE TEMPERATURE AND ITS SHOT AND THERMAL 984. 

C COMPONENTS. 985. 

C THE SUBSCRIPT NOTATION USED IN THE PROGRAM TO IDENTIFY THE NH + 1 986. 
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C SMALL-SIGNAL SIDEBANDS IS THAT OF A.A,M. SALEH, 'THEORY OF RESISTIVE 987. 

C MIXERS' ,M. I. T. PRESS,CAMBRIDGE ,MASS . , 1971 . SIDEBAND FREQUENCY 988. 

C ( IF+N*LO ) IS DENOTED BY THE ARRAY SUBSCRIPT ! NH/2 + 1 - N). THE 989. 

C LOWER SIDEBANDS ARE TREATED AS NEGATIVE FREQUENCIES CONSIDERABLY 990. 

C SIMPLIFYING THE EQUATIONS IN THE ANALYSIS. 991. 

C IF ARRAY DIMENSIONS ARE ALTERED. THEY MUST BE CHANGED HERE, IN 992. 

C SUBROUTINE LGSIG AND IN THE BLOCK DATA PROGRAM, IN ADDITION THE 993. 

C PRINT FORMAT OF THE CONVERSION LOSS MATRIX MUST BE ALTERED IF 994. 

C OTHER THAN EIGHT LO HARMONICS’ ARE CONSIDERED. 995. 

C 996. 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 997. 

C FOR COMMON/CONST/: 998, 

REAL*8 QEL,BOLTZ,PI ,TEQ 999. 

C— -FOR COMMON/DIODE/: 1000. 

REAL *8 ALP , PH I ,GAM , C0, IS , RS , RSKIN , FP , IF , IGO , ICO , GO , CO 1001. 

C— -FOR COMMON/FORITS/: 1002. 

REALMS GOCQS! 9 ) ,GOSIN( 9 ) ,COCOSC 9 ) ,COS IN! 9 ) , VDCOS! 9 ) , VDSIN! 9 ) 1003. 

REAL*8 IDCOS< 9),IDSIN(9) 1004. 

INTEGER IER 1005. 

C FOR COMMON/ IMPED/: 1006. 

COMPLEX*16 ZEMBSB! 9) 1007. 

REAL*8 ZER( 8 ) ,ZE I ( 8 ) , ZEMBDC ,RSLO(8),RSSB(9) 1008. 

C— FOR COMMON/LOOPS/: 1009. 

INTEGER NH ,NLO,NVLO, NPTS ,NPRI NT, NITER 1010. 

C— FOR VARIABLES NOT IN ANY COMMON BLOCKS: 1011. 

COMPLEX*l 6 Ai9,9),COR<9,9>,FG!’9>,FC<9) 1012. 

C0MPLEX*16 T!9),ZIN!9), ZIFOUT , DET 1013. 

REAL*8 XLMAT!9, 9), TM, THERM, SHOT, LIO 1014. 

REAL*8 GO MAG ( 9 ) ,GOPHA! 9 ) , COMAG! 9 ) , C JPHA! 9 ) 1015. 

INTEGER 0H,NHPI,NHD2P1 , NHD2,WK1< 9 ) , WK2! 9 ) , 1 ,0 1016. 

C— THE COMMON BLOCKS USED ARE: 1017. 

COMMON/CON ST/QEL , BOLTZ, P I , TEQ 1018. 

COMMON/ D I ODE /ALP , PHI , GAM,C0, IS, RS , RSKIN , FP , IF , IGJ , ICO ,G0 ,C0 1019. 

COMMON/FORITS/GOCOS ,GJSIN ,C0 COS, COSIN, VDCOS,VDSIN,IDCGS,IDSIN, IER 1020. 
COMMON/ IMPED/ZER.ZE I, ZEMBDC, ZEMBSB , RSLO, RSSB 1021 . 

COMMON/ I VMAG/GOMAG,GOPHA, COMAG, COPHA 1022. 

COMMON/LOOPS/NH,NLO,NVLO, NPTS, NPRI NT, NITER 1023. 

C DEFINE SOME USEFUL CONSTANTS 1024. 

NHP 1=NH+ 1 1025. 

NHD2=NH/2 1026. 

NHD2P1=NHD2+1 1027. 

C FORM THE COMPLEX FOURIER COEFFICIENTS OF THE DIODE CONDUCTANCE 1028. 

C AND CAPACITANCE 1029. 

DO 10 0H=2,NHP1 1030. 

FG( OH )=DCMPLX! GOCOS(OH),-GJSIN!OH > )*0.5D0 1031. 

10 FC!OH)=DCMPLX!COCOS!OH),-COSIN!OH ) >*0.5D0 1032. 

FGt 1 >=DCMPLX!GOCOS< 1 ) ,0.000) 1033. 

FCt 1 >=DCMPLX( COCOS! 1 ),0.0D0) 1034. 

C CALL PRINT3 TO WRITE THE FOURIER COEFFICIENTS 1035. 

CALL PRI NT3( FG ,FC , GO MAG , GJPHA, COMAG ,C0PHA,NH,NHP1 ) 1036. 

c OPEN CIRCUIT THE IF LOAD TO FIND THE IF PORT IMPEDANCE 1037. 

ZEMBSB! NHD2P1 >=DCMPLX( 1 .0D10.0.0D0) 1038. 

C FORM THE V* MATRIX WITH THE OPEN CIRCUITED IF BY CALLING YPRIME 1039. 

CALL YPRIME!FG,FC,NHD2,NHD2P1,NHP1 ,FP, IF, A, ZEMBSB, RSSB) 1040. 

C— TAKE THE INVERSE OF THE Y' MATRIX TO FIND THE OUTPUT IMPEDANCE 1041. 

CALL CMINV! A,NHP 1 , DET, WK1 , WK2 , NHP1*NHP 1 > 1042. 

C THE IF OUTPUT IMPEDANCE IS THE CENTER ELEMENT OF THE Z' MATRIX+RS 1043. 

ZIFOUT=A! NHD2P1 ,NHD2P1 )+RSSB(NHD2Pl ) 1044. 
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C— CONJUGATE MATCH THE IF LOAD IMPEDANCE TO THE IF PORT IMPEDANCE 1045. 

ZEMBSBt NHD2P1 )=DCONJG< ZIFOUT > 1046. 

C— FORM THE Y* MATRIX WITH A MATCHED IF LOAD 1047. 

CALL YPRIME(FG,FC,NHD2,NHD2P1,NHP1,FP,IF,A,ZEMBSB,RSSB) 1048. 

C INVERT THE V’ MATRIX TO OBTAIN THE Z’ MATRIX 1049. 

CALL CMINV(A,NHP1)DET,WK1,WK2,NHP1*NHPI) 1050. 

C FORM THE CONVERSION LOSS MATRIX AND INPUT IMPEDANCE AT EACH SIDEBAND1051 . 

DO 50 1=1 ,NHPl 1052. 

ZINt I )=RSSB( I )+A< I , I )*t RSSBt I )+ZEMBSB{ I ) >/<RSSB( I )+ZEMBSB( I ) 1053. 

l-AM.in 1054. 

DO 40 J = 1 , NHP 1 1055. 

IFtl-J) 20.30,20 1056. 

20 LI0 = t f CDABSt RSSBt I I+ZEMBSBl I ) >*CDABSt RSSBt J )+ZEMBSBtO ) )/ 1057. 

1 ( 2 ,0D0*CDABS( A( 1,0))) )**2 )/( DREALt ZEMBSBt I ) )*DREAL t ZEMBSBt J )) > 1058. 

C— CONVERT TO DB WHEN FORMING THE LOSS MATRIX 1059. 

XLMATt I ,0 ) = 10.0D0*DLOG10( LIO ) 1060. 

GOTO 40 1061 . 

C THE DIAGONAL ELEMENTS HAVE NO OBVIOUS MEANING AND ARE ZEROED FOR 1062. 

C CONVENIENCE 1063, 

30 XLMATt 1,0 )=0.0D0 1064. 

40 CONTINUE 1065, 

50 CONTINUE 1066. 

C BEGIN THE NOISE ANALYSIS BY FORMING THE SHOT NOISE CORRELATION 1067. 

C— MATRIX 1068. 

CALL CORREL t ALP , FG ,COR ,NHP 1 ) 1069. 

C— -CALCULATE THE EQUIVALENT INPUT NOISE TEMPERATURE CONSIDERING THE 1070. 

C SHOT NOISE ONLY 1071. 

CALL TMIXt NHD2 , SHOT, T , COR , A , RSSB,ZEMBSB,NHP1 ,NHD2P1 ) 1072. 

C ADO THE THERMAL NOISE COMPONENT TO THE SHOT NOISE CORRELATION MATRIX1073. 

CALL TNOI SE ( COR, RSSB , ZEMBSB , NHP 1 .NHD2P1 ) 1074. 

C— -RECALCULATE THE EQUIVALENT NOISE TEMPERATURE NOW INCLUDING BOTH 1075. 

C— SHOT AND THERMAL NOISE CONTRIBUTIONS 1076. 

CALL TMIXtNHD2,TM,T, COR, A, RSSB, ZEMBSB, NHP1 ,NHD2P1 ) 1077. 

C— FIND THE THERMAL NOISE COMPONENT OF THE EQUIVALENT INPUT NOISE 1078. 

C— TEMPERATURE OF THE MIXER 1079. 

THERM=TM-SHOT 1080. 

C-— CALL PRINT4 TO PRINT THE RESULTS OF THE SMALL-SIGNAL AND NOISE 1081. 

C ANALYSIS 1082. 

CALL PRI NT4tXLMAT,ZIN,ZIF0UT,TM, THERM, SHOT. NHD2, NHP 1.NHD2P1 ) 1083. 

RETURN 1084. 

END 1085. 

1086. 

1087. 

1088. 

SUBROUTINE YPRIME t FG , FC , NHD2 , NHD2P 1 , NHP 1 , FP , I F , A , ZEMBSB , RSSB ) 1089. 

1090. 

YPRIME FORMS THE AUGMENTED ADMITTANCE MATRIX Y 1 <A> OF THE MIXER 1091. 

1092. 

THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1093. 

- — FOR COMMON/CONST/: 1094. 

REAL»8 QEL ,BOLTZ,PI ,TEQ 1095. 

C FOR VARIABLES NOT IN ANY COMMON BLOCKS: 1096. 

COMPLEX*! 6 ZEMBSBt NHP1 ), At NHP1 ,NHP1 ),FGtNHPl) ,FCt NHP1 > 1097. 

REAL *8 FP , IF , WP , WI F , RSSB t NHP 1 ) 1098. 

INTEGER NHD2,NHD2P1,NHP1,I,0 1099. 

C— THE COMMON BLOCKS USED ARE: 1100, 

COMMON/ CON ST/QEL , BOLTZ , P 1 , TEQ 1101 . 

WI F =2 .0D0*P I * I F 1102. 
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WP=2.0D0*PI*FP 1103. 

C FORM THE ADMITTANCE MATRIX Y OF THE INTRINSIC DIODF. 1104. 

DO 60 1=1 ,NHP1 1105. 

DO 50 0=1 ,NHP1 1105. 

20 IF(U-I) 30,40,40 1107. 

C— FIND THE LOWER HALF OF THE Y MATRIX 1108. 

30 A( 1 ,0 )=DCONOG(FG( 1-0 + 1 ) 1+DCMPLX(0.0D0,WIF+WP*{ NHD2P1-I > ) 1109. 

1*DCON0G(FC( 1-0+11) 1110. 

GOTO 50 1111. 

C FIND THF UPPER HALF OF THE Y MATRIX 1112, 

40 A( 1 ,0 )=FG(0-I + 1 >+DCMPLX<0.0D0,WIF+WP*CNHD2Pl-I >)*FC( 0-1 + 1 ) 1113. 

50 CONTINUE 1114. 

60 CONTINUE 1115. 

C— ADD l/( RS+ZEMBSB ) TO THE DIAGONAL ELEMENTS OF Y TO FORM THE 1116. 

C AUGMENTED ADMITTANCE MATRIX Y' OF THE MIXER 1117. 

DO 70 1=1 ,NHP 1 1118. 

70 A( I , I )=A< 1 , 1 1 + 1 .0D0/( ZEMBSBt I )+RSSB( 1 1 ) 1119. 

RETURN 1120. 

END 1121. 

1122 . 

1123. 

1124. 

SUBROUTINE CORRELl ALP , FG ,COR ,NHP 1) 1125. 

C 1126. 

C CORREL FORMS THE NOISE CURRENT CORRELATION MATRIX FOR THE SHOT 1127. 

C— NOISE. THE THERMAL NOISE COMPONENTS ARE ADDED IN SUBROUTINE 1128. 

C— TNOISE. 1129. 

C 1130. 

C— THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1131. 

COMPLEX*! 6 FGINHPl >,COR<NHPI ,NHPl ) 1132. 

REAL*8 QEL.BOLTZ.PI ,TEQ,ALP 1133. 

INTEGER NHP1,I,0 1134. 

C— THE COMMON BLOCKS USED ARE: 1135. 

COMMON/CONST/QEL , BOLTZ , P I , TEQ 1136. 

C— FORM THE SHOT NOISE CORRELATION MATRIX USING I=FG/ALP 1137. 

DO 10 1=1 „NHP1 1138. 

DO 20 0=1,1 1139. 

COR(0,T)=2.0D0*QEL*FG( I-J + l l/ALP 1140. 

20 COR( 1 ,0 >=DCON0G(COR(0, I ) ) 1141. 

10 CONTINUE 1142. 

RETURN 1143. 

END 1144. 

1145. 

1146. 

1147. 

SUBROUTINE TNOISECCOR.RSSB .ZEMBSB ,NHP1 ,NHD2P1 ) 1148. 

1149. 

-—TNOISE FORMS THE THERMAL NOISE CURRENT CORRELATION MATRIX AND 1150. 

—ADDS IT TO THE SHOT NOISE CORRELATION MATRIX. 1151. 

1152. 

THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1153. 

COMPLEX*16 COR( NHP 1 , NHP 1 > , ZEMBSBt NHP1 ) 1154. 

REAL *8 QEL.BOLTZ.PI , TEQ, RSSBINHPD 1155. 

C— -THE COMMON BLOCKS USED ARE: 1156. 

COMMON/CONST/QEL, BOLTZ.P I ,TEQ 1157. 

C— SINCE THE THERMAL NOISE MATRIX IS DIAGONAL ADD THESE ELEMENTS TO 1158. 

C THE DIAGONAL TERMS OF THE SHOT NOISE CORRELATION MATRIX 1159. 

DO 35 1=1 , NHP 1 1160. 
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IFU.EQ.NHD2P1 > GOTO 30 1161. 

C0R( I , I )=C0R( 1 , n+4 .0D0*BOLTZ*TEQ*RSSBi I >/ 1162. 

1( CDABSl ZEMBSBt 1 >+RSSB( I > ) )**2 1163. 

GOTO 35 1164. 

C AT THE IF THE THERMAL NOISE TERM IS GIVEN BY: 1165, 

30 COR( I ,!>=COR{ I , I )+4.0D0*BOLT2*TEQ*RSSB< I )/ 1166. 

1 ( CDABS( ZEMBSBC I >-RSSB( I ) ) 1**2 1167. 

35 CONTINUE 1168. 

RETURN 1169. 

END 1170. 

1171 . 

1172. 

1173. 

SUBROUTINE TMIXt NSB , TM , T.COR.A, RSSB , 2EMBSB ,NHP 1 , NHD2P 1 > 1174. 

C 1175. 

C TMIX COMPUTES THE EQUIVALENT SINGLE SIDEBAND INPUT NOISE 1176. 

C TEMPERATURE OF THE MIXER REFERRED TO SIDEBAND NSB. 1177. 

C— NOTE THAT SIDEBAND NSB IS ARRAY SUBSCRIPT NH/2 + 1 - NSB. 1178. 

C 1179. 

C— THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1180. 

COMPLEX* 16 COR{ NHP1 ,NHP1 1 , A( NHP 1 ,NHP 1 ) .ZEMBSBC NHP1 ) , VSQ, T( NHP I ) 1181 . 

REAL*8 QEL , BOLTZ ,PI ,TEQ,TM,RSSB<NHP1 ) 1182. 

INTEGER I,0,NSB,NHP1 .NHD2P1 1183. 

C— THE COMMON BLOCKS USED ARE: 1184. 

COMMON/CONST/QEL , BOLTZ , P I ,TEQ 1185. 

C— POST MULTIPLY COR BY THE CONJUGATE TRANSPOSE OF THE CENTER ROW 1186. 

C — -OF THE Z' MATRIX (ROW 0) 1187. 

DO 10 1=1 ,NHP 1 1188. 

T( I )=0.0D0 1189. 

DO 20 0=1, NHP 1 1190. 

20 T( I)=T( I l+CORU ,0 )*DCONOG(A(NHD2P1 ,0 )> 1191. 

10 CONTINUE 1192. 

C— “PREMULTIPLY COR BY THE CENTER ROW OF THE MIXER Z' MATRIX 1193. 

VSQ=0.0D0 1194. 

DO 30 1=1, NHP 1 1195. 

30 VSQ=VSQ+A(NHD2P1 ,1)*T< I) 1196. 

C— COMPUTE THE EQUIVALENT INPUT NOISE TEMPERATURE OF THE MIXER 1197. 

TM=( DREAL(VSQ)*( CDABS< ZEMBSBC NSB l+RSSBl NSB ) ) >**2 )/ 1198. 

1 ( BOLTZ*4 .0D0*DREAL< ZEMBSBC NSB ) )*( CDABS( A< NHD2P 1 ,NSB ))**2 ) ) 1199. 

RETURN 1200. 

END 1201. 

1202. 
1203. 
I 204 . 

SUBROUTINE CM I NV ( A , N , D , L , M , NSO ) 1205. 

1206, 

CMINV IS A SLIGHTLY MODIFIED VERSION OF THE IBM SSP ROUTINE MINV 1207. 

--FOR INVERTING A COMPLEX MATRIX, ONLY THE FIRST TWO STATEMENTS AND 1208. 

THOSE NUMBERED 10 AND 45 HAVE BEEN ALTERED. 1209. 

1210. 

CQMPLEX*16 A , D , B IGA , HOLD 1211. 

DIMENSION A< NSQ ) , L ( N > , Ml N ) 1212. 

0=1. 0D0 1213. 

N.K=“N - 1214. 

DO 80 K=1,N 1215. 

NK=NK+N 1216. 

L( K > = K 1217. 

M(K)=K 1218. 



KK=NK+K 

1215, 


BIGA-A(KK) 

1220. 


DO 20 0=K,N 

1221. 


IZ=N*(0-1 ) 

1222. 


DO 20 I=K, N 

1223. 


I0=IZ+I 

1224. 

10 

IFtCDABSt BIGAJ-CDABSCAC 10 ) ) } 15,20,20 

1225. 

15 

B IGA=A( Id ) * 

1226. 


L£K)=I 

1227. 


M( K )=0 

1228. 

20 

CONTINUE 

1229. 


0=L< K ) 

1230. 


IFtO-K) 35,35,25 

1231 , 

25 

KI=K-N 

1232. 


DO 30 1=1 , N 

1233. 


KI=KI 

1234. 


HOLD=-ACKI) 

1235. 


JI=KI-K+0 

1236. 


A(KI)=A{JI) 

1237, 

30 

At 0 1 ) =HOLD 

1238. 

35 

I=M(K> 

1239. 


IF(I-K) 45,45,38 

1240. 

38 

op=N*n-n 

1241, 


DO 40 0 = 1 , N 

1242. 


0K=NK+0 

1243. 


0I=0P+0 

1244. 


H0LD=-A< OK) 

1245. 


At OK )=A< 01 ) 

1246. 

40 

At 0 I ) =ROLD 

1247. 

45 

I F < CDABS( BIGA ) > 48,46,48 

1248. 

46 

D=0.0D0 

1249. 


RETURN 

1250. 

48 

DO 5S 1=1, N 

1251 . 


IFt I"K > 50,55,50 

1252. 

50 

IK=NK+I 

1253. 


At IK )=At IK )/( “BIGA > 

1254. 

55 

CONTINUE 

1255. 


DO 65 1=1 f N 

1 256. 


IK=NK+I 

1257. 


H0LD=A( IK) 

1258. 


I0=I-N 

1259. 


DO 65 0=1, N 

1260, 


10= IO+N 

1261. 


IF(I-K) 60,65,60 

1262. 

60 

IFtO-K) 62,85,62 

1263. 

62 

K0=I0-I+K 

1264; 


A( 10 ) = HOLD*A( KO >+A< 10) 

1265. 

65 

CONTINUE 

1266. 


K0=K-N 

1267. 


DO 75 0=1, N 

1268. 


K0=K0+N 

1269. 


IF(O-K) 70,75,70 

1270, 

70 

At KO )=At KO l/BIGA 

1271, 

75 

CONTINUE 

1272, 


D=D*BIGA 

1273. 


At KK ) = 1 *000/ BIGA 

1274, 

80 

CONTINUE 

1275. 


K=N 

1276. 
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100K=<K-1) 1277. 

IFt K ) 150 , 150, 105 1278. 

105 I =L ( K ) 1279. 

IF( I-K) 120,120,108 1280. 

108 JQ=N*t K— 1 } 1281. 

OR=N*( 1-1 ' 1282. 

DO 110 JO , N 1283. 

JK=0Q+0 1284. 

HOLD=Af OK ) 1285. 

0I=JR+0 1286. 

A(JK)=-A<01) 1287. 

110 At 0 I > =HOLD 1288. 

120 0=M( K ) 1289. 

IFtJ-K) 100,100,125 1290. 

125 KI=K-N 1291. 

DO 130 1 = 1, N 1292. 

KI=KI+N 1293. 

HOLD=A< KI ) 1294. 

JI=KI-l<+0 1295. 

ACKI)=-A<0n 1296. 

130 At JI ) =HOLD 1297. 

GO TO 100 1298. 

150 RETURN 1299. 

-END 1300. 

1301 . 

1302. 

1303. 

SUBROUTINE PR1NT1 ( ZEMB .ZEMBDC ,ZER,ZEI ,ZEMBSB , PRMT, Y , DERY , VLO, 1304. 

1VDBIAS,IDBIAS,RSSB,RSL0, NHARM, NHP 1 , NHD2 ) 1305. 

C 1306. 

C — PRINT1 WRITES THE VALUES OF THE INPUT VARIABLES AND THE INITIAL 1307. 
C— - CONDITIONS FOR THE NONLINEAR ANALYSIS SECTION OF THE PROGRAM. 1308. 

C- — 1309. 

C— -THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1310. 

C — FOR COMMON/ADJLO/: 1311. 

REAL*8 LOVLO, UPVLO, VLOINC , IDCACC 1312. 

INTEGER LOFLAG ,UPFLAG 1313. 

C FOR COMMON/CONST/: 1314. 

REALMS QEL , BOLTZ , P I ,TEQ 1315. 

C FOR COMMON/DIODE/: 1316. 

REALMS ALP , PHI , GAM, C0, IS,RS,RSKIN,FP,IF, IGJ , ICO , GJ ,CJ 1317. 

C FOR COMMON/LOOPS/: 1318. 

INTEGER NH,NLO,NVLO,NPTS,NPRI NT, NITER 1319. 

C FOR COMMON/RKG/: 1320. 

REAL*8 ACC , VDINIT 1321. 

INTEGER NDIM 1322. 

C FOR COMMON/TLINE/: 1323. 

REAL*8 Z0.ZQACC 1324. 

INTEGER ZQFLAG 1325. 

C FOR VARIABLES NOT IN ANY COMMON BLOCKS: 1326. 

C0MPLEX*16 ZEMBtNHARM) .ZEMBSBtNHPl ) 1327. 

REAL*B VDBIAS, IDB IAS , VLO ,ZQACC ,Z0,ZEMBDC , PRMTt 5 ) ,Yf 1 ) , DERYf 1) 1328. 

REAL *8 ZERt NHARM) ,ZEI ( NHARM ) ,RSSB( NHP1 ), RSLOt NHARM > 1329. 

INTEGER NHARM, NHP1 ,NHD2, I ,K,0 1330. 

C THE COMMON BLOCKS USED ARE: 1331. 

COMMON/AD JVLO/LOVLO, UP VLO, LOFLAG. UPFLAG, VLOINC, IDCACC 133P. 

COMMON/CONST/QEL ,BOLTZ,P I ,TEQ 1333, 

COMMON/D IODE/ALP, PHI , GAM , C0, IS ,RS , RSKIN , FP , IF , IGO , ICO ,G0 ,C0 1334. 
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COMMON /LOOPS/N H , NLO, NVLO, NPTS , NPRINT, NITER 1335. 

COMMQN/RKG/ACC ,VDINIT,NDIM 1336. 

COMMON/TLINE/Z0,ZQACC,2QFLAG 1337. 

C-*— -PRINT THE* TITLE 1338. 

WRITE! 6,50) FP 1333. 

SJ2f FORMAT! 1H1 , IX , ' ANALYSIS OF A V-SPF6.2.' GN2 MICROWAVE MIXER*/ 1340. 
*1X , 50! ' - ' ) ) 1341. 

C— WRITE THE VALUES OF THE RELEVENT VARIABLES. 1342. 

WRITE -(6,75) 1343. 

75 FORMAT! /IX, 'INPUT DATA* ) 1344. 

WRITE! 6, 100) 1345. 

100 FORMAT! //IX, ’ DIODE PARAMETERS: * ,T25 , ‘ALP ' ,T4l PHI ' ,T56 GAM ‘ , 1346. 

1T71 , ' RS ‘ ,T85 , 1 C0 1 ,T100, ' IS ' , TI13 , * RSKIN ' ) 1347. 

WRITE! 6,110} ALP, PHI, GAM, RS,C0, IS, RSKIN 1348. 

110 FORMAT! 19X,4(F10.3,5X),2!1PE10.3,5X),1PE10.3) 1349. 

WRITE! 6 , 1 20) FP , IF ,TEQ 1350. 

120 FORMAT! /IX. * OPERATING FREQUENCIES AND TEMPERATURE: ' .TS0, ' FP ' , 1351. 

1T65, 1 IF ' , T79 , *TEQ'/T45,2( 1PE 10. 3 , 5X ) .0PF10.3 ) 13EL. 

WRITE! 6, 130) VDBIAS , IDBIAS 1353. 

130 FORMAT! /IX, * BIAS SETTI NGS : 1 ,T24 VDB IAS * ,T30 IDBIAS * /T20, 1354. 

1F10. 3 ,5X , F10.6 ) 1355. 

WRITE (6,140) VLO, VLOINC , I DC ACC 1356. 

140 FORMAT! /IX, 1 ADJVLO VARIABLES :* ,T25 ,* VLO * ,T39 ,* VLOINC’ , T54, 1357. 

1 'IDCACC'/T2I,3!F10.6,5X)) 1358. 

WRITE! 6, 150) PRMT! 1 > , PRMT! 2 ) , PRMT! 3 ) ,PRMT( 4 > , Y(1 ) ,DERV( 1 > , NDIM 1359. 
150 FORMAT! /IX, ' DRKGS VARIABLES : ' ,T2 1 ,' PRMT! 1 ) ' ,T35 , 1 PRMT! 2 >' ,T50, 1360. 

1 ' PRMT! 3 ) ' ,T65, 'PRMT! 4 ) 1 ,T80, ‘ Y( 1 >' ,T95, ' DERY! 1 >' ,TU0, ‘NDIM’/ 1361 . 

2T20 , ' ! LOW L 1M ) 1 , T35 , ' ( UP L 1M ) ' ,T50, ' ! INCR ) ' ,T66 , ' ! ACC > ' , T80, 1362. 

3’ (VD) ’ ,T95, '(DV/DT) ’ ,T109, * ( NEQS > ’ / 1363. 

4T22 , F 10. 8 , IX , 2! F10.8 , SX > , 1 PE10. 3 ,2X , 2! 0PF10.3 , 6X > , 4X , 12 > 1364. 

WRITE! 6,160) NITER, NLO, NVLO, NPTS, NHARM, NPRINT 1365. 

160 FORMAT! /IX, ' LOOP L IMITS : ‘ , T2 1 NITER T31 NLO 1 ,T40, ' NVLO ' 1366, 

1T51, * NPTS * ,T62 , ' NHARM 1 ,T72, ' NPRINT ' /T21 - 14 , 6X , 2! 12 , SX ) , IX , 1367. 

212 ,09X , 12 ,8X, 13 ) 1368. 

WRITE! 6 , 1 70) 20,ZQACC 1369. 

170 FORKAT! /IX, ’ CONVERGENCE PARAMETERS: ' ,T40, '20' ,T57, 'ZQACC'/ 1370. 

1T34, F 10.2, 10X,1PE10.3) 1371. 

WRITE ('6, 2 80) ZEMBDC,ZEMBDC,NHD2,ZEMBSB! 1 ) L372. 

180 FORMAT! ///IX, 1 EMBEDDING IMPEDANCES : 1 , T48 HARMONICS OF THE L0‘ 1373. 

1TI05, ’HARMONIC SIDEBANDS ‘ /T25 , ' HARM#' ,T37 , ’ZER’ ,T50, 1374. 

*'2EI',T71, 'ZEMB 1 ,T92, 1 SIDEBAND#' ,TU2, ' ZEMBSB ' /T26 , 'DC ,133, 1375. 

*1PE10~3 ,T61 , 1PE10. 3 , T95, 12 , T1J&3 ,1PE10.3,T116,1PE10.3) 1376. 

DO 10 1=1, NHARM 1377, 

K=NHD2-I 1378. 

0 - 1*1 1379. 

10 WRITE! 6, 190) 1 , 2ER! I ) ,ZE I ( I) ,2EMB< I) , K , ZEMBSB! J ) 1380. 

190 FORMAT! 1X,T26, 12,T33,2! 1PE 10. 3 , 3X ) , T6 1 , 2! 1 PE 10. 3 , 3X ) , T95, 12 , 1381. 

*T103 , 1PE 10. 3 , T116 , 1 PE 10. 3 ) 1382. 

WRITE! 6,200} RS,NHD2,RSSB< 1) 1383. 

200 FORMAT! ///IX, 1 DIODE SERIES RESISTANCES ,T32, ' HARMONICS OF THE 1384. 
1,' LO’ ,T64, 'HARMONIC SI DEBANDS ' /T34 HARM# ' ,T44 RS! F )' ,T62 , 1385. 

2 'SIDEBAND#' ,T77, 'RS(F ) VT35, 'DC* ,T41 ,F3 . 4 ,T65 , 12 , T74 , FB . 4 } 1386. 

DO 20 1=1, NHARM 1387. 

K=NHD2-I 1388. 

0 = 1 + 1 . . 1389. 

20 WRITE! 6,210) I ,RSLO! I) , K,RSSB(0 ) 1390. 

210 FORMAT! 1~! , T36 , 12,T4 l , F8 . 4, T65, 12,T74,F8,4> 1391. 

C — -WRITE THE HEADING FOR THE. NEXT SECTION OF PRINTOUT 1392. 
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WRITE! 6 ,220) 1393. 

220 FORMAT! I HI, 'RESULTS OF THE VLO ADJUSTMENTS* // 3 1394. 

RETURN 139S. 

END 1396. 

1397. 

1398. 

1399. 

SUBROUTINE PRINT2! RHO, VI , VR, DVR, VDCOS , VDSIN , IDCOS, IDSIN , ZQMAG , 1400. 

1 ZQPHA, VRDC, DVRDC, VI DC, RHODC,AV,BV, ITER ,ZQF LAG ,0 VLO , NH , NHP 1 > 1401 . 

1402. 

PRINT2 WRITES THE RESULTS OF EACH REFLECTION CYCLE OF THE LOOP 1403. 

ITER IN SUBROUTINE LGSIG . 1404. 

1405. 

—THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1406. 

COMP LEX* 16 RHOt NH 1 ,VR{ NH ) . VI ( NH ) ,DVR! NH ) Hfi 7. 

REALMS VDCOS! NHP 1 ), VDSIN! NHP 1 ), IDCOS! NHP 11 , IDSIN! NHP 1 ) , ZQMAG! NH ) 1408. 

REAL*8 ZQPHA! NH ) , VIDC , VRDC , DVRDC , RHODC , AV! NH ) , BV( NH ) 1409. 

INTEGER ITER , J VLO, NH , NHP 1 1410. 

—WRITE THE RESULTS OF THE REFLECTION CYCLE 1411. 

WRITE! 6 , 100) ITER, JVLO 1412. 

100 FORMAT! ////IX, 'NONLINEAR ANALYSIS RESULTS: REFLECTION CYCLE #’ 1413. 

1,14,' IN VLO ADJUSTMENT LOOP NUMBER', 13/) 1414. 

WRITE! 6 , 110) 1415. 

110 FORMAT! /2X, ‘AVtD.BV! I)' ) 1416. 

WR I TE (6,120) 1 1 ,AVC I ) , BV( I) , I ■= 1 , NH ) 1417. 

120 FORMAT! 1H*,6!BX,4(I7,1PE12.3,1PE12.3)/1X>) 1418. 

WR I TE ! 6 , 1 30 ) 1419. 

130 FORMAT! /2X, 'VR' ) 1420. 

WRITE! 6 , 1 20) ( I , VR! I ) , 1 = 1 , NH ) 1421. 

WRITE ! 6 , 1 40) 1422. 

140 FORMAT! /2X, 'DVR' ) 1423. 

WRITE! 6 , 120) 1 1 ,DVR( 1 ) , 1 = 1 , NH ) 1424. 

WRITE! 6 , 1 50) 1425. 

150 FORMAT! /2X, 'VDCOS, VDSIN ' ) 1426. 

WRITE! 6 , 1 20 ) ! I , VDCOS! 1 + 1 ), VDSIN! I + 1 >,I = 1 ,NH) 1427. 

WRITE! 6 , 1 60 ) 142B. 

160 FORMAT! /2X,‘ IDCOS, IDS I N* ) 1429. 

WRITE! 6 , 120) U,IDC0SU + 1>,IDSINU + 1),I<*1,NH> 1430. 

WRITE! 6 , 1 70) 1431. 

170 FORMAT! /2X, ’ZQMAG, ZQPHA ' ) 1432. 

WRITE (6, 180) { I , ZQMAG! I ) , ZQPHA! 1 ) , 1=1 , NH ) 1433. 

180 FORMAT! IH+,6(8X,4(I7,1PE12.3,0PF7.0,5X)/1X)) 1434. 

WRITE ( 5 , 1 90 ) VDCOS! I), IDCOS! 1), VRDC, DVRDC, ZQFLAG 1435. 

190 FORMAT! //2X, ’DC TERMS: VDCOS= * , 1 PE 10. 3 , T35 , ’ IDCOS= ’ , 1PE 10. 3, 1436. 

1T54 , VRDC = ' , 1PE 10. 3 , T7S , ' DVRDC=' , 1PE10.3///2X, ' ZQF LAG= ' , 12) 1437. 

RETURN 1438. 

END 1439. 

1440. 

1441. 

1442. 

SUBROUTINE PRINT3C FG , FC ,GOMAG , GOPHA.CJMAG , CJPHA , NH , NHP I ) 1443. 

1444, 

PRINT3 WRITES THE FOURIER COEFFICIENTS OF THE DIODE CONDUCTANCE 1445. 

—AND CAPACITANCE WHICH ARE USED IN THE SMALL-SIGNAL ANALYSIS. 1446. 

1447. 

-—THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1448. 

COMPLEX* 16 FG! NHP 1 ) ,FC{ NHP 1 ) 1449. 

REAL*8 GJMAG! NHP 1 ) ,GJPHA( NHP1 ) , COMAG! NHP 1 ) , CJPHA! NHP 1 ) 1450. 
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INTEGER NHPl.NH 1451. 

C— TRANSFORM THE FOURIER COEFFICIENTS TO MAGNITUDES (HALF AMPLITUDES) 1452. 
C— AND PHASES f IN DEGREES) 1453. 

DO 10 I-l.NHPl 1454. 

GOMAG! I >«CDABS! FG! I ) ) 1455. 

C0MAG1 I )=>CDABS( FC< I ) ) 1456. 

GOPHA! I )»DATAN2< DINAG! FG! I ) ) ,DREAL(FG( I ) ) >*57.2957795100 1457. 

10 COPHA! I )=DATAN2( D IMAG! FC( I ) ),DREAL!FC( I ) > >*57. 2957795 1D0 1458. 

WRITE! 6,50) 1459. 

50 FORMAT! 1H1 , 'RESULTS OF THE SMALL-SIGNAL ANALYSIS'/) 1460. 

WRITE! 6 , 100 > 1461. 

100 FORMAT! /IK, 1 FOURIER COEFFICIENTS OF THE DIODE*, 1462. 

1* CONDUCTANCE AND CAPACITANCE WAVEFORMS*/) 1463. 

WRITE! 6 , 1 10) 1464. 

110 FORMAT! /2X, ‘GOMAG .GJPHA ‘ > 1465. 

WRITE t 6 , 120) (I, GOMAG! 1+1), GOPHA! 1 + 1 >,I-1,NH> 1466, 

120 FORMAT(lH-i-,6<8X,4U7,lPE12.3,0PF7.0,5X>/lX>> 1467. 

WRITE ( 6 , 1 30 ) 1468. 

130 FORMAT! /2X, 'COMAG, COPHA' ) 1469 . 

WRITE! 6 , 120 > ( I, COMAG! 1 + 1 >, COPHA! 1 + 1 ), 1 = 1, NH> 1470. 

WRITE! 6, 140) GOMAG! 1 ), COMAG! 1 ) 1471. 

140 FORMAT! //2X, 'DC TERMS: GOMAG = ',1PE 10. 3 , 4X , * COMAG * * 1472. 

1 , 1 PE 10.3/ > 1473. 

RETURN 1474. 

END 1475. 

1476. 

1477. 

1478. 


SUBROUTINE PRINT4! XLMAT.ZI N ,ZIFOUT,TM, THERM , SHOT, NHD2 , NHP l , NHD2P 1 ) 1 479 . 


C 1480. 

C— SUBROUTINE PRINT4 WRITES THE RESULTS OF THE SMALL-SIGNAL 1481. 

C— AND NOISE ANALYSIS. THESE INCLUDE THE CONVERSION LOSS MATRIX, 1482. 

C— THE INPUT IMPEDANCE AT EACH SIDEBAND, THE IF OUTPUT IMPEDANCE, 1483. 

C AND THE EQUIVALENT INPUT NOISE TEMPERATURE OF THE MIXER WITH ITS 1484. 

C— THERMAL AND SHOT NOISE COMPONENTS. 1485 . 

C 1486. 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1487. 

COMPLEX*l 6 ZIN! NHP 1 ) .ZIFOUT 1488. 

REAL*8 XLMAT ( NHP 1 , NHP 1 > , TM .THERM , SHOT , REZIN! 9 ) , IMZIN! 9 ) 1489. 

INTEGER NHP1 , 1 ,0.,K,NHD2P1 , NHD2P2 ,'NHD2 1490. 

C— DEFINE THE ARRAY ELEMENT WHICH CORRESPONDS TO THE LOWER SIDEBAND 1491. 
NHD2P2 a NHD2Pl+l 1492. 

WRI TE < 6 , 100 ) 1493. 

100 FORMAT! //44X. 'CONVERSION LOSS MATRIX ( DB ) ' / ) 1494. 

WR ITE t 6 , 200 ) 1495. 

200 FORMAT! IX. T25, ' 4 ' ,T35 , * 3 ' ,T45 , ' 2 ’ ,T55 , * 1 * ,T65, *0* ,T74, '-1 ' ,T84, 1496 . 

I'-2' , T94 , ' -3 ' ,T104 , *-4 * // ) 1497. 

DO ! 0 I “ 1 , NHP 1 1498. 

K=NHD2P 1 - i 1499. 

10 WRITE! 6,300) K, < XLMAT! I,0),0=1,NHPI) 1500. 

300 FORMAT! 9X,I4,4X,9F10.2 ) 1501. 

WRITE < 6 , 350 ) XLMAT! NHD2P1 ,NHD2) , XLMAT! NHD2P 1 , NHD2P2 ) 1502. 

350 FORMAT! //2X, 'UPPER SIDEBAND CONVERSION LOSS: L!0,1) =>*,F7.2, 1503. 

*' DB' //2X. 'LOWER SIDEBAND CONVERSION LOSS: L!0,-1) »*,F7.2, 1504. 

*' DB ' / ) 1505. 

WRITE ! 6,390) 1506. 

390 FORMAT! 2X, 'INPUT IMPEDANCES T25 ,' 4 * ,T35 ,' 3 *, T45 ,* 2 ' ,T55 , 1 1 ' , 1507. 

*T65, '0' ,T74, ' -1 ' ,T84 , ' -2 ' ,T94, ' -3' ,T104, '“4* ) 1508. 


DO 20 1 = 1 , NHP1 1BJ0T9. 

REZIN! I )=DREAL(ZIN( I ) ) 1510. 

20 IMZIN! I )=DIMAG(ZIN( I ) ) 1511. 

WRITE! G ,400} CREZIN(I),I=1,NHP1) 1512. 

400 FORMAT !/3X, ' REAL(ZIN): ' ,3X,9F10.2) 1513. 

WRITE! 6 ,500) {IMZIN! I),I=1,NHP1) 1514. 

500 FORMAT !/3X, ’ IMAG! ZIN ) : ' , 3X , 9F10, 2 ) 1515. 

WRITE (6,550) ZIFOUT 1516. 

550 FORMAT! //2X, ’ IF OUTPUT IMPEDANCE: ZIFOUT = '.F8.2,' + 0’ 1517. 

*F8 . 2 ) 1518. 

WRITE! 6,600) TM, THERM, SHOT 1519. 

600 FORMAT! //2X, 'EQUIVALENT INPUT NOISE TEMPERATURES: ' ,T52, ’TM* 1520. 

1TG7, 'THERM' ,T82, 'SHOT' ,T98/T46,F10. I ,T62,F10. 1 ,T77,F10. 1 ) 1521. 

RETURN 1522. 

END 1523. 

1524. 

1525. 

1526. 

SUBROUTINE PLOTt I GO DAT, CO DATA, VDDATA, IDDATA , NPTS , ITER.C0) 1527. 

C 1528. 

C SUBROUTINE PLOT GRAPHS THE CURRENT THROUGH THE DIODE CONDUCTANCE 1529. 

C ( IGO ) ,THE DIODE CAPACITANCE (CO), THE TOTAL CURRENT THROUGH THE 1530. 

C DIODE (IGO+ICO), AND THE VOLTAGE ACROSS THE INTRINSIC DIODE TER - 1531.. 

C M INALS ! V! 1 ) > (WHICH DOES NOT INCLUDE THE DIODE SERIES RESISTANCE) 1532. 

C AS FUNCTIONS OF TIME, OVER ONE LOCAL OSCILLATOR CYCLE. 1533. 

C— - 1534. 

C THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 1535. 

REAL*8 IGODAT! NPTS ) , CODATA! NPTS ), VDDATA! NPTS ) , IDDATA! NPTS ) 1536 . 

REAL*8 MAX I GO ,MAXCO ,MAXVD , MAX ID ,MINCO ,C0 1537. 

INTEGER ITER, OPT, YPT, NPTS, IGJPOS .COPOS ,C0POS , VDPOS , I DPOS,C, ZERO 1538. 

INTEGER BLANK, DOT,STAR,YGPOS(50),YCPOS(50),YIDPOS! 50 ) , YVDPOS! 50) 1539. 

C DEFINE THE NUMERICS USED IN THE GRAPHS 1540. 

DATA BLANK, DOT, STAR, C, ZERO/' C '0' / 1541. 

C DETERMINE THE GRAPH SCALES 1542. 

MAXIGO=DABS( IGODAT! 1 ) ) 1543. 

MAXCO=DABS( CODATA! 1 ) ) 1544. 

MI NCO=CODATA( 1 ) 1545. 

MAXVD=DABS( VDDATA! 1 ) ) 1546. 

MAXID=PABS( IDDATA! 1 ) > 1547. 

DO 10 0PT=2, NPTS 1548. 

IF! MAX I GO. LT. DABS! IGODAT! OPT >) ) MAXIGO=DABS( IGODAT(OPT)) 1549. 

IF (MAXCO . LT. DABS! CO DATA! OPT ) ) ) MAXCO = DABS( CODATA! OPT) ) 1550. 

IFfMlNCO.GT. CO DATA! OPT ) ) MINCO=CJDATA( OPT ) 1551. 

I F ( MAXVD . LT. DABS! VDDATA! OPT ) ) ) MAXVD=DABS( VDDATA! OPT ) ) 1552. 

IF ( MAX I D , LT. DABS! IDDATA! OPT ) ) ) MAXI D=DABS( I DDATA! OPT ) ) 1553. 

10 CONTINUE 1554. 

C THE GRAPH HEADINGS 1555. 

WRITE! 6 , 100) ITER 1556. 

100 FORMAT! INI , IX, 'GRAPHS FOR REFLECTION CYCLE NUMBER ',14/) 1557. 

WRITE! 6 , 1 10) 1558. 

110 FORMAT! /3X, ' IGO(MA)’ ,5X, ' DIODE CONDUCTANCE CURRENT VS TIME FOR', 1559. 
1' ONE LO CYCLE', T67, ' CO! PF ) ' , 5X , ’ D IODE CAPACITANCE VS', 1560. 

2' TIME FOR ONE LO CYCLE’/) 1561. 

C THE LOOP FOR THE POINTS TO PLOTTED VERTICALLY DOWN THE PAGE 1562. 

C— PREVENT A DIVISION BY 0 IF THE CAPACITANCE DOES NOT VARY 1563. 

IF! MAXCO. EQ.MINCJ) MAXCO=MAXCO+l .0D0 1564. 

C LET C0 BE THE Y AXIS IF IT IS NOT IN THE RANGE MINCJ TO MAXCO 1565. 

I F ( C0. LT . MINCJ ) MI NCJ=C0 1566. 


IF<C0.GT.MAXCJ> MAX CO =C0 1567. 

C0POS=01+DINT( 50.0D0/! MAXCO-MINCO )*( C0-MINCJ )+DSIGN(0.5D0,C0) ) 1568. 

DO 2 OPT= 1 , NPTS 1569. 

IGOPOS=01+DINT( 50.0D0/MAXIGO*IGODAT! OPT)+DSIGN(0.5D0, I GO DAK OPT) > >1570. 
COPOS=01 +DINT! 50.000/ ! MAXCO-MINCO )*< C0DATA( OPT J-MINCO )+DSIGN< 1571. 

10.5D0,CODATA(OPT) > ) 1572. 

C SET THE GRAPH LIMITS 1573. 

IF( IGOPOS.LT. I) I GOPOS= 1 1574. 

IF! IGOPOS.GT.50J IGOPOS=50 1575. 

IF < COPOS . LT. 1 > COPOS=1 1576. 

I F( COPOS . GT. 49 ) COPOS=49 1577. 

C— CLEAR THE HORIZONTAL LINE 1578. 

DO 1 YPT=1 ,50 1579. 

YCPOS! YPT >=BLANK 1580. 

1 YGPOS! YPT)=BLANK 1581. 

C SET THE GRAPH'S Y AXIS 1582. 

YGPOS! 1 )=DOT 1583. 

YCPOS{ C0POS ) = DOT 1584. 

C— THE PLOTTED POINTS ARE REPRESENTED AS ASTERIKS 1585. 

YGPOS! IGOPOS)=STAR 1586. 

YCPOS!COPOS)=STAR 1587. 

C— -WRITE 'C0' ON THE Y AXIS OF THE CAPACITANCE GRAPH 1588. 

I F { C0POS . EQ. 50) GOTO 6 1589. 

IF i OPT . EQ . 1 ) YCPOS< C0POS )=C 1590. 

IFCOPT.EQ. 1 > YCPOS! C0POS+1 )=ZERO 1591. 

GOTO 7 1592. 

6 IFIOPT.EQ.l) YCPOS(C0POS-n=C 1593. 

IF(OPT.EQ.l) YCPOS( C0POS )=ZERO 1594- 

7 CONTINUE 1595. 

C— PRINT THIS LINE OF THE GRAPHS 1596. 

WRITE! 6 , 120) I GO DAT! OPT-) ,( YGPOS! YPT) , YPT=1 , 50 ), CODATA! OPT > , 1597. 

I! YCPOS! YPT) , YPT= I ,50) 1598. 

120 FORMAT! 3PF9.3.2X.50A1 ,3X, 12PF9.4.2X.50A1 ) 1599. 

2 CONTINUE 1600. 

WRITE! 6, 100) ITER 1601. 

WRITE! 6,130) 1602. 

130 FORMAT! //3X, 1 ID(MA) ' ,5X, ’TOTAL DIODE CURRENT VS TIME FOR ONE L0\ 1603. 
1* CYCLE ' ,T67 , 1 VD< VOLTS )', 8X , ' DIODE VOLTAGE VS TIME FOR', 1604. 

2' ONE LO CYCLE'/) 1605. 

C THE DO LOOP FOR THE POINTS TO BE PLOTTED VERTICALLY DOWN THE PAGE 1606. 

DO 4 0PT=1 ,NPTS 1607. 

IDP0S = 25+DINT! 25 .0D0/MAX I D*I DD'aTA! OPT) + DSIGN(0.5D0, IDDATA! OPT) ) ) 1608. 

VDPOS=25+DINT! 25 .0D0/MAXVD*VDDATA! OPT >+DS IGN! 0. 5D0, VDDATAC OPT) > ) 1609. 

C SET THE GRAPH LIMITS 1610. 

IF! IDPOS.LT. 1 ) IDP0S=1 1611. 

IF! IDPOS.GT.50) IDPOS=50 1612. 

IF! VDPOS . LT. 1 ) VDP0S=1 1613. 

IF! VDPOS.GT .50) VDPOS=O0 1614. 

C— CLEAR THE HORIZONTAL LINE 1615. 

DO 5 YPT=1 ,50 1616. 

YIDPOS! YPT)=BLANK 1617. 

Y VDPOS! YPT >=B LANK 1618. 

5 CONTINUE 1619. 

C SET THE V AXIS 1620. 

YIDPOS!'25 )=DOT 1621. 

YVDPOS! 25 )=DOT 1622. 

C— THE PLOTTED POINTS ARE REPRESENTED AS ASTERIKS 1623. 

YIDPOS! IDPOS)=STAR 16Z4. 
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YVOPOS{ VDPOS )«STAR 1625. 

C- — PRINT THIS LINE OF THE GRAPHS 1626. 

WRITE! 6, 14J0T) IDDATAtOPT), ( V IDPOS! YPT > , YPT=1 , 50) , VDDATA! JPT), 1627. 

H YVOPOSt YPT) , YPT-1 ,50) 1528. 

14J0T FORMAT! 3PF9 , 3 , 2K » SXTA1 ,3X,0PF9.3,2X, 50A1 ) 1629. 

4 CONTINUE 1630, 

3 RETURN 1631. 

END 1632. 

1633. 

1634. 

1635. 

BLOCK DATA 1636. 

1637, 

—FOR COMMON/ADOLO/: 1638. 

R£AL*8 LOVLO , UPVLO , VLOINC , IDCACC 1639. 

INTEGER LOFLAG.UPFLAG 1640. 

C — FOR COMMON/CONST/: 1641. 

REAL*8 QEL , BOLTZ,PI , TEQ 1642. 

C FOR COMMON/DIODE/: 1643. 

REAL *8 ALP,PHI,GAM,C0,IS,RS,RSKIN,FP,IF,IGO,ICJ,GJ,CO 1644. 

C— FOR COMMON/IMPED/: 1645. 

C0MPLEX*16 ZEMBSB( 9 ) 1646. 

REALMS ZER! 8 ) ,ZE I ( 8 } .ZEMBDC ,RSL0{8),RSSB(9) 1647. 

C FOR COMMON/LOOPS/: 1648. 

INTEGER NH , NLO, NVLO.NPTS, NPRINT, NITER 1649. 

C FOR COMMON/RKG/: 1650. 

REAL*8 ACC.VDINIT 1651. 

INTEGER NDIM 1652. 

C— FOR COMMON/TLINE/: 1653. 

REAL*8 Z0.ZQACC 1654. 

INTEGER ZQFLAG 1655. 

C — -FOR COMMON/VOLTS/: 1656. 

REAL*8 AV! 8 ) , BV! 8 ), VI DC , VLO, VDB IAS , IDB IAS 1657. 

C THE COMMON BLOCKS USED AREs 1658. 

COMMON/ADOVLO/LOVLO, UPVLO, LOFLAG,UPFLAG, VLOINC, IDCACC 1659. 

COMMON/CONST/QEL .BOLTZ.PI ,TEQ 1660. 

COMMON/DIODE/ALP, PHI, GAM, C0, IS, RS,RSKIN,FP,rF,IGJ,ICJ,GJ, CO 1661 . 

COMMON/ IMPED/ZER.ZEl , ZEMBDC, ZEMBSB ,RSLO , RSSB 1662. 

COMMON /L OOPS /N H, NLO, N VLO, NPTS, NPR I NT, NITER 1663. 

COMMQN/RKG/ACC , VDINIT, NDIM 1664. 

COMMON/TL INE/Z0,ZQACC , ZQFLAG 1665 . 

COMMON/ VOL TS/AV , BV , VIDC , VLO , VDB I AS , IDBIAS 1666. 

C VARIABLES ARE INITIALIZED AS FOLLOWS: 1667. 

C— COMMON/ADOVLO/ VARIABLES: 1668. 

DATA VLOINC, IDCACC/0. 500, 0.0200/ 1669. 

C— COMMON/CONST/VARIABLES: 1670. 

DATA QEL , BOLTZ, P 1/ 1 . 602 192D-I9, 1 . 3806 2D -2 3 , 3 . 14159265 3589 79D0/ 1671 . 

DATA TEQ/296.0D0/ 1672. 

C COMMON/DIODE/VARIABLES: 1673. 

DATA ALP, PHI, GAM/40. 0D0,1 . 1D0.0.5D0/ 1674. 

DATA C0, IS, RS/2. 00-13,5. 00-9,5.000/ 167F. 

DATA FP, IF/15. 009,15.008/ 1676. 

DATA RSKIN/4.7434D-6/ 1677. 

C COMMON/IMPED/VARIABLES: 1678. 

DATA ZER( 1 ),ZER<2),ZER{3>/48. 000,64. 4800,29. 02D0/ 1679. 

DATA ZER! 4 ) ,ZER( 5) ,ZER{ 6 )/59. 18D0.43. 5D0, 68 . 3100/ 1680. 

DATA ZER! 7) ,ZER!8>/61 .4500,82 .4900/ 1681. 

DATA ZEK 1 ) ,ZE I C 2 ) ,Z£ I ( 3 ) / 1 8 . 74P0, 99 .1500, 166. 6 D0/ 1682. 
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DATA ZEI(4),ZEIt5),ZEI(6)/232. 4D0,299 .2D0, 367 . 8D0/ 1683- 

DATA ZEI<7),ZEI(8)/438.5D0,511. 4D0/ 1684. 

DATA ZEMBSBt 1 ) .ZEMBSBt 2 )/( 59 . 18D0, 232 . 4 DJ0T ) , ( 29 .02D0, 1 66 . 6D0 >/ 1685. 

DATA ZEMBSBt 3 ) , ZEMBSBt 4 >/( 64 .48D0, 99 . 1 5D0 ) , ( 48 .0D0, 1 8 . 74DJET ) / 1686, 

DATA ZEMBSBt 5 ) ,ZEMBSB( 6 )/ 1 BUS .80S ,0 .000) , ( 48 . JSTDJ0T , -18.774D0)/ 1687. 

DATA ZEMBSBt 7), ZEMBSBt 8 )/(64.48D0, -99. 15D0>, ( 29 .02D0, - 1 66 . 6D0 >/ 1688. 

DATA ZEMBSB(9)/(59. 1800,-232.400)/ 1689. 

DATA ZEMBDC/50.0D0/ 1690. 

C COMMON/LOOPS/VARIABLES: 1691. 

DATA NH,NLO,NPTS, NVLO, N ITER , NPRI NT/8, 1,51 ,50, 500 , 100/ 1692. 

C COMMOM/RKG/VARIABLES : 1693. 

DATA VDINIT, ACC, NDIM/0, 000, 1.00-6,1/ 1694. 

C COMMON/TL IKE /VARIABLES: 1695. 

DATA Z0 , ZQAC C / 200 . 0D0 , 0 , 0 1 D0/ 1696. 

C— COMMON/VOLTS/VARIABLES: 1697. 

DATA VDBIAS, IDBIAS/0.0D0,0.002D0/ 1698. 

DATA VLO/2.50D0/ 1699. 

END 1700. 
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ANALYSIS OF a 15.00 ChZ MICRO * *VE MIXER j 


INPUT DATA 


CIQDE PARAMETERS 

: ALP 

PHI 

GAM 

RS 

CO 

IS 

R5KIN 


40.000 

1. 100 

G. 500 

5.000 

2. 000 D— 13 

S. 0000-09 

4.7420-06 

CFERATING FREQUENCIES AND TEMPERATURE ! 

FP 

IF 

TEQ 






I.EOOD+IO 1 

• 5 GOD +09 

296.000 



EIAS settings: 

VCBIAS 

IDBIAS 







a .a 

O.C0200Q 






ACJVLC VARIABLES 

: vlc 

VLO INC 

I DCACC 






2 .500000 

a.saaoao 

0. C200 00 





CHKGS VARIAALES! 

PR MT { 1 ) 

PRMT (2 ) 

PRMT (3) 

PRMT (4 J 

YU 1 

DERY (1 ) 

NDXM 


{LOW LIH) 

(UP LIM) 

C I NCR ) 

{ ACC) 

CVD) 

CDV/DT) 

(NEQS ) 


0 .0 

0.12319971 

C.l 2219971 

1.0 COO— 06 

0.0 

1*000 

1 


LCOF LIMITS: NITER NLC NVLO NP75 NHARM SPRINT 

500 I SC 51 a 100 


CONVERGENCE PARAMETERS 5 20 20 ACC 

2C0.0Q 1 ■ 00 OD — 02 


EMBEDDING IMPEDANCES t 

HARM# ZER 

DC 2.00QD+0 2 

1 4*aOOD+Cl 

2 6.448D+01 

2 £ .9020+0 I 

4 5*9180+01 

5 4 .350D+O I 

6 6*8310+01 

7 6.14SD+01 

£ 6.149D+0I 


HARMONICS OF THE LO 

ZE I ZEMB 

2 . 0000*02 

1.274D+01 2.0000*92 0*0 

9*51 5D+0 I 7.0300*01 9.91SD+01 

1. £660+02 3. 5C3 0 +01 1 .6660+02 

2.224D+02 6.5340*01 2. 3240+02 

£.9920+02 4* S BO 0+01 2.9920+02 

2*6760+02 7. 4730+01 3.6760+02 

4*3850+02 6.799D+01 4.38SD+02 

5. 1140+02 8.813D+0I 5.U4D+02 


HARMONIC SIDEBANDS 
SIDEBAND# ZEMBSB 

4 S.91ED+01 2.324 0+02 

3 2.9C2D+01 1.G66D+02 

2 6.44ED+01 9.915D+01 

1 4.8C0D+0I 1. £740+01 

0 S.OCCD+Ol 0.0 

-1 4.8CCD+0 I -I.S77D+02 

-2 6 .4420*0 1 —9 * S15D*01 

-3 2.9C2D+01 -1.666D+02 

-4 5*9 IED+0 1 -2.3240+02 


CIODE SERIES RESISTANCES! HARMONICS OF Tl-E LO 

HARM# RSCF) * 

DC 5 . 0 CC 0 

1 5. 5CCS 

2 S.GEU 

3 £.QC££ 

4 6.1519 

5 6.25SC 

6 £.4£3C 

7 6.527C 

8 G.G<2£ 


HARMONIC 
SI DEBAND# 

4 

3 

2 

1 

0 

-1 

-2 

-3 

-4 


SIDEBANDS 
RS.(F) 
6.1763 
6*0229 
5.8419 
5.6093 
5.1 837 
5.5511 
5.8008 
5*9893 
6. 1473 


i 



POULTS OF THE VLC ADJUSTMENTS 



VALUES DF THE DC CURRENT AND LO VOLTAGE FOR CKLE 1 OF THE LCCP FOR ADJUSTING VLD TC GIVE 0*002000 AMPS ARE: 

iocast ii=n . ra2 ms vlo before adjustment: e.sqoqc 

VLO AFTER ADJUSTMENT: 2.00000 

VALUES OF THE DC CURRENT AND LO VOLTAGE FOR CKLS 3 OF THE LCCP FOR ADJUSTING VuO TO GIVE 0*002000 AHPS ARE: 
I OCCSt 1>=T .001 7 <!9 VLD BEFORE ADJUSTMENTS £.00000 

VLO AFTER ADJUSTMENT : 2.250C0 

VALUES OF THE DC CURRENT AND LO VCLTAGE FOR CYCLE 3 OF THE LCCP FOR ADJUSTING VLO TC GIVE G. 002000 AHPS ARES 
IDCOSt 1)=3 .001945 VLO BEFORE AC JUS 7MENT : 2.25000 

VLO AFTER ADJUSTMENT: £.37500 

VALUES QF THE DC CURRENT AND LO VOLTAGE FOR CTCLE 4- OF THE LCCP FOR ADJUSTING VLD TO GIVE 0*002000 AHPS ARE* 
IDCOSt 13=0 .002044 VLD BEFORE ADJUSTMENT: £.37500 

VLO AFTER ADJUSTMENT: S. 31250 

VALUES OF THE DC CURRENT AND LO VOLTAGE FOR CYC LE S QF THE LCCP FOR ADJUSTING VLD TO GIVE 0.002000 AHPS ARE: 
IDCOSt 1 ) =0 .001554 VLO BEFORE ADJUSTMENT: £.312£0 

VLD AFTER ADJUSTMENT: £.31250 


Ni C NL I NEAR ANALYSIS RESULTS: 

REFLECTION 

CYCLE * 

IS IN VLQ 

ADJUSTMENT 

LCOP NUMEER 5 




| 

A V( I J i £V { I) 

I 

1 *1 SfiD+CO 

0.0 

2 

-3.7D7D-0 2 

-5.629D-02 

3 

5 *90 00—0 3 

-1 *1420-02 

A 

3.095D-03 

-3. 9530-03 a 


5 

3+377 C— 0 3 

-2.250D-C3 

e 

1 .8040-02 

-l .155D-0 3 

7 

1 .2670—03 

-6.206D-04 

e 

9.29BO-04 

—3. SI ED— 04 |j 

\R 

I 

-0.3230-0 1 

-4 . 1GGP-0I 

2 

l « 19 CD— 0 1 

3.07CO-03 

3 

l .2660-02 

-9 .451 0-03 

4 

5.8730-03 

-2. 987D-03 * 


5 

3.6 1 2D— 0 1 

- : .4 A 8D— 03 

6 

2.3310-02 

-9.8630-04 

7 

1 *467D-D3 

-5.955D-04 

a 

1 .0700—03 

-3.ES4D-04 i 

LVR 

1* 

U60SD-C7 

— 2 * 7 5 7 D— C 7 

2 

2.S4 6D-0 7 

-8. 7040-07 

3 

5.S63C-06 

-7.002C-06 

4 

-4.904D-07 

-1.4160-06 


5 

1 **'930-06 

-4*6770-07 

6 

-l .095D-05 

-4*2920-0 7 

7 

-l .40 5 D— 05 

1.0310-05 

e 

-2.1B4D-05 

7.26SD-0 6 

VDCCS +VDSIN 

t 

2*2400—0 1 

4. 1C6D-C1 

2 

8+ 173D-02 

-5.936D-0 2 

3 

1 .6660-02 

-1.965C-03 

4 

8.S73D-03 

! 

—9. £660—04 t 



5.9880-02 

-9.0I7D-04 

6 

4.14CD-03 

-1 .7590-04 

7 

2.7490-03 

-2.036D-O5 

e 

2*0190—03 

—3. 2250-05 i 

t 

ICCCSf IDSIN 

1 

1 .D 44D— 0 2 

-2.0E3D-C3 

2 

— 7 • 67 3D— 0 4 

— 2 *66 ID— 04 

3 

-3 *379C— 05 

— 1 .0438—04 

4 

-1.2920-05 

1 

—3. 47QD-tOS i 


s 

-6.1 SOD— 0 6 

-l .890D-C5 

6 

-2.60 8D-06 

- J -075D-05 

7 

-9.241 C-n 7 

-6.096C-06 

a 

-6. 0820-07 

-3. €310-06 * 

I CHAG « 2GPH A 

1 

2.1 97C-01 

-72. 

2 

l.OOCD+OG 

1 SC* 

3 

9*9990—0 1 

180. 

4 

9.99BD-01 

IBS. ! 


5 

i .oooDvoa 

iec. 

6 

9.9S4D-G 1 

-I so. 

7 

1 .0050+00 

-179. 

8 

1 • 003D+00 

i 

cc terms: vdcos 

S-B.Q44C-03 

IDCQS= 

1*59*C- 

C3 VR0C=-2 

*03eo-oi 

DVRDC 

=-1.4990-07 




r 


2CFLAG=> 0 


REQUIRED LO PCUEfi: 1.474C-03 


‘ VJD 

H 


CffAPHS FOR REFLECTION CYCLE NUMEER 
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IGJ(PA) 


CIO CE CONDUCTANCE CURRENT VS TIFE FQF ONE LO CYCLE CJ(PF) 


Cl ODE CAPACITANCE VS TINE FCR ONE UP CYCLE 


S.147 


10.216 


10 .756 


10*310 


9. 72£ 


9.074 


8.371 


7.606 


6.773 


5.261 


4.54S 


4.004 


3.076 


2.200 


1 *441 


0.027 


0.401 


0.158 


0*049 

♦ 

0 .012 

* 

0 .002 

* 

0.000 


0.000 

* 

-0-000 

* 

-0.000 

* 

-0.000 

* 

-0.000 

* 

-0*000 

* 

-C.000 

* 

-0.000 

* 

-0.000 

* 

-0.000 

♦ 

-0 .0 00 

* 

-0.000 

* 

-0 -000 

* 

^0.000 


-0.000 

* 

-0.000 

* 

-0*000 

* 

-0.000 

* 

-0.000 

* 

-0 *000 

* 

-o.aoo 

* 

-0.000 

* 

-0*000 

* 

-a. oao 

* 

0.000 

* 

0.000 

* 

0*002 

* 

0*044 

* 

9-656 

* 



0*2416 

0*2444 

0.2446 

0*2444 

0*2442 

0*2439 

0*2436 

0*2432 

0.2427 

0*2421 

0*2414 

0*2406 

0*2396 

0.2363 

0.2366 

0.2346 

0*2320 

0.2207 

0.2240 

0.2204 

0.2104 

0*2102 

0.2040 

0.1994 

0*1941 


0*1799 
0.1 759 
0.1725 
0.1695 
0.1 671 
0*1652 
0.1 639 
0.1631 
0.1620 
0*1632 
0.1 641 
0.1656 
0.1677 
0*1 704 
0*1737 
0.1776 
0.1023 
0.1576 
n.l 935 
0.2002 
0.2076 
0.2157 
0.2245 
0.2337 


* 

* 

* 

* 

* 


* 

* 

* 

* 



* 


* 


♦ 


* 


CO 


* 


* 


* 


* 

* 


* 


* 


* 

* 

♦ 

* 

■*- 

* 

* 


4 


* 


* 


□ 



ittjwWVh i ii 
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GRAPHS FDR REFLECTION CYCLE NUMBER 
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ID C M A } TOTAL DIODE CURRENT VS TIME FOR ONE LQ CYCLE 


VOt VOLTS) 


DIODE VOLTAGE VS TIME FDR ONE LO CYCLE 


11*557 
11*143 
10*652 
10*069 
9.424 
S. 727 
7*964 
7 . ue 
6.181 
5 • 159 
A* 068 
2*923 
1.731 
0.501 
-0 * 754 
-2*013 
-3,246 
-4.427 
-S.53C 
-6*535 
-7.420 

-a. 161 

- 8 . 732 
-9*117 
-9.304 
-9.290 
-9*076 
-B.664 
-S 'i 056 
-7.265 
-6.306 
-5.206 
-3.990 
-2*681 
-1.301 
0*125 
1 .574 
3.007 
4.396 
5.722 
6.9 57 
8.117 
9.152 
10.051 
10.797 
11.385 

u.Bie 

12.101 

12.225 

12.179 

11.940 



0.346 
3.363 
0.365 
0*363 
0*362 
0.360 
0.358 
0.356 
0*353 
0*349 
0*345 
0.340 
0*333 
0*325 
0*314 
0.300 
0*252 
0*259 
0.230 
O* 194 
0.152 
0* 104 
0*051 
-O. 007 
-0.068 
-0. 132 
-0* 196 
-0.260 
-0*322 
-D.379 
-0.431 
-0*476 
- 0 .S 12 
-0.539 
-0*555 
-0-560 
— 0 * 553 
-0*534 
-0.505 
-0.465 
-0.416 
-0.359 
-0.294 
-0*225 
-0.151 
-0.075 
0.003 
0. 079 
O. 155 
0.227 
0.295 




* 


* 


* 



* 



* 


#*#**♦*** 
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RESULTS OF THE SMALL-SIGNAL ANALYSIS 


FCURIER COEFFICIENTS OF THE DIODE CONDUCTANCE AND CAPACITANCE ttAvEFDRUS 


CJMAG.GJPHA 

1 

7.1 3Sti3a2 

-30* 

2 

5.07CD-02 

— 72. 

3 

2.8S3D-02 

-95* 

4 

1*6100-02 

-92* 

r* 


5 

1 • 3570- C 2 

-se* 

6 

1 * 135D— 0 2 

-92. 

7 

8.810D-03 

-92. 

E 

7.359D-03 

-92. 

* ’ 3 

« 

C JMAG *CJPHA 

1 

2*1460-14 

—60 * 

c 

2-448D-15 

5. 

3 

1 .9200-15 

-12. 

4 

7*7720-16 

-1 7. 



5 

4.8880-16 

-£. 

6 

2*50 4D-1 e 

— 6. 

7 

2*330 0-16 

-10. 

e 

1 • 6Q3D— 16 

-e. 



CC TERKS: GJWAG = 7^976C-02 CJHAG = 6.CC2D-13 


CCN^ERS ION LOSS MATRIX { DB1 

4 2 2 1 0 —1 —2 —3 ->4 



4 


0 .0 

36.24 

36.71 

32*35 

32*05 

36*30 

44.04 

59*41 

67*42 


3 


39.37 

C .0 

22. £5 

32*77 

34.12 

37*00 

47.1? 

61 .69 

66*77 


2 


36.72 

20.33 

fl.C 

14.07 

27*47 

28*41 

36*23 

63.34 

53.21 


1 


35.10 

2C.77 

15. £1 

0.0 

5*87 

27*91 

27.28 

37.84 

42 • 7 € 


0 


34.56 

27.54 

22* 24 

6*51 

0*0 

5*95 

22*42 

36.88 

34.21 


-1 


42. 5 8 

37 ,44 

27. £? 

28.72 

£.47 

0*0 

14.60 

30 n 25 

34.26 


-2 


52 .42 

55.71 

36* 17 

20*03 

27.37 

13*01 

0.0 

29*24 

35.66 


-3 


66.72 

62.22 

4 £ * 12 

37.20 

33. B4 

32*43 

27.55 

0*0 

38. 7£ 


-4 


67.34 

£5.91 

44.72 

36.34 

31*90 

31.70 

35, 6 E 

37.55 

0.0 

LFPER 

SIDEBAND 

CONVERSION 

LCSS: L{ 0 t. 

1) = €.51 

05 






LCWER 

SIDEBAND 

CCNVERSICA 

Less: L < 0 , 

-1) “ £ .55 

D 8 






INPUT 

impedances 

4 

3 

2 

1 

0 

-1 

—2 

-3 

-4 

REALIZIN1 : 


9.14 

£.01 

10.16 

19*35 

94.14 

23.38 

10.97 

8 . 10 

9.26 

IPAG(ZIN) I 


-5.57 

-12 *SH 

-2C.££ 

-38*64 

-1 6.36 

44*21 

22*60 

14*43 

9-95 

IF output impedance: 

Z IFOUT = 94 

.14 + J -18 

.36 







ECUIVALENT INPUT NOISE TEMPERATURES! 7M THERM SHOT 

456.3 146.2 310.1 


■I 
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APPENDIX H 


Modifications to the Mixer Analysis Program for 
Running the Examples in Section 4 


The mixer analysis program appearing in Appendix I of this report must be 
altered slightly in order to use it for the examples described in section 4. The appro 
priate statement modifications for running each example are listed on the following 
pages. Line numbers indicate whether the statements are to be inserted between or 
used as replacements for those in the program of Appendix I. 
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A n. 1: Alterations to the program in Appendix I for studying the effects of series 

inductance on the performance of the simple mixer circuit of Pig. 9 using a con- 
ventional Schottky-barrier diode. To represent the circuit of Pig. 8 replace 
line 622 with DATA ZR/50.0D0, 7*0,. 01D0 , and substitute DATA Z0, ZQACC/ 
50. 0D0, 0. 01D0/ for line 1696. Por the Schottky diode whose junction capaci- 
tance is independent of voltage set GAM to zero in line 1674. 

Ls: THE DIODE SERIES INDUCTANCE WHICH WHEN COMBINED WITH ZR AND 190.1 

ZI PRODUCES THE EMBEDDING IMPEDANCES ZER AND ZEI. 100.2 

LSDAT 5 ARRAY HOLDING THE VALUES OF LS FOR A COMPLETE RUN OF THE 190.3 

PROGRAM. 190.4 

NRUN: COUNTS THE NUMBER OF COMPLETE RUNS OF THE PROGRAM. 209.1 

RESULT: ARRAY FOR STORING RESULTS OF COMPLETE RUNS OF THE PROGRAM. 221.1 

ZI: IMAGINARY PART OF THE IMPEDANCE OF THE DIODE MOUNT WHICH V/HEN 283.1 

COMBINED WITH LS PRODUCES THE EMBEDDING IMPEDANCE ZEI. 283.2 

ZR: REAL PART OF THE IMPEDANCE OF THE DIODE MOUNT WHICH WHEN 294.1 

COMBINED WITH LS PRODUCES THE EMBEDDING IMPEDANCE ZER. 294.2 

—THE MAIN DRIVER PROGRAM 351,01 

351.02 

THE VARIABLE TYPES USED IN THIS ROUTINE ARE AS FOLLOWS: 351.03 

FOR COMMON/CONST/: 351.04 

REAI.*8 QEL ,BOLTZ,PI ,TEQ 351.05 

C— FOR COMMON/DIODE/: 351.06 

REAL *8 ALP , PHI , GAM , C0, 1 S , RS , RSKIH , FP , IF , IGU , ICO ,G0 ,C<3 351.07 

C FOR COMMON/TLINE/: 351.08 

REAL*8 Z0, ZQACC 351.09 

INTEGER ZQFLAG 351.1 

C FOR COMMON/VOLTS/: 351.11 

REAL*8 VL0,VDSIAS,IDBIAS,AV!8),BV'P},VIDC 351.12 

C- — FOR VARIABLES NOT IN ANY COMMON BLOCKS: 351.13 

COMPLEX* 1 6 ZIFOUT 351.14 

REAL *8 LSDAT! 35 }, LS , RESULT! 20, 20) , LOPWR , SHOT , THERM ,TM , XLMAT( 9*9) 351.15 

INTEGER NRUN 351.16 

C— -THE COMMON BLOCKS USED ARE: 351.17 

COMMON/CONST/QEL ,BOLTZ,PI,TEQ 351.18 

COMMON /D I ODE /ALP , PH I , GAM, C0, 1 S , RS, RSKIN , FP , IF , IGJ , ICG , GO , CO 351.19 

COMMON/TLINE/Z0, ZQACC, ZQFLAG 351.2 

COMMON / VOL TS/AV ,BV,VIDC,VLO,VDBIAS, IDBIAS 351.21 

C— THE VALUES OF LS FOR THE VARIOUS PROGRAM RUNS 351.22 

DATA LSDAT/ 1.0D-11 ,1 . 5 D - 1 1 , 2 -0D-1 1 , 2 . 5D-11 ,3 .0D-11 , 3 .5D-11 , 351.23 

*4.00-11,4.50-11,5.00-11,5.50-11,6.00-11,6.50-11,7.00-11,7.50-11 351.24 

*,8.0D-11 ,8.5D-H‘,9.0D-11 ,9.50-11 ,10. 0D-11, 10.5D-I1 ,11 .00-11 , 351.25 

*12.0D-11,13.0D-1I ,14.00-11 , 1S.0D-11 , 16.0D-I1 , 17.0D-11 ,18.0D-11 , 351.26 

*19 .0D-1 1 , 20.0D-1 1 ,21 .0D-1 1 ,22.00-11,23.00-11 ,24.0D-11 ; 25.0D-11/ 351.27 

NRUN=0 351.28 

C RUN THE FROGRAM AS MANY TIMES AS THERE IS DATA 351.29 

DO 10 1=1,35 351.3 

C— SET LS TO ITS VALUE FOR THIS RUN 351.31 

LS=LSDAT( I ) 351.32 

CALL LGS IG( LS , LOPWR) 353. 

CALL SMS IG< TM , SHOT, XL MAT , ZI FOUT ) 355. 

C INCREMENT THE VARIABLE WHICH COUNTS THE NUMBER OF RUNS 355.01 

NRUN=NRUN+1 355.02 

C STORE THE RESULTS OF EACH RUN IN AN ARRAY 355.03 

RESULT! NRUN, 1 >=C0 355.04 

RESULT!NRUN,2)=RS 355.05 

RESULT! NRUN , 3 )=IS 355.06 

RESULT! NRUN,4)=FP 355.07 

RESULT! NRUN , 5 )=I F 355.08 

RESULT(NRUN,6 )=ALP 355.09 

RESULT! NRUN,7)=PHI 355.1 

•RESULT! NRUN, 8 )=GAM 355.11 

RESULT! NRUN, 9 >=I0BIAS 355.12 

RESULT! NRUN, 10)= VDB I AS 355.13 

RESULT! NRUN , 1 1 )=TEQ 355.14 
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RESULT(NRUN, 12 )=Z0 355.15 

RESULT! NRUN , 13 )=VLO 355.16 

RESULT! NRUN, 14 >=LOPWR 355.17 

RESULT! NRUN , 1 5 )=LS 355.18 

RESULT! NRUN, 16 )=DREAL( ZIFOUT ) 355.19 

RESULT! NRUN , 17 l-DIMAG! ZI FOUT ) 355.2 

RESULT! NRUN , 1 8 >=XLMAT( 5,4) 355 2 1 

RESULT! NRUN , 19 )=TM 355.22 

RESULT! NRUN , 20) “SHOT 355.23 

10 CONTINUE 355.24 

C— WRITE THE RESULTS OF ALL THE RUNS 355.25 

WR I T£ ( 6 , 1 50 ) 355.26 

150 FORMAT! IH1 , IX, 'RESULTS OF THE RUNS ON THIS PRINTOUT' //59X, ' DATA* / 355.27 

1/T5, ‘ C0 1 , T1 6 , ' RS ' ,T27, ’IS’ ,T38, ' FP ' , T49 , ' IF 1 , T59 , ’ALP' ,T70, 'PHI ’ , 355.28 

2T81 . ’GAM’ ,T91 , ’ IDBIAS' ,T102, ’ VDBIAS ’ ,T1 14 , ’TEQ’ ,T126, *Z0’ ) 355.29 

WRITE ! 6,1 80X RESULT! 1 ,J), 0=1,12) 355.3 

180 FORMAT ! / 12! 1PE 10. 3 , 1 X ) > 355.31 

WRITE! 5 , 1 90 ) 355.32 

190 FORMAT! ///58X, ' RESULTS ' //T6, ’VLO' ,T19, 'LOPWR' ,T32, 'LS' ,T41, 355.33 

1 'RE!ZIFOUT>' ,T54 , ’ IM< ZIFOUT ) ' , T69 , 'XL!0, 1 ) ’ ,T84, 'TM’.T9B, 355.34 

2' SHOT ' ,T109, 'THERM' ) 355.35 

DO 20 1 = 1, NRUN 355.36 

THERM=RESULT! 1,19 l-RESULT! 1,20) 355.37 

C-— WRITE THE RESULTS OF ALL THE RUNS ON A DISK FILE 355.38 

WRITE! 12 ) ! RESULT! I ,0 ) , J=1 ,20) 355.39 

WR I TE { 6 , 200 ) (RESULT! 1,0 ) , J= 13 ,20) .THERM 355.4 

200 FORMAT! /IX , 9! IPE 10. 3 , 3X ) > 355.41 

20 CONTINUE 355.42 

SUBROUTINE LGSIGUS, LOPWR) 361. 

REAL*8 LS 418.1 

CALL ZEMBED(ZER ,ZE I , ZEMBDC .ZEMBSB , LS , RS , RSLO , RSSB , RSKIN , FP , I F 437. 

1 , IDBIAS, RSSB, RSLO, LS.NH, NHP1 ,NHD2 ) 487. 

SUBROUTINE ZEMBED < ZER ,ZE 1 , ZEMBDC .ZEMBSB , LS , RS , RSLO , RSSB , RSKI N , 600. 

REALMS ZER! NH ) ,ZE I ( NH ) , ZEMBDC ,LS,RS,RSKIN,FP,IF , P I ,ZR(8) ,ZI ( 8 ) 618. 

C THE IMPEDANCES OF THE DIODE MOUNT ARE GIVEN: 621. 

DATA ZR/50.000, 7* 1000.0D0/ 622. 

DATA ZI/8 ,V 0.0D0/ 623. 

C FORM THE EMBEDDING IMPEDANCES AT EACH HARMONIC OF THE LO 623.1 

DO 10 1=1, NH 623.2 

ZER! I )=2R( I ) 623.3 

10 ZEI! I)=ZI< 1 )+2,0D0*PI*FP*DFLOAT( I )*LS 623.4 

C FORM THE SIDEBAND IMPEDANCES (SIDEBAND I = ELEMENT NH/2 +1-1) 623.5 

DO 30 1=1, NHP 1 623.6 * 

K = N H D 2 P 1 - 1 623.7 

IF ( K . EO.0) GOTO 20 623.8 

ZEMBSB! I )=DCMPLX(ZR( IABS(K) },DSIGN(ZI( IABS! K > ) , D FLOAT! K ) )+ 623.9 

1LS'- V 2 .0D0 ft PI*( DFLOAT! K >*FP+IF ) ) 624 . 

GOTO 30 624.1 

20 ZEMBSB! I ) = DCMPLX(ZEMBf)C,0.0D0> 624.2 

30 CONTINUE 624.3 

SUBROUTINE SMSIG! TM.SHOT.XLMAT.ZIFOUT) 969. 

I VDB IAS , IDBIAS, RSSB, RSLO, LS.NHARM, NHP 1 ,NHD2> 1305. 

REAL*8 ZER! NHARM > , ZE 1 ( NHARM > , RSSB! NHP1 ), RSLO! NHARM) ,LS 1329. 

WRITE! 6.175) LS 1371.1 

175 FORMAT! /IX,* ADDITIONAL INPUT DATA: ' ,T50, ’ LS ' /T46, 1PE10. 3 ) 1371.2 

DATA VI.OINC, IDCACC/0. 1D0.0.02D0/ 1669. 

U.'-.TA TEO/295.0D0/ 1672. 

DATA ALP , PH I , GAM/34 . 78D0,0. 90D0, 0 . 5D0/ 1674. 

DATA C0, IS.RS/l . 18D-I4, 1 . 4D-15 , 4 . 4D0/ 1675. 

DATA FP, IF/115. 0D9, 4. 0D9/ 1676. 

DATA R3KIN/0.0D0/ 1677. 

C*** DELETE LINES 1683 TO 1693 SINCE THE EMBEDDING IMPEDANCES ARE 1683. 

C"** FORMED IN SUBROUTINE 2EMBED. 1684. 

DATA NH, HLO.NPTS.NVLU, NITER, NPRINT/8, 1 ,51 ,50,500, 100/ 1692. 

DATA VDINI T, ACC, NDIM/0.0D0, 1.00-6,1/ 1694. 

DATA Z0.ZQACC/200.0D0, 0.0100/ 1696. 

DATA VDBIAS. I UB1AS/0. 400,0.00200/ 1698. 

DATA VLO/ I . 30D0/ 1699. 
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A H. 2: Alterations to the program in Appendix I for studying the effects of the series 

inductance on the performance of the simple mixer circuits of Figs. 8 and 9, using 
a Mott diode with the C-V relation given in Fig. 11. Only changes which are 
additional to those in AH. 1 are listed. 


B: THE V INTERCEPTS OF THE C-V PIECEWISE LINEAR APPROXIMATION 117.1 

USED FOR THE MOTT DIODE EXAMPLES. 117.2 

C: THE CAPACITANCE AT VOLTAGE Y(l) FOR THE MOTT DIODE. 121.1 

C(VD>: FUNCTION WHICH RETURNS A VALUE OF CAPACITANCE FOR A GIVEN 302 . 1 

VOLTAGE AT THE TERMINALS OF THE MOTT DIODE. 302.2 

REAL * 8 X,Y<1 >,DERV(1>,V,C 8154. 

CO=2.0D0’*PI*FP*C(DMIN1(Y< 1 > ,0. 999D0*PHI > ) 8114. 

REAL FUNCTION C*0<VD> 890.01 

C ' 898.02 

C C( VU > IS USED FOR THE MOTT DIODE EXAMPLES ONLY. 8S8.03 

C— -THIS FUNCTION ROUTINE RETURNS THE PROPER VALUE OF CAPACITANCE FOR 098.04 

C A GIVEN VALUE OF THE DIODE VOLTAGE AS CALCULATED IN DRKGS. 898.05 

C — IT USES DATA FOR THE SLOPE AND C INTERCEPT OF THE PIECEWISE 890.06 

C— -LINEAR APPROXIMATION TO THE C-V CURVE OF THE MOTT DIODE. THE DOPING 090.07 

C- — PROFILE OF THE MOTT DIODE V/AS KINDLY SUPPLIED DY M.V. SCHNEIDER, 890.00 

c — bell laboratory , holmdel » n . o . 090.09 

c — 090.1 

C— THE VARIABLE TYPES USED IN THIS SUBROUTINE ARE AS FOLLOWS: 890.11 

REAL *8 M( 17>,B< 17>,VLIM< 18) ,VD 898.12 

C— THE SLOPE AND INTERCEPT DATA 898.13 

DATA M/1 .20-16,6.0-16, 9.0-16,1. 1030-1 5 , 1 . 783D- 1 5 ,2 . 375D-15 , 3 . D-l 5 890.14 

1 ,3.9640-15,6.5150-15,9.20-15, 10. 4550- 15,1 5. 833D-1 5, 2 1.750-1 5 090.15 

2,37. 5D-1S ,74.0-15, 154, 5D-1 5, 4 132. 6570-15/ 890. 1C 

DATA R/9. 78D-I5.il. ID- 15, 11. 4D-15.il . 542D- 1 5 , 11 . 753D-15 , 1 1 . 8D-1 5 098.17 

1,11.00-15,1 1 .742D-15, 10. 84 ID-15, 9. 62CD- 15,9. 036D- 15,5.9170-1 5 890.10 

2, 2. 13D-15, -8. 500-15,-36. 32D- 1 5, -100. 720-15,-3442. 372D-15/ 898.19 

C— THE VOLTAGE LIMITS ON THE REGIONS OF CAPACITANCE VARIATION 890.2 

DATA VLIM/-4.D0,-2.D0,-1 .D0.-7.D-1 , - . 31 D0 , -8 . D-2 ,0 . D0 , 6 . D-2 , .3400 098.21 

1 , .47D0, .52D0, .5800, .G4D0, .68D0, .76D0, ,8D0, . 84D0, . 0991 00/ 898.22 

1=1 898.23 

IFIVD.LT.VLIMI 1 ) ) GOTO 15 098.24 

IFIVD.EQ.VLIMU )> GOTO 10 090.25 

IF<VD.GT.VL1M(18)> GOTO 20 098.26 

C— FIND THE REGION IN WHICH THE VOLTAGE FROM DRKGS FALLS AND CALC, C 090.27 
5 IF < VD . GT . VL I M < 1)1 GOTO 10 898.28 

C=M< I-I >*VD+B{ 1-1 ) 898.29 

GOTO 25 898.3 

10 1=1+1 898.31 

GOTO 5 890,32 

15 WRITE < 6 , 100 ) VD.VLIMU) 898.33 

100 FORMAT! ////IX , ' UNDERFLOW: VD “\F8.4,' IS LESS THAN '.F7.4// 1 898.34 

GOTO 25 898.35 

20 WRITE (6, 200) VD , VL IM( 18) 890.36 

200 FORMAT! ////IX, 'OVERFLOW: VD =',F8.5,' IS GREATER THAN \F8.5//> 898.37 

25 RETURN 098.38 

END 898.39 

890.4 

898.41 

890.42 


o n 


99 


A n. 3: Alterations to the program in Appendix I for studying the performance of the 

Scliottky diode mixers when the zero voltage junction capacitance CO is varied. \ 

One value of series inductance is used for each run. Only changes which are 
additional to those listed in A n. 1 are shown here. 


C0DAT: ARRAY HOLDING THE VALUES OF C0 FOR EACH RUN OF THE ENTIRE 131.1 

ANALYSIS. 131.2 

REAL*8 C0DAT! 20 ),LS, RESULT! 2.0, 20 ) .LOPWR, SHOT, THERM, TM.XLMAT! 9,9 ) 351.15 

C-— THE VALUES OF C0 FOR THE VARIOUS PROGRAM RUNS 351.22 

DATA C0DAT/1 .0D-15, ,2.00-15, , 3 .0D-1 5 , 4 .0D-1 5 , 5 .00- 1 5 , 6 .0D- 1 5 , 351.23 

*7.00-15,0.00-15,9.00-15, 10. 0D- 15, 11 .0D-1 5, 12.0D-15, 13.00-15, 351.24 

*14.0D-15, 15.0D-15, 16.00-15, 17.0D-15, 18.0D-15, 19.0D-15, 351.25 

» *20.00-15/ 351.26 

DATA LS/.04D-9/ 351.27 

DO 10 1=1,20 351.3 

C-— SET C0 TO ITS VALUE FOR THIS RUN 351,31 

C0=C0DAT( I ) 351.32 

RESULT! NRUN,1)*LS 355.04 

RESULT! NRUN, 15 )=C0 355.10 

1/T5 , 'LS',T16, * RS ' ,T27, ' I S 1 , T38 , ' FP * ,T49, ' I F ' , T59 , ' ALP ‘ , T70, 'PHI ' , 355.20 

190 FORMAT! ///58X, 1 RESULTS ' //T6 , 1 VLO ' , T1 9 , ' LOPWR 1 ,T32, ‘C0 1 ,T41 , 355.33 
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APPENDIX m 

Results of Program Runs to Study the Effects of Series Inductance and 
Junction Capacitance on the Performance of Some Simple 
Mixer Circuits Described in Section 4 


This appendix contains graphs of the equivalent input noise temperature, upper 
sideband conversion loss and IP output impedance as a function of series inductance LS 
and junction capacitance CO. 

The mixer analysis program, with the modifications in Appendix n, was used 
for each of three diodes in the mixer circuits of Figs. 8 (short circuited harmonics) 
and 9 (open circuited harmonics). The diodes were forward biased at 0.4 V in all 
cases, and the LO power adjusted to give a rectified dc current of 2. 0 mA. The 
signal, LO, and intermediate frequencies were 119 GHz, 115 GHz, and 4 GHz re- 
spectively. 


PRECEDING PAGE SIANK NOT FILMED 
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A m. 1: Mixer performance data plotted as a function of series inductance for three 

different diodes in the simple mixer circuits of (a) Fig. 8 and (b) Fig. 9. In 
all cases the diodes were forward biased at 0.4, and the LO level adjusted to 
give a rectified current of 2. 0 mA. The signal, LO and intermediate frequencies 
were 1X9 GHz, 115 GHz and 4 GHz respectively. 
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CO = 11.8 fFd. 

•f = Schottky diode with varying capacitance (GAM=0.5) 

O = Schottky diode with constant capacitance (GAM=0) 
x = Mott diode with a realistic C-V variation 
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HI. 2 Mixer perfoiunance data plotted as a function of zero voltage junction capaci- 
tance for a conventional Ga As Schottky diode (GAM = * 5) and a Ga As Schottky 
diode with no capacitance variation (GAM=0) in the mixer circuits of (a) Mg. 8 
and (b) Mg. 9. The diodes were forward biased to 0.4 V and the LO level ad- 
justed to give a rectified current of 2. 0 mA. The signal, LO and intermediate 
frequencies were 119 GHz, 115 GHz and 4 GHz respectively. Graphs appear for 
nine different values of series inductance LS. 
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(b) 

o/e harmonics 


LS = 0. 04 nH. 


+ = varying capacitance (GAM=0.5) 
0 ~ constant capacitance (GAM=0) 
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s/c harmonics 
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o/c harmonics 


LS = 0.06 nH. 
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(a) 

s/c harmonics 
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+ = varying capacitance (GAM=0.5) 
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(b) 

o/c harmonics 


(a) 

s/c harmonics 

LS = 0.2 nH. 

+ = varying capacitance (GAM=0.5) 
0 * constant capacitance (GAM=0) 


113 



' K 

•u 


■ ■ \- 





ACKNOWLEDGEMENT 


The authors wish to thank M. V* Schneider of Bell Telephone Laboratories, 
Hokndel, N. for supplying the data on the doping profile of his Mott diodes. 



* 


PRECEDING page blank not filmed 




i ij l j ta «g «gjlgflF ■ - 






115 


REFERENCES 

[1] Held, D. N. and Kerr, A. R. , ’’Conversion loss and noise of microwave and 
millimeter-wave mixers: Part 1 - Theory” and ’’Part 2 - Experiment”, IEEE 
Trans, on Microwave Theory and Tech. , vol. MTT-26, no. 2, pp. 49-61, Feb. 

1978. 

[2] Held, D. N. , "Analysis of room temperature millimeter- wave mixers using GaAs 
Schottky barrier diodes", ScD dissertation, Columbia Univ. , New York, N. Y. , 
1976. 

[3 ] Torrey, H. C. and Whitmer, C. A. , Crystal Rectifiers , MIT Radiation Lab. 
Series, vol. 15, New York, McGraw Hill, 1948. 

[4] Fieri, D. A. and Cohen, L. D. , "Nonlinear analysis of the Schottky- Barrier 
mixer diode", IEEE Trans, on Microwave Theory and Tech. , vol. MTT-21, 
pp. 39-43, Jan. 1973. 

[5] Egami, S. , "Nonlinear, ’inear analysis and computer aided design of resistive 
mixers", IEEE Trans, on Microwave Theory and Tech. , vol. MTT-22, pp. 270- 
275, March 1974. 

[6] Gwarek, W. K. , "Nonlinear analysis of microwave mixers", M.S. thesis, M.I.T., 
Cambridge, Mass. , 1974. 

[7] Kerr, A. R. , "A technique for determining the local oscillator waveforms in a 
microwave mixer", IEEE Trans, on Microwave Theory and Tech. , vol. MTT-23, 
pp. 828-831, Oct. 1975. 

[8] Strutt, M. J. O. , "Noise figure reduction in mtxer stages", Proc. IRE, vol. 34, 
no. 12, pp. 942-950, Dec. 1946. 

[9] van der Ziel, A., Noise: Sources, Characterization, and Measurement , 

Englewood Cliffs, N. J. , Prentice Hall, 1970. 


PRECEDING PAGE NOT 


116 


[10] Kim, C. S. , "Tunnel diode converter analysis", IRE Trans, on Electron 
Devices, vol. ED-8, no. 5, pp. 394-406, Sept. 1961. 

[11] Dragone, C. , "Analysis of thermal and shot noise in pumped resistive diodes", 
Bell System Tech. Journal, vol. 47, pp. 1883-1902, 1968. 

[12] Uhlir, A,, Jr., "Shot noise in p-n junction frequency converters", Bell Syst, 
Tech. J. , vol. 37, no. 4, pp. 951-988, July 1958. 

[13] Saleh, A. A. M. , Theory of Resistive Mixers , Cambridge, Mass., MIT Press, 
1971, 


